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SPECIFICATION 
D EVICE AND METHOD OF CONTROLLING EXHAUST GAS SENSOR 
TEMPERATURE, AND RECORDING MEDIUM FOR EXHAUST GA S SENSOR 
TEMPERATURE CONTROL PROG RAM 

Technical field: 

,The present invention relates to an apparatus 
for and a method of controlling the temperature of an ex- 
haust gas sensor disposed in the exhaust passage of an in- 
ternal combustion engine, and a recording medium storing a 
program for controlling the temperature of such an exhaust 
gas sensor. 
Background art : 

Exhaust gas sensors are often disposed in the 
exhaust passages of internal combustion engines for detect- 
ing a physical quantity as to an exhaust gas component 
state, such as an exhaust gas component concentration, for 
the purpose of controlling the operation of the internal 
combustion engine or monitoring the status of an exhaust 
gas purifying system. Specifically, the exhaust gas sensor 
is disposed at a certain location in the exhaust gas pas- 
sage and has an element sensitive to an exhaust gas compo- 
nent state to be detected, the element being positioned for 
contact with an exhaust gas flowing through the exhaust 
passage. For example, an air- fuel ratio sensor such as an 
0 2 sensor or the like is disposed as an exhaust gas sensor 



upstream or downstream of an exhaust gas purifying catalyst 
disposed in the exhaust passage for the purpose of control- 
ling the air-fuel ratio of the internal combustion engine 
in order to keep well the purifying ability of the cata- 
lyst . 

Some air-fuel ratio sensors have a built-in 
heater for heating the active element thereof for increas- 
ing the temperature of the element and activating the ele- 
ment to enable the element to perform its essential func- 
tions and also removing foreign matter deposited on the 
element. For example, an air-fuel ratio sensor such as an 
O2 sensor or the like usually has an electric heater for 
heating the active element thereof. After the internal 
combustion engine has started to operate, the electric 
heater is energized to increase the temperature of the ac- 
tive element of the 0 2 sensor to activate the active ele- 
ment and keep the active element active. 

As shown in FIG. 3 of the accompanying drawings, 
the 0 2 sensor produces an output voltage Vout which changes 
with a large gradient with respect to a change in the air- 
fuel ratio of an exhaust gas, i.e., which is highly sensi- 
tive to a change in the air-fuel ratio, only in a small 
range A (near a stoichiometric air-fuel ratio) of values of 
the air-fuel ratio that is represented by an oxygen concen- 
tration in the exhaust gas to which the active element is 
sensitive. A change in the output voltage Vout of the 0 2 
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sensor, i.e., a gradient of the output voltage Vout with 
respect to the air-fuel ratio, is smaller in air-fuel ratio 
ranges that are richer and leaner than the highly sensitive 
range A. The output characteristics of the 0 2 sensor, 
i.e., the gradient of the highly sensitive range A, etc., 
vary depending on the temperature of the active element. 
When the air-fuel ratio is to be controlled using the out- 
put voltage from the 0 2 sensor, therefore, it is desirable 
to keep the output characteristics of the 0 2 sensor in a 
desired range as much as possible and hence to keep the 
temperature of the active element of the 0 2 sensor in a de- 
sired temperature range as stably as possible for better 
air-fuel ratio control. 

Not only 0 2 sensors but also many exhaust gas 
sensors have their output characteristics affected by the 
temperature of the active element. If the internal combus- 
tion engine is to be controlled using the output signal 
from the 0 2 sensor, then it is preferable to keep the tem- 
perature of the active element of the exhaust gas sensor in 
a desired temperature range as stably as possible for bet- 
ter engine control. When the active element of the ex- 
haust gas sensor is heated to clean the active element, it 
is also preferable to keep the temperature of the active 
element of the exhaust gas sensor in a desired temperature 
range for a better cleaning effect. 
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As disclosed in Japanese laid-open patent publi- 
cation No. 2000-304721 by the applicant of the present ap- 
plication, it is known to estimate the temperature of the 
active element of an exhaust gas sensor (an air- fuel ratio 
sensor in the publication) and control the energization of 
a heater (an electric heater) based on the estimated tem- 
perature for thereby keep the temperature of the active 
element in a desired temperature range to obtain appropri- 
ate output characteristics from the exhaust gas sensor. 
According to the disclosed arrangement, the resistance of 
the heater is recognized from detected values of a current 
flowing through the heater and a voltage applied across the 
heater, and the temperature of the active element is esti- 
mated based on the detected resistance of the heater. 

According to the disclosure of the above publi- 
cation, however, since a duty cycle which determines the 
electric power to be supplied to the heater is uniquely de- 
termined by a table from an estimated value of the tempera- 
ture of the active element of the exhaust gas sensor, it is 
difficult to control the temperature of the active element 
of the exhaust gas sensor stably at a desired temperature 
because of a change in the temperature of the exhaust gas , 
an ambient air temperature, etc. Furthermore, if the tem- 
perature of the active element of the exhaust gas sensor is 
low, then in order to keep the duty cycle of the heater at 
a certain maximum level, the heater tends to consume more 



electric power than required and the actual temperatures of 
the heater and the active element are liable to be exces- 
sively high, causing damage to the heater and the active 
element - 

The present invention has been made in view of 
the above background. It is an object of the present in- 
vention to provide an apparatus for and a method of con- 
trolling the temperature of the active element of an ex- 
haust gas sensor stably at a desired temperature. Another 
object of the present invention is to provide a recording 
medium storing a temperature control program for control- 
ling the temperature of the active element of an exhaust 
gas sensor stably at a desired temperature. 

Disclosure of the invention: 

An apparatus for controlling the temperature of 
an exhaust gas sensor according to the present invention is 
an apparatus for controlling the temperature of an exhaust 
gas sensor disposed in an exhaust passage of an internal 
combustion engine and having an active element for contact- 
ing an exhaust gas flowing through the exhaust passage and 
a heater for heating the active element. To achieve the 
above object, according to a first aspect of the present 
invention, an apparatus for controlling the temperature of 
an exhaust gas sensor is characterized by comprising means 
for sequentially acquiring element temperature data repre- 



senting the temperature of the active element , means for 
sequentially acquiring heater temperature data representing 
the temperature of the heater, and heater control means for 
sequentially generating a control input which defines an 
amount of heat generating energy supplied to the heater so 
as to equalize the temperature of the active element repre- 
sented by the element temperature data to a predetermined 
target temperature, and controlling the heater depending on 
the control input, and characterized in that the control 
input generated by the heater control means includes at 
least an input component depending on the difference be- 
tween the temperature of the active element represented by 
the element temperature data and the target temperature and 
an input component depending on the temperature of the 
heater represented by the heater temperature data. 

Similarly, a method of controlling the tempera- 
ture of an exhaust gas sensor according to the present in- 
vention is a method of controlling the temperature of an 
exhaust gas sensor disposed in an exhaust passage of an in- 
ternal combustion engine and having an active element for 
contacting an exhaust gas flowing through the exhaust pas- 
sage and a heater for heating the active element. Accord- 
ing to the first aspect of the present invention, a method 
of controlling the temperature of an exhaust gas sensor is 
characterized by comprising the steps of sequentially ac- 
quiring element temperature data representing the tempera- 



ture of the active element and heater temperature data rep- 
resenting the temperature of the heater, sequentially gen- 
erating a control input which defines an amount of heat 
generating energy supplied to the heater so as to equalize 
the temperature of the active element represented by the 
element temperature data to a predetermined target tempera- 
ture, and controlling the heater depending on the control 
input, and characterized in that when the control input is 
generated, the control input is generated so as to include 
at least an input component depending on the difference be- 
tween the temperature of the active element represented by 
the element temperature data and the target temperature and 
an input component depending on the temperature of the 
heater represented by the heater temperature data. 

Furthermore, a recording medium storing a tem- 
perature control program for an exhaust gas sensor accord- 
ing to the present invention is a recording medium readable 
by a computer and storing a temperature control program for 
enabling the computer to perform a process of controlling 
the temperature o an active element of an exhaust gas sen- 
sor disposed in an exhaust passage of an internal combus- 
tion engine and having the active element for contacting an 
exhaust gas flowing through the exhaust passage and a 
heater for heating the active element. According to the 
first aspect of the present invention, a recording medium 
storing a temperature control program for an exhaust gas 



sensor is characterized in that the temperature control 
program includes a program for enabling the computer to 
perform a process of sequentially acquiring element tem- 
perature data representing the temperature of the active 
element and heater temperature data representing the 
temperature of the heater, a control input generating 
program for enabling the computer to perform a process of 
sequentially generating a control input which defines an 
amount of heat generating energy supplied to the heater so 
as to equalize the temperature of the active element 
represented by the element temperature data to a 
predetermined target temperature, and a program for 
enabling the computer to perform a process of controlling 
the heater depending on the control input, wherein the 
control input generating program has an algorithm for 
enabling the computer to generate the control input so as 
to include at least an input component depending on the 
difference between the temperature of the active element 
represented by the element temperature data and the target 
temperature and an input component depending on the 
temperature of the heater represented by the heater tem- 
perature daAacording to the first aspect of the present in- 
vention, since the heater is controlled based on a control 
input (a manipulated variable for an object to be con- 
trolled) including an input component depending on the dif- 
ference between the temperature of the active element rep- 
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resented by the element temperature data and the target 
temperature and an input component depending on the tem- 
perature of the heater represented by the heater tempera- 
ture data, or stated otherwise, a control input generated 
by combining at least the above input components, when the 
temperature of the active element varies with respect to 
the target temperature, it is possible to converge the tem- 
perature of the active element to the target temperature 
while suppressing excessive variations of the control input 
which defines an amount of heat generating energy supplied 
to the heater. As a result, the temperature of the active 
element of the exhaust gas sensor can stably be controlled 
at the target temperature. 

According to the first aspect of the present in- 
vention, the element temperature data representative of the 
temperature of the active element may directly be detected 
and acquired using a temperature sensor, or may be esti- 
mated and acquired based on a suitable parameter or a 
model. The heater temperature data representative of the 
temperature of the heater may also be acquired in the same 
manner. 

The input component depending on the difference 
between the temperature of the active element and the tar- 
get temperature in the control input may be a component 
proportional to the difference, a component proportional to 
an accumulative sum (integrated value) of values of the 



difference at respective times, or a sum of these compo- 
nents. This also applies to a second aspect, a third as- 
pect, and a fourth aspect of the present invention which 
will be described later. 

According to a second aspect of the present in- 
vention, an apparatus for controlling the temperature of an 
exhaust gas sensor is characterized by comprising means for 
sequentially acquiring element temperature data represent- 
ing the temperature of the active element, means for se- 
quentially acquiring exhaust gas temperature data repre- 
senting the temperature of the exhaust gas, and heater con- 
trol means for sequentially generating a control input 
which defines an amount of heat generating energy supplied 
to the heater so as to equalize the temperature of the ac- 
tive element represented by the element temperature data to 
a predetermined target temperature, and controlling the 
heater depending on the control input, and characterized in 
that the control input generated by the heater control 
means includes at least an input component depending on the 
difference between the temperature of the active element 
represented by the element temperature data and the target 
temperature and an input component depending on the tem- 
perature of the exhaust gas represented by the exhaust gas 
temperature data. 

Similarly, according to the second aspect of the 
present invention, a method of controlling the temperature 



of an exhaust gas sensor is characterized by comprising the 
steps of sequentially acquiring element temperature data 
representing the temperature of the active element and ex- 
haust gas temperature data representing the temperature of 
the exhaust gas, sequentially generating a control input 
which defines an amount of heat generating energy supplied 
to the heater so as to equalize the temperature of the ac- 
tive element represented by the element temperature data to 
a predetermined target temperature, and controlling the 
heater depending on the control input, and characterized in 
that when the control input is generated, the control input 
is generated so as to include at least an input component 
depending on the difference between the temperature of the 
active element represented by the element temperature data 
and the target temperature and an input component depending 
on the temperature of the exhaust gas represented by the 
exhaust gas temperature data. 

Furthermore, according to the second aspect of 
the present invention, a recording medium storing a tem- 
perature control program for an exhaust gas sensor is char- 
acterized in that the temperature control program includes 
a program for enabling the computer to perform a process of 
sequentially acquiring element temperature data represent- 
ing the temperature of the active element and exhaust gas 
temperature data representing the temperature of the ex- 
haust gas, a control input generating program for enabling 
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the computer to perform a process of sequentially generat- 
ing a control input which defines an amount of heat gener- 
ating energy supplied to the heater so as to equalize the 
temperature of the active element represented by the ele- 
ment temperature data to a predetermined target tempera- 
ture, and a program for enabling the computer to perform a 
process of controlling the heater depending on the control 
input , wherein the control input generating program has an 
algorithm for enabling the computer to generate the control 
input so as to include at least an input component depend- 
ing on the difference between the temperature of the active 
element represented by the element temperature data and the 
target temperature and an input component depending on the 
temperature of the exhaust gas represented by the exhaust 
gas temperature data. 

According to the second aspect of the present 
invention, since the heater is controlled based on a con- 
trol input including an input component depending on the 
difference between the temperature of the active element 
represented by the element temperature data and the target 
temperature and an input component depending on the tem- 
perature of the exhaust gas represented by the exhaust gas 
temperature data, or stated otherwise, a control input gen- 
erated by combining at least the above input components, it 
is possible to converge the temperature of the active ele- 
ment to the target temperature while compensating for 
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variations of the temperature of the exhaust gas with re- 
spect to the temperature of the active element. Stated 
otherwise, it is possible to converge the temperature of 
the active element to the target temperature while sup- 
pressing variations of the temperature of the active ele- 
ment due to variations of the temperature of the exhaust 
gas. As a result, the temperature of the active element of 
the exhaust gas sensor can stably be controlled at the tar- 
get temperature (desired temperature) . 

According to the second aspect of the present 
invention, the element temperature data representative of 
the temperature of the active element may directly be de- 
tected and acquired using a temperature sensor, or may be 
estimated and acquired based on a suitable parameter or a 
model, as with the first aspect. The exhaust gas tempera- 
ture data representative of the temperature of the exhaust 
gas may also be acquired in the same manner. 

According to a third aspect of the present in- 
vention, an apparatus for controlling the temperature of an 
exhaust gas sensor is characterized by comprising means for 
sequentially acquiring element temperature data represent- 
ing the temperature of the active element, and heater con- 
trol means for sequentially generating a control input 
which defines an amount of heat generating energy supplied 
to the heater so as to equalize the temperature of the ac- 
tive element represented by the element temperature data to 



a predetermined target temperature, and controlling the 
heater depending on the control input, and characterized in 
that the control input generated by the heater control 
means includes at least an input component depending on the 
difference between the temperature of the active element 
represented by the element temperature data and the target 
temperature and an input component depending on the target 
temperature . 

Similarly, according to the third aspect of the 
present invention, a method of controlling the temperature 
of an exhaust gas sensor is characterized by comprising the 
steps of sequentially acquiring element temperature data 
representing the temperature of the active element, sequen- 
tially generating a control input which defines an amount 
of heat generating energy supplied to the heater so as to 
equalize the temperature of the active element represented 
by the element temperature data to a predetermined target 
temperature, and controlling the heater depending on the 
control input, and characterized in that when the control 
input is generated, the control input is generated so as to 
further include at least an input component depending on 
the difference between the temperature of the active ele- 
ment represented by the element temperature data and the 
target temperature and an input component depending on the 
target temperature . 



Furthermore, according to the third aspect of 
the present invention, a recording medium storing a tem- 
perature control program for an exhaust gas sensor is char- 
acterized in that the temperature control program includes 
a program for enabling the computer to perform a process of 
sequentially acquiring element temperature data represent- 
ing the temperature of the active element, a control input 
generating program for enabling the computer to perform a 
process of sequentially generating a control input which 
defines an amount of heat generating energy supplied to the 
heater so as to equalize the temperature of the active ele- 
ment represented by the element temperature data to a pre- 
determined target temperature, and a program for enabling 
the computer to perform a process of controlling the heater 
depending on the control input, wherein the control input 
generating program has an algorithm for enabling the com- 
puter to generate the control input so as to include at 
least an input component depending on the difference be- 
tween the temperature of the active element represented by 
the element temperature data and the target temperature and 
an input component depending on the target temperature. 

According to the third aspect of the present in- 
vention, since the heater is controlled based on a control 
input including an input component depending on the differ- 
ence between the temperature of the active element repre- 
sented by the element temperature data and the target tern- 
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perature and an input component depending on the target 
temperature, or stated otherwise, a control input generated 
by combining at least the above input components, it is 
possible to converge the temperature of the active element 
quickly to the target temperature. As a result, the tem- 
perature of the active element of the exhaust gas sensor 
can stably be controlled at the target temperature. 

According to the third aspect of the present in- 
vention, the element temperature data representative of the 
temperature of the active element may directly be detected 
and acquired using a temperature sensor, or may be esti- 
mated and acquired based on a suitable parameter or a 
model, as with the first aspect. 

According to a fourth aspect of the present in- 
vention, an apparatus for controlling the temperature of an 
exhaust gas sensor is characterized by comprising means for 
sequentially acquiring element temperature data represent- 
ing the temperature of the active element, and heater con- 
trol means for sequentially generating a control input 
which defines an amount of heat generating energy supplied 
to the heater so as to equalize the temperature of the ac- 
tive element represented by the element temperature data to 
a predetermined target temperature, and controlling the 
heater depending on the control input, and characterized in 
that the control input generated by the heater control 
means includes at least an input component depending on the 
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difference between the temperature of the active element 
represented by the element temperature data and the target 
temperature and an input component depending on the tem- 
perature of the active element. 

Similarly, according to the fourth aspect of the 
present invention, a method of controlling the temperature 
of an exhaust gas sensor is characterized by comprising the 
steps of sequentially acquiring element temperature data 
representing the temperature of the active element, sequen- 
tially generating a control input which defines an amount 
of heat generating energy supplied to the heater so as to 
equalize the temperature of the active element represented 
by the element temperature data to a predetermined target 
temperature, and controlling the heater depending on the 
control input, and characterized in that when the control 
input is generated, the control input is generated so as to 
include at least an input component depending on the dif- 
ference between the temperature of the active element rep- 
resented by the element temperature data and the target 
temperature and an input component depending on the tem- 
perature of the active element. 

Furthermore, according to the fourth aspect of 
the present invention, a recording medium storing a tem- 
perature control program for an exhaust gas sensor is char- 
acterized in that the temperature control program includes 
a program for enabling the computer to perform a process of 



sequentially acquiring element temperature data represent- 
ing the temperature of the active element, a control input 
generating program for enabling the computer to perform a 
process of sequentially generating a control input which 
defines an amount of heat generating energy supplied to the 
heater so as to equalize the temperature of the active ele- 
ment represented by the element temperature data to a pre- 
determined target temperature, and a program for enabling 
the computer to perform a process of controlling the heater 
depending on the control input, wherein the control input 
generating program has an algorithm for enabling the com- 
puter to generate the control input so as to include at 
least an input component depending on the difference be- 
tween the temperature of the active element represented by 
the element temperature data and the target temperature and 
an input component depending on the temperature of the ac- 
tive element. 

According to the fourth aspect of the present 
invention, since the heater is controlled based on a con- 
trol input including an input component depending on the 
difference between the temperature of the active element 
represented by the element temperature data and the target 
temperature and an input component depending on the target 
temperature for the active element, or stated otherwise, a 
control input generated by combining at least the above in- 
put components, it is possible to reduce overshooting of 



the temperature of the active element with respect to the 
target temperature, and to allow the temperature of the ac- 
tive element to smoothly follow the target temperature. As 
a result, the temperature of the active element of the ex- 
haust gas sensor can stably be controlled at the target 
temperature. 

According to the fourth aspect of the present 
invention, the element temperature data representative of 
the temperature of the active element may directly be de- 
tected and acquired using a temperature sensor, or may be 
estimated and acquired based on a suitable parameter or a 
model, as with the first aspect. 

In the first through fourth aspects of the in- 
vention, the temperature of the active element is directly 
controlled at the target temperature. Generally, the tem- 
perature of the heater and the temperature of the active 
element are highly correlated to each other in a steady 
state where those temperatures are substantially constant. 
Therefore, when the temperature of the heater is controlled 
at a certain target temperature, the temperature of the ac- 
tive element can indirectly be controlled at a temperature 
corresponding to the target temperature for the heater. 

In view of the foregoing, according to a fifth 
aspect of the present invention, an apparatus for control- 
ling the temperature of an exhaust gas sensor is character- 
ized by comprising means for sequentially acquiring element 
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temperature data representing the temperature of the active 
element, means for sequentially acquiring heater tempera- 
ture data representing the temperature of the heater, and 
heater control means for sequentially generating a control 
input which defines an amount of heat generating energy 
supplied to the heater so as to equalize the temperature of 
the heater represented by the heater temperature data to a 
predetermined target temperature, and controlling the 
heater depending on the control input, and characterized in 
that the control input generated by the heater control 
means includes at least an input component depending on the 
difference between the temperature of the heater repre- 
sented by the heater temperature data and the target tem- 
perature, and an input component depending on the tempera- 
ture of the active element represented by the element tem- 
perature data. 

Similarly, according to the fifth aspect of the 
present invention, a method of controlling the temperature 
of an exhaust gas sensor is characterized by comprising the 
steps of sequentially acquiring element temperature data 
representing the temperature of the active element and 
heater temperature data representing the temperature of the 
heater, sequentially generating a control input which de- 
fines an amount of heat generating energy supplied to the 
heater so as to equalize the temperature of the heater rep- 
resented by the heater temperature data to a predetermined 
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target temperature, and controlling the heater depending on 
the control input, and characterized in that when the con- 
trol input is generated, the control input is generated so 
as to include at least an input component depending on the 
difference between the temperature of the heater repre- 
sented by the heater temperature data and the target tem- 
perature, and an input component depending on the tempera- 
ture of the active element represented by the element tem- 
perature data. 

Furthermore, according to the fifth aspect of 
the present invention, a recording medium storing a tem- 
perature control program for an exhaust gas sensor is char- 
acterized in that the temperature control program includes 
a program for enabling the computer to perform a process of 
sequentially acquiring element temperature data represent- 
ing the temperature of the active element and heater tem- 
perature data representing the temperature of the heater, a 
control input generating program for enabling the computer 
to perform a process of sequentially generating a control 
input which defines an amount of heat generating energy 
supplied to the heater so as to equalize the temperature of 
the heater represented by the heater temperature data to a 
predetermined target temperature, and a program for ena- 
bling the computer to perform a process of controlling the 
heater depending on the control input, wherein the control 
input generating program has an algorithm for enabling the 
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computer to generate the control input so as to include at 
least an input component depending on the difference be- 
tween the temperature of the heater represented by the 
heater temperature data and the target temperature, and an 
input component depending on the temperature of the active 
element represented by the element temperature data. 

The fifth aspect of the invention corresponds to 
the first aspect of the invention, and provides advantages 
similar to those of the first aspect in relation to con- 
trolling the temperature of the heater at the target tem- 
perature. Specifically, according to the fifth aspect of 
the invention, since the heater is controlled based on a 
control input (a manipulated variable for an object to be 
controlled) including an input component depending on the 
difference between the temperature of the heater repre- 
sented by the heater temperature data and the target tem- 
perature (a target temperature for the heater) and an input 
component depending on the temperature of the active ele- 
ment represented by the element temperature data, or stated 
otherwise, a control input generated by combining at least 
the above input components, when the temperature of the 
heater varies with respect to the target temperature, it is 
possible to converge the temperature of the heater to the 
target temperature while suppressing excessive variations 
of the control input. As a result, the temperature of the 
heater of the exhaust gas sensor can stably be controlled 
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at the target temperature, and hence the temperature of the 
active element can stably be controlled at a temperature 
corresponding to the target temperature for the heater. 

According to the fifth aspect of the present in- 
vention, the element temperature data representative of the 
temperature of the active element and the heater tempera- 
ture data representative of the temperature of the heater 
may be detected and acquired by a temperature sensor, or 
may be estimated and acquired, as with the first aspect. 

The input component depending on the difference 
between the temperature of the heater and the target tem- 
perature in the control input may be a component propor- 
tional to the difference, a component proportional to an 
accumulative sum (integrated value) of values of the dif- 
ference at respective times, or a sum of these components. 
This also applies to a sixth aspect, a seventh aspect, and 
an eighth aspect of the present invention which will be de- 
scribed later. 

According to a sixth aspect of the present in- 
vention, an apparatus for controlling the temperature of an 
exhaust gas sensor is characterized by comprising means for 
sequentially acquiring heater temperature data representing 
the temperature of the heater, means for sequentially ac- 
quiring exhaust gas temperature data representing the tem- 
perature of the exhaust gas, and heater control means for 
sequentially generating a control input which defines an 
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amount of heat generating energy supplied to the heater so 
as to equalize the temperature of the heater represented by 
the heater temperature data to a predetermined target tem- 
perature, and controlling the heater depending on the con- 
trol input, and characterized in that the control input 
generated by the heater control means includes at least an 
input component depending on the difference between the 
temperature of the heater represented by the heater tem- 
perature data and the target temperature, and an input com- 
ponent depending on the temperature of the temperature of 
the exhaust gas represented by the exhaust gas temperature 
data . 

Similarly, according to the sixth aspect of the 
present invention, a method of controlling the temperature 
of an exhaust gas sensor is characterized by comprising the 
steps of sequentially acquiring heater temperature data 
representing the temperature of the heater and exhaust gas 
temperature data representing the temperature of the ex- 
haust gas, sequentially generating a control input which 
defines an amount of heat generating energy supplied to the 
heater so as to equalize the temperature of the heater rep- 
resented by the heater temperature data to a predetermined 
target temperature, and controlling the heater depending on 
the control input, and characterized in that when the con- 
trol input is generated, the control input is generated so 
as to include at least an input component depending on the 
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difference between the temperature of the heater repre- 
sented by the heater temperature data and the target tem- 
perature, and an input component depending on the tempera- 
ture of the temperature of the exhaust gas represented by 
the exhaust gas temperature data. 

Furthermore, according to the sixth aspect of 
the present invention, a recording medium storing a tem- 
perature control program for an exhaust gas sensor is char- 
acterized in that the temperature control program includes 
a program for enabling the computer to perform a process of 
sequentially acquiring heater temperature data representing 
the temperature of the heater and exhaust gas temperature 
data representing the temperature of the exhaust gas, a 
control input generating program for enabling the computer 
to perform a process of sequentially generating a control 
input which defines an amount of heat generating energy 
supplied to the heater so as to equalize the temperature of 
the heater represented by the heater temperature data to a 
predetermined target temperature, and a program for ena- 
bling the computer to perform a process of controlling the 
heater depending on the control input, wherein the control 
input generating program has an algorithm for enabling the 
computer to generate the control input so as to include at 
least an input component depending on the difference be- 
tween the temperature of the heater represented by the 
heater temperature data and the target temperature and an 
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input component depending on the temperature of the exhaust 
gas represented by the exhaust gas temperature data. 

The sixth aspect of the invention corresponds to 
the second aspect of the invention, and provides advantages 
similar to those of the second aspect in relation to con- 
trolling the temperature of the heater at the target tem- 
perature. Specifically, according to the sixth aspect of 
the invention, since the heater is controlled based on a 
control input including an input component depending on the 
difference between the temperature of the heater repre- 
sented by the heater temperature data and the target tem- 
perature and an input component depending on the tempera- 
ture of the exhaust gas represented by the exhaust gas tem- 
perature data, or stated otherwise, a control input gener- 
ated by combining at least the above input components, it 
is possible to converge the temperature of the heater to 
the target temperature while compensating for variations of 
the temperature of the exhaust gas with respect to the tem- 
perature of the active element and the temperature of the 
heater. Stated otherwise, it is possible to control the 
temperature of the heater to the target temperature while 
suppressing variations of the temperature of the heater due 
to variations of the temperature of the exhaust gas. As a 
result, the temperature of the heater of the exhaust gas 
sensor can stably be controlled at the target temperature, 
and the temperature of the active element of the exhaust 
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gas sensor can stably be controlled at a temperature corre- 
sponding to the target temperature for the heater. 

According to the sixth aspect of the present in- 
vention, the heater temperature data representative of the 
temperature of the heater and the exhaust gas temperature 
data representative of the temperature of the exhaust gas 
may be detected and acquired or estimated and acquired, as 
with the first aspect. 

According to a seventh aspect of the present in- 
vention, an apparatus for controlling the temperature of an 
exhaust gas sensor is characterized by comprising means for 
sequentially acquiring heater temperature data representing 
the temperature of the heater, and heater control means for 
sequentially generating a control input which defines an 
amount of heat generating energy supplied to the heater so 
as to equalize the temperature of the heater represented by 
the heater temperature data to a predetermined target tem- 
perature, and controlling the heater depending on the con- 
trol input, and characterized in that the control input 
generated by the heater control means includes at least an 
input component depending on the difference between the 
temperature of the heater represented by the heater tem- 
perature data and the target temperature and an input com- 
ponent depending on the target temperature. 

Similarly, according to the seventh aspect of 
the present invention, a method of controlling the tempera- 
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ture of an exhaust gas sensor is characterized by compris- 
ing the steps of sequentially acquiring heater temperature 
data representing the temperature of the heater, sequen- 
tially generating a control input which defines an amount 
of heat generating energy supplied to the heater so as to 
equalize the temperature of the heater represented by the 
heater temperature data to a predetermined target tempera- 
ture, and controlling the heater depending on the control 
input, and characterized in that when the control input is 
generated, the control input is generated so as to include 
at least an input component depending on the difference be- 
tween the temperature of the heater represented by the 
heater temperature data and the target temperature and an 
input component depending on the target temperature. 

Furthermore, according to the seventh aspect of 
the present invention, a recording medium storing a tem- 
perature control program for an exhaust gas sensor is char- 
acterized in that the temperature control program includes 
a program for enabling the computer to perform a process of 
sequentially acquiring heater temperature data representing 
the temperature of the heater, a control input generating 
program for enabling the computer to perform a process of 
sequentially generating a control input which defines an 
amount of heat generating energy supplied to the heater so 
as to equalize the temperature of the heater represented by 
the heater temperature data to a predetermined target tern- 



- 29 - 



perature, and a program for enabling the computer to per- 
form a process of controlling the heater depending on the 
control input, wherein the control input generating program 
has an algorithm for enabling the computer to generate the 
control input so as to include at least an input component 
depending on the difference between the temperature of the 
heater represented by the heater temperature data and the 
target temperature and an input component depending on the 
target temperature . 

The seventh aspect of the invention corresponds 
to the third aspect of the invention, and provides advan- 
tages similar to those of the third aspect in relation to 
controlling the temperature of the heater at the target 
temperature. Specifically, according to the seventh aspect 
of the invention, since the heater is controlled based on a 
control input including an input component depending on the 
difference between the temperature of the heater repre- 
sented by the heater temperature data and the target tem- 
perature and an input component depending on the target 
temperature, or stated otherwise, a control input generated 
by combining at least the above input components, it is 
possible to converge the temperature of the heater quickly 
to the target temperature. As a result, the temperature of 
the heater of the exhaust gas sensor can stably be con- 
trolled at the target temperature and the temperature of 
the active element can stably be controlled at a tempera- 
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ture corresponding to the target temperature for the 
heater. 

According to the seventh aspect of the present 
invention, the heater temperature data representative of 
the temperature of the heater may be detected by a tempera- 
ture sensor and acquired or may be estimated and acquired, 
as with the third aspect. 

According to an eighth aspect of the present in- 
vention, an apparatus for controlling the temperature of an 
exhaust gas sensor is characterized by comprising means for 
sequentially acquiring heater temperature data representing 
the temperature of the heater, and heater control means for 
sequentially generating a control input which defines an 
amount of heat generating energy supplied to the heater so 
as to equalize the temperature of the heater represented by 
the heater temperature data to a predetermined target tem- 
perature, and controlling the heater depending on the con- 
trol input, and characterized in that the control input 
generated by the heater control means includes at least an 
input component depending on the difference between the 
temperature of the heater represented by the heater tem- 
perature data and the target temperature and an input com- 
ponent depending on the temperature of the heater. 

Similarly, according to the eighth aspect of the 
present invention, a method of controlling the temperature 
of an exhaust gas sensor is characterized by comprising the 
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steps of sequentially acquiring heater temperature data 
representing the temperature of the heater, sequentially 
generating a control input which defines an amount of heat 
generating energy supplied to the heater so as to equalize 
the temperature of the heater represented by the heater 
temperature data to a predetermined target temperature, and 
controlling the heater depending on the control input, and 
characterized in that when the control input is generated, 
the control input is generated so as to Include at least an 
input component depending on the difference between the 
temperature of the heater represented by the heater tem- 
perature data and the target temperature and an input com- 
ponent depending on the temperature of the heater. 

Furthermore, according to the eighth aspect of 
the present invention, a recording medium storing a tem- 
perature control program for an exhaust gas sensor is char- 
acterized in that the temperature control program includes 
a program for enabling the computer to perform a process of 
sequentially acquiring heater temperature data representing 
the temperature of the heater, a control input generating 
program for enabling the computer to perform a process of 
sequentially generating a control input which defines an 
amount of heat generating energy supplied to the heater so 
as to equalize the temperature of the heater represented by 
the heater temperature data to a predetermined target tem- 
perature, and a program for enabling the computer to per- 
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form a process of controlling the heater depending on the 
control input, wherein the control input generating program 
has an algorithm, for enabling the computer to generate the 
control input so as to include at least an input component 
depending on the difference between the temperature of the 
heater represented by the heater temperature data and the 
target temperature and an input component depending on the 
temperature of the heater. 

The eighth aspect of the invention corresponds 
to the fourth aspect of the invention, and provides advan- 
tages similar to those of the fourth aspect in relation to 
controlling the temperature of the heater at the target 
temperature. Specifically, according to the eighth aspect 
of the invention, since the heater is controlled based on a 
control input including an input component depending on the 
difference between the temperature of the heater repre- 
sented by the heater temperature data and the target tem- 
perature and an input component depending on the tempera- 
ture of the heater, or stated otherwise, a control input 
generated by combining at least the above input components, 
it is possible to reduce overshooting of the temperature of 
the heater with respect to the target temperature, and to 
allow the temperature of the heater to smoothly follow the 
target temperature. As a result, the temperature of the 
heater of the exhaust gas sensor can stably be controlled 
at the target temperature and the temperature of the active 
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element can stably be controlled at a temperature corre- 
sponding to the target temperature for the heater. 

According to the eighth aspect of the present 
invention, the heater temperature data representative of 
the temperature of the heater may directly be detected and 
acquired using a temperature sensor, or may be estimated 
and acquired based on a suitable parameter or a model, as 
with the fourth aspect. 

It is preferable to combine two or more of the 
first through fourth aspects of the present invention with 
respect to any of the temperature control apparatus, the 
temperature control method, and the recording medium. If 
the temperature control apparatus according to the first 
and second aspects are combined, then the temperature con- 
trol apparatus according to the first aspect is character- 
ized by means for sequentially acquiring exhaust gas tem- 
perature data representing the temperature of the exhaust 
gas, wherein the control input generated by the heater con- 
trol means includes an input component depending on the 
temperature of the exhaust gas represented by the exhaust 
gas temperature data. Similarly, if the temperature con- 
trol methods according to the first and second aspects are 
combined, then the temperature control method according to 
the first aspect is characterized by the step of sequen- 
tially acquiring exhaust gas temperature data representing 
the temperature of the exhaust gas, wherein when the con- 
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trol input is generated, the control input is generated so 
as to further include an input component depending on the 
temperature of the exhaust gas represented by the exhaust 
gas temperature data. Furthermore, if the recording medi- 
ums according to the first and second aspects are combined, 
then the recording medium according to the first aspect is 
characterized in that the temperature control program fur- 
ther includes a program for enabling the computer to per- 
form a process of sequentially acquiring exhaust gas tem- 
perature data representing the temperature of the exhaust 
gas, wherein the control input generating program has an 
algorithm for enabling the computer to generate the control 
input so as to further include an input component depending 
on the temperature of the exhaust gas represented by the 
exhaust gas temperature data. 

By thus combining the first and second aspects, 
the advantages of those aspects are added together for con- 
trolling the temperature of the active element more stably 
at the target temperature. 

In the temperature control apparatus according 
to the first aspect or the combination thereof with the 
temperature control apparatus according to the second as- 
pect, the control input generated by the heater control 
means may include an input component depending on the tar- 
get temperature, providing an invention based on the combi- 
nation of the temperature control apparatus according to 
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the first and third aspects or an invention based on the 
combination of the temperature control apparatus according 
to the first through third aspects. Similarly, in the tem- 
perature control method according to the first aspect or 
the combination thereof with the temperature control method 
according to the second aspect, when the control input is 
generated, the control input may be generated so as to fur- 
ther include an input component depending on the target 
temperature, providing an invention based on the combina- 
tion of the temperature control methods according to the 
first and third aspects or an invention based on the combi- 
nation of the temperature control methods according to the 
first through third aspects. Furthermore, in the recording 
medium according to the first aspect or the combination 
thereof with the recording medium according to the second 
aspect , the control input generating program may have an 
algorithm for enabling the computer to generate the control 
input so as to further include an input component depending 
on the target temperature, providing an invention based on 
the combination of the recording mediums according to the 
first and third aspects or an invention based on the combi- 
nation of the recording mediums according to the first 
through third aspects. 

If the first and third aspects are thus com- 
bined, the advantages of those aspects are added together 
for controlling the temperature of the active element more 
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stably at the target temperature. Particularly, if the 
first through third aspects are combined, the advantages of 
those aspects are added together for appropriately increas- 
ing the stability of the temperature of the active element 
with respect to the target temperature. 

Moreover, in the temperature control apparatus 
according to the first aspect or the combination thereof 
with the temperature control apparatus according to one or 
more of the second and third aspects, the control input 
generated by the heater control means may include an input 
component depending on the temperature of the active ele- 
ment represented by the element temperature data, providing 
an invention based on the combination of the temperature 
control apparatus according to the first and fourth aspects 
or an invention based on the combination thereof with the 
temperature control apparatus according to one or more of 
the second and third aspects. Similarly, in the tempera- 
ture control method according to the first aspect or the 
combination thereof with the temperature control method ac- 
cording to one or more of the second and third aspects, 
when the control input is generated, the control input may 
be generated so as to further include an input component 
depending on the temperature of the active element repre- 
sented by the element temperature data, providing an inven- 
tion based on the combination of the temperature control 
methods according to the first and fourth aspects or an in- 
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vention based on the combination thereof with the tempera- 
ture control method according to one or more of the second 
and third aspects. Furthermore, in the recording medium 
according to the first aspect or the combination thereof 
with the recording medium according to one or more of the 
second and third aspects, the control input generating pro- 
gram may have an algorithm for enabling the computer to 
generate the control input so as to further include an in- 
put component depending on the temperature of the active 
element represented by the element temperature data, pro- 
viding an invention based on the combination of the re- 
cording mediums according to the first and fourth aspects 
or an invention based on the combination thereof with the 
recording medium according to one or more of the second and 
third aspects. 

By thus combining two or more of the first 
through fourth aspects, the advantages of those aspects are 
added together for increasing the stability of the tempera- 
ture of the active element with respect to the target tem- 
perature . 

In the temperature control apparatus according 
to the second aspect, the control input generated by the 
heater control means may include an input component depend- 
ing on the target temperature, providing an invention based 
on the combination of the temperature control apparatus ac- 
cording to the second and third aspects. Similarly, in the 
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temperature control method according to the second aspect, 
when the control input is generated, the control input may 
be generated so as to further include an input component 
depending on the target temperature, providing an invention 
based on the combination of the temperature control methods 
according to the second and third aspects. Furthermore, in 
the recording medium according to the second aspect, the 
control input generating program may have an algorithm for 
enabling the computer to generate the control input so as 
to further include an input component depending on the tar- 
get temperature, providing an invention based on the combi- 
nation of the recording mediums according to the second and 
third aspects. 

With the above arrangement, the advantages of 
the second and third aspects are added together for con- 
trolling the temperature of the active element more stably 
at the target temperature. 

In the temperature control apparatus according 
to the second aspect or the combination thereof with the 
temperature control apparatus according to the third as- 
pect, the control input generated by the heater control 
means may include an input component depending on the tem- 
perature of the active element represented by the element 
temperature data, providing an invention based on the com- 
bination of the temperature control apparatus according to 
the second and fourth aspects or an invention based on the 
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combination of the temperature control apparatus according 
to the second through fourth aspects. Similarly, in the 
temperature control method according to the second aspect 
or the combination thereof with the temperature control 
method according to the third aspect, when the control in- 
put is generated, the control input may be generated so as 
to further include an input component depending on the tem- 
perature of the active element represented by the element 
temperature data, providing an invention based on the com- 
bination of the temperature control methods according to 
the second and fourth aspects or an invention based on the 
combination of the temperature control methods according to 
the second through fourth aspects. Furthermore, in the re- 
cording medium according to the second aspect or the combi- 
nation thereof with the recording medium according to the 
third aspect, the control input generating program may have 
an algorithm for enabling the computer to generate the con- 
trol input so as to further include an input component de- 
pending on the temperature of the active element repre- 
sented by the element temperature data, providing an inven- 
tion based on the combination of the recording mediums ac- 
cording to the second and fourth aspects or an invention 
based on the combination of the recording mediums according 
to the second through fourth aspects. 

By thus combining those aspects, the advantages 
of the aspects are added together for increasing the sta- 



- 40 - 



bility of the temperature of the active element with re- 
spect to the target temperature. 

In the temperature control apparatus according 
to the third aspect, the control input generated by the 
heater control means may include an input component depend- 
ing on the temperature of the active element represented by 
the element temperature data, providing an invention based 
on the combination of the temperature control apparatus ac- 
cording to the third and fourth aspects. Similarly, in the 
temperature control method according to the third aspect, 
when the control input is generated, the control input may 
be generated so as to further include an input component 
depending on the temperature of the active element repre- 
sented by the element temperature data, providing an inven- 
tion based on the combination of the temperature control 
methods according to the third and fourth aspects. Fur- 
thermore, in the recording medium according to the third 
aspect, the control input generating program may have an 
algorithm for enabling the computer to generate the control 
input so as to further include an input component depending 
on the temperature of the active element represented by the 
element temperature data, providing an invention based on 
the combination of the recording mediums according to the 
third and fourth aspects. 

With the above combination, the advantages of 
the third and fourth aspects are added together for in- 
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creasing the stability of the temperature of the active 
element with respect to the target temperature. 

According to the present invention, it is pref- 
erable to combine all the first through third aspects or 
all the first through fourth aspects. If all the first 
through third aspects are combined, then the input compo- 
nent depending on the difference between the temperature of 
the active element and the target temperature should pref- 
erably comprise a component proportional to an accumulative 
sum (integrated value) of values of the difference at re- 
spective times and a component proportional to the differ- 
ence, for example. If all the first through fourth aspects 
are combined, then the input component depending on the 
difference between the temperature of the active element 
and the target temperature should preferably comprise a 
component proportional to an accumulative sum (integrated 
value) of values of the difference at respective times. 

The combination of the inventions according to 
the first through fourth aspects is also applicable to the 
inventions according to the fifth through eighth aspects, 
and two or more of the fifth through eighth aspects should 
preferably be combined. Specifically, if the temperature 
control apparatus according to the fifth and sixth aspects 
are combined, then the temperature control apparatus ac- 
cording to the fifth aspect is characterized by means for 
sequentially acquiring exhaust gas temperature data repre- 
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senting the temperature of the exhaust gas, wherein the 
control input generated by the heater control means in- 
cludes an input component depending on the temperature of 
the exhaust gas represented by the exhaust gas temperature 
data. Similarly, if the temperature control methods ac- 
cording to the fifth and sixth aspects are combined, then 
the temperature control method according to the fifth as- 
pect is characterized by the step of sequentially acquiring 
exhaust gas temperature data representing the temperature 
of the exhaust gas, wherein when the control input is gen- 
erated, the control input is generated so as to further in- 
clude an input component depending on the temperature of 
the exhaust gas represented by the exhaust gas temperature 
data. Furthermore, if the recording mediums according to 
the fifth and sixth aspects are combined, then the re- 
cording medium according to the fifth aspect is character- 
ized in that the temperature control program further in- 
cludes a program for enabling the computer to perform a 
process of sequentially acquiring exhaust gas temperature 
data representing the temperature of the exhaust gas, 
wherein the control input generating program has an algo- 
rithm for enabling the computer to generate the control in- 
put so as to further include an input component depending 
on the temperature of the exhaust gas represented by the 
exhaust gas temperature data. 
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By thus combining the fifth and sixth aspects, 
the advantages of those aspects are added together for con- 
trolling the temperature of the active element more stably 
at the target temperature and for controlling the tempera- 
ture of the active element more stably at a temperature 
corresponding to the target temperature for the heater. 

In the temperature control apparatus according 
to the fifth aspect or the combination thereof with the 
temperature control apparatus according to the sixth as- 
pect, the control input generated by the heater control 
means may include an input component depending on the tar- 
get temperature (a target value for the temperature of the 
heater) , providing an invention based on the combination of 
the temperature control apparatus according to the fifth 
and seventh aspects or an invention based on the combina- 
tion of the temperature control apparatus according to the 
fifth through seventh aspects. Similarly, in the tempera- 
ture control method, according to the fifth aspect or the 
combination thereof with the temperature control method ac- 
cording to the sixth aspect, when the control input is gen- 
erated, the control input may be generated so as to further 
include an input component depending on the target tempera- 
ture, providing an invention based on the combination of 
the temperature control methods according to the fifth and 
seventh aspects or an invention based on the combination of 
the temperature control methods according to the fifth 
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through seventh aspects. Furthermore, in the recording me- 
dium according to the fifth aspect or the combination 
thereof with the recording medium according to the sixth 
aspect , the control input generating program may have an 
algorithm for enabling the computer to generate the control 
input so as to further include an input component depending 
on the target temperature, providing an invention based on 
the combination of the recording mediums according to the 
fifth and sixth aspects or an invention based on the combi- 
nation of the recording mediums according to the fifth 
through seventh aspects. 

If the fifth and seventh aspects are thus com- 
bined, the advantages of those aspects are added together 
for controlling the temperature of the heater more stably 
at the target temperature and for controlling the tempera- 
ture of the active element more stably at a temperature 
corresponding to the target temperature for the heater. 
Particularly, if the fifth through seventh aspects are com- 
bined, the advantages of the fifth through seventh aspects 
are added together for appropriately increasing the stabil- 
ity of the temperature of the heater and the temperature of 
the active element. 

Moreover, in the temperature control apparatus 
according to the fifth aspect or the combination thereof 
with the temperature control apparatus according to one or 
more of the sixth and seventh aspects, the control input 
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generated by the heater control means may include an input 
component depending on the temperature of the heater repre- 
sented by the heater temperature data, providing an inven- 
tion based on the combination of the temperature control 
apparatus according to the fifth and eighth aspects or an 
invention based on the combination thereof with the tem- 
perature control apparatus according to one or more of the 
sixth and seventh aspects. Similarly, in the temperature 
control method according to the fifth aspect or the combi- 
nation thereof with the temperature control method accord- 
ing to one or more of the sixth and seventh aspects, when 
the control input is generated, the control input may be 
generated so as to further include an input component de- 
pending on the temperature of the heater represented by the 
heater temperature data, providing an invention based on 
the combination of the temperature control methods accord- 
ing to the fifth and eighth aspects or an invention based 
on the combination thereof with the temperature control 
method according to one or more of the sixth and seventh 
aspects. Furthermore, in the recording medium according to 
the fifth aspect or the combination thereof with the re- 
cording medium according to one or more of the sixth and 
seventh aspects, the control input generating program may 
have an algorithm for enabling the computer to generate the 
control input so as to further include an input component 
depending on the temperature of the heater represented by 
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the heater temperature data, providing an invention based 
on the combination of the recording mediums according to 
the fifth and eighth aspects or an invention based on the 
combination thereof with the recording medium according to 
one or more of the sixth and seventh aspects. 

By thus combining those aspects, the advantages 
of the aspects are added together for increasing the sta- 
bility of the temperature of the heater with respect to the 
target temperature and for controlling the temperature of 
the active element more stably at a temperature correspond- 
ing to the target temperature for the heater. 

In the temperature control apparatus according 
to the sixth aspect, the control input generated by the 
heater control means may include an input component depend- 
ing on the target temperature (a target value for the tem- 
perature of the heater) , providing an invention based on 
the combination of the temperature control apparatus ac- 
cording to the sixth and seventh aspects. Similarly, in 
the temperature control method according to the sixth as- 
pect, when the control input is generated, the control in- 
put may be generated so as to further include an input com- 
ponent depending on the target temperature, providing an 
invention based on the combination of the temperature con- 
trol methods according to the sixth and seventh aspects. 
Furthermore, in the recording medium according to the sixth 
aspect, the control input generating program may have an 
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algorithm for enabling the computer to generate the control 
input so as to further include an input component depending 
on the target temperature, providing an invention based on 
the combination of the recording mediums according to the 
sixth and seventh aspects. 

With the above arrangement, the advantages of 
the sixth and seventh aspects are added together for con- 
trolling the temperature of the heater more stably at the 
target temperature and for controlling the temperature of 
the active element more stably at a temperature correspond- 
ing to the target temperature for the heater. 

Moreover, in the temperature control apparatus 
according to the sixth aspect or the combination thereof 
with the temperature control apparatus according to the 
seventh aspect, the control input generated by the heater 
control means may include an input component depending on 
the temperature of the heater represented by the heater 
temperature data, providing an invention based on the com- 
bination of the temperature control apparatus according to 
the sixth and eighth aspects or an invention based on the 
combination thereof with the temperature control apparatus 
according to the sixth through eighth aspects. Similarly, 
in the temperature control method according to the sixth 
aspect or the combination thereof with the temperature con- 
trol method according to the seventh aspect, when the con- 
trol input is generated, the control input may be generated 
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so as to further include an input component depending on 
the temperature of the heater represented by the heater 
temperature data, providing an invention based on the com- 
bination of the temperature control methods according to 
the sixth and eighth aspects or an invention based on the 
combination thereof with the temperature control method ac- 
cording to the sixth through eighth aspects. Furthermore, 
in the recording medium according to the sixth aspect or 
the combination thereof with the recording medium according 
to the seventh aspect, the control input generating program 
may have an algorithm for enabling the computer to generate 
the control input so as to further include an input compo- 
nent depending on the temperature of the heater represented 
by the heater temperature data, providing an invention 
based on the combination of the recording mediums according 
to the sixth and eighth aspects or an invention based on 
the combination thereof with the recording medium according 
to the sixth through eighth aspects. 

By thus combining those aspects, the advantages 
of the aspects are added together for increasing the sta- 
bility of the temperature of the heater with respect to the 
target temperature and for controlling the temperature of 
the active element more stably at a temperature correspond- 
ing to the target temperature for the heater. 

Moreover, in the temperature control apparatus 
according to the seventh aspect, the control input gener- 
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ated by the heater control means may include an input com- 
ponent depending on the temperature of the heater repre- 
sented by the heater temperature data, providing an inven- 
tion based on the combination of the temperature control 
apparatus according to the seventh and eighth aspects. 
Similarly, in the temperature control method according to 
the seventh aspect, when the control input is generated, 
the control input may be generated so as to further include 
an input component depending on the temperature of the 
heater represented by the heater temperature data, provid- 
ing an invention based on the combination of the tempera- 
ture control methods according to the seventh and eighth 
aspects- Furthermore, in the recording medium according to 
the seventh aspect, the control input generating program 
may have an algorithm for enabling the computer to generate 
the control input so as to further include an input compo- 
nent depending on the temperature of the heater represented 
by the heater temperature data, providing an invention 
based on the combination of the recording mediums according 
to the seventh and eighth aspects. 

With this combination, the advantages of the 
seventh and eighth aspects are added together for increas- 
ing the stability of the temperature of the heater with re- 
spect to the target temperature and for controlling the 
temperature of the active element more stably at a tempera- 
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ture corresponding to the target temperature for the 
heater. 

It is preferable to combine all the fifth 
through seventh aspects or all the fifth through eighth as- 
pects in particular. If all the fifth through seventh as- 
pects are combined, then the input component depending on 
the difference between the temperature of the heater and 
the target temperature should preferably comprise a compo- 
nent proportional to an accumulative sum (integrated value) 
of values of the difference at respective times and a com- 
ponent proportional to the difference, for example. If all 
the fifth through eighth aspects are combined, then the in- 
put component depending on the difference between the tem- 
perature of the heater and the target temperature should 
preferably comprise a component proportional to an accumu- 
lative sum (integrated value) of values of the difference 
at respective times. 

In the temperature control apparatus according 
to the second aspect wherein the control input includes an 
input component depending on the temperature of the exhaust 
gas or in the invention based on the combination of the 
temperature control apparatus according to one or more of 
the first, third, and fourth aspects and the temperature 
control apparatus according to the second aspect, the input 
component depending on the temperature of the exhaust gas, 
in the control input sequentially generated by the heater 
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control means, comprises an input component depending on 
time-series data of the temperature of the exhaust gas in- 
cluding a present value of the temperature of the exhaust 
gas and a future value of the temperature of the exhaust 
gas after a first predetermined time, wherein the heater 
control means preferably generates the control input in- 
cluding the input component according to a predictive con- 
trol algorithm. Similarly, in the temperature control 
method according to the second aspect wherein the control 
input includes an input component depending on the tempera- 
ture of the exhaust gas or in the invention based on the 
combination of the temperature control method according to 
one or more of the first, third, and fourth aspects and the 
temperature control method according to the second aspect, 
the input component depending on the temperature of the ex- 
haust gas, included in the control input, comprises an in- 
put component depending on time-series data of the tempera- 
ture of the exhaust gas including a present value of the 
temperature of the exhaust gas and a future value of the 
temperature of the exhaust gas after a first predetermined 
time, wherein the control input including the input compo- 
nent is preferably generated according to a predictive con- 
trol algorithm. Furthermore, in the recording medium ac- 
cording to the second aspect wherein the control input in- 
cludes an input component depending on the temperature of 
the exhaust gas or in the invention based on the combina- 
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tion of the recording medium according to one or more of 
the first, third, and fourth aspects and the recording me- 
dium according to the second aspect, the input component 
depending on the temperature of the exhaust gas, included 
in the control input, comprises an input component depend- 
ing on time-series data of the temperature of the exhaust 
gas including a present value of the temperature of the ex- 
haust gas and a future value of the temperature of the ex- 
haust gas after a first predetermined time, wherein the 
control input generating program for enabling the computer 
to generate the control input including the input component 
preferably has a predictive control algorithm. 

The foregoing also applies to the invention ac- 
cording to the sixth aspect or the invention based on the 
combination of one or more of the fifth, seventh, and 
eighth aspects and the sixth aspect, with respect to any of 
the temperature control apparatus, the temperature control 
method, and the recording medium according to the present 
invention for generating a control input to equalize the 
temperature of the heater to the target temperature. 

With the above arrangement, the input component 
included in the control input which depends on the tempera- 
ture of the exhaust gas is an input component depending on 
not only the present value of the temperature of the ex- 
haust gas, but also the time-series data of the temperature 
of the exhaust gas including at least the future value af- 
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ter the first predetermined time. Therefore, the invention 
according to the second aspect or an invention including 
the same is capable of minimizing a change of the tempera- 
ture of the active element with respect to a change of the 
temperature of the exhaust gas. Therefore, the stability 
with which to control the temperature of the active element 
at the target temperature can effectively be increased ir- 
respective of variations of the temperature of the exhaust 
gas. Similarly, the invention according to the sixth as- 
pect or an invention including the same is capable of mini- 
mizing a change of the temperature of the heater with re- 
spect to a change of the temperature of the exhaust gas . 
Therefore, the stability with which to control the tempera- 
ture of the heater at the target temperature can effec- 
tively be increased irrespective of variations of the tem- 
perature of the exhaust gas, and the temperature of the ac- 
tive element can more stably be controlled at a temperature 
corresponding to the target temperature for the heater. 
The time -series data may include not only the present value 
of the temperature of the exhaust gas and the future value 
thereof after the first predetermined time, but also a plu- 
rality of future values of the temperature of the exhaust 
gas from present to a time after the first predetermined 
time . 

If the control input is generated according to 
the predictive control algorithm taking into account the 
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future value of the temperature of the exhaust gas, then 
the future value of the temperature of the exhaust gas 
which is required to generate the control input may be es- 
timated according to a suitable algorithm. However, in 
general, the temperature of the exhaust gas does not change 
appreciably abruptly. In the temperature control apparatus 
according to the second aspect or the sixth aspect, the 
heater control means may generate the control input by de- 
termining the future value of the temperature of the ex- 
haust gas until after the first predetermined time as being 
equal to the present value of the temperature of the ex- 
haust gas. Similarly, in the temperature control method 
according to the second aspect or the sixth aspect, the 
predictive control algorithm may comprise an algorithm for 
generating the control input by determining the future 
value of the temperature of the exhaust gas until after the 
first predetermined time as being equal to the present 
value of the temperature of the exhaust gas. Furthermore, 
in the recording medium according to the second aspect or 
the sixth aspect, the algorithm of the control input gener- 
ating program may enable the computer to generate the con- 
trol input by determining the future value of the tempera- 
ture of the exhaust gas until after the first predetermined 
time as being equal to the present value of the temperature 
of the exhaust gas. 
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With the above arrangement, since the present 
value of the temperature of the exhaust gas may sequen- 
tially be grasped, the algorithm for controlling the tem- 
perature of the exhaust gas sensor, including a process of 
grasping the temperature of the exhaust gas (acquiring the 
data representing the temperature of the exhaust gas), can 
easily be constructed. 

In the temperature control apparatus according 
to the third aspect wherein the control input includes an 
input component depending on the target temperature for the 
active element or in the invention based on the combination 
of the temperature control apparatus according to one or 
more of the first, second, and fourth aspects and the tem- 
perature control apparatus according to the third aspect, 
the input component depending on the target temperature 
which is sequentially generated by the heater control means 
comprises an input component depending on time -series data 
of the target temperature including a present value of the 
target temperature and a future value of the target tem- 
perature after a second predetermined time, wherein the 
heater control means preferably generates the control input 
including the input component according to a predictive 
control algorithm. Similarly, in the temperature control 
method according to the third aspect or in the invention 
based on the combination of the temperature control method 
according to one or more of the first, second, and fourth 
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aspects and the temperature control method according to the 
third aspect, the input component depending on the target 
temperature, included in the control input, comprises an 
input component depending on time- series data of the target 
temperature including a present value of the target tem- 
perature and a future value of the target temperature after 
a second predetermined time, wherein the control input in- 
cluding the input component is preferably generated accord- 
ing to a predictive control algorithm. Furthermore, in the 
recording medium according to the third aspect or in the 
invention based on the combination of the recording medium 
according to one or more of the first, second, and fourth 
aspects and the recording medium according to the third as- 
pect , the input component depending on the target tempera- 
ture, included in the control input, comprises an input 
component depending on time-series data of the target tem- 
perature including a present value of the target tempera- 
ture and a future value of the target temperature after a 
second predetermined time, wherein the control input gener- 
ating program for enabling the computer to generate the 
control input including the input component preferably has 
a predictive control algorithm. 

The foregoing also applies to the invention ac- 
cording to the seventh aspect or the invention based on the 
combination of one or more of the fifth, sixth, and eighth 
aspects and the seventh aspect, with respect to any of the 
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temperature control apparatus, the temperature control 
method, and the recording medium according to the present 
invention for generating a control input to equalize the 
temperature of the heater to the target temperature. 

With the above arrangement, the input component 
included in the control input which depends on the target 
temperature is an input component depending on not only the 
present value of the target temperature, but also the time- 
series data of the target temperature including at least 
the future value after the second predetermined time. 
Therefore, the invention according to the third aspect or 
an invention including the same is capable of preventing 
the temperature of the active element from overshooting the 
target temperature therefor when the target temperature for 
the active element is changed. The rate at which the tem- 
perature of the active element converges to the target tem- 
perature can be increased. As a result, the temperature of 
the active element can quickly and smoothly follow the tar- 
get temperature. In particular, immediately after the in- 
ternal combustion engine has started to operate, the tem- 
perature of the active element can quickly be converted to 
and stabilized at the target temperature. Therefore, imme- 
diately after the internal combustion engine has started to 
operate, the output characteristics of the active element 
can quickly be stabilized at desired characteristics. As 
the ability of the temperature of the active element to 
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follow a change of the target temperature for the active 
element is increased, the target temperature for the active 
element can variably be set to a value that is suitable for 
the operating state of the internal combustion engine. 

Similarly, the invention according to the sev- 
enth aspect or an invention including the same is capable 
of preventing the temperature of the heater from overshoot- 
ing the target temperature therefor when the target tem- 
perature for the heater is changed. The rate at which the 
temperature of the heater converges to the target tempera- 
ture can be increased. As a result, the temperature of the 
heater can quickly and smoothly follow the target tempera- 
ture. In particular, immediately after the internal com- 
bustion engine has started to operate, the temperature of 
the heater can quickly be converted to and stabilized at 
the target temperature, and the temperature of the active 
element can quickly be converted and stabilized at a tem- 
perature corresponding to the target temperature for the 
heater. Therefore, immediately after the internal combus- 
tion engine has started to operate, the output characteris- 
tics of the active element can quickly be stabilized at de- 
sired characteristics. As the ability of the temperature 
of the heater to follow a change of the target temperature 
for the heater is increased, the target temperature for the 
heater and a desired temperature of the active element can 
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variably be set to a value that is suitable for the operat- 
ing state of the internal combustion engine. 

The time-series data of the target temperature 
may include not only the present value of the target tem- 
perature and the future value thereof after the second pre- 
determined time, but also a plurality of future values of 
the target temperature from present to a time after the 
second predetermined time. 

In the temperature control apparatus according 
to the first through fourth aspects of the present inven- 
tion, the heater control means preferably generates the 
control input according to an optimum control algorithm. 
Similarly, in temperature control methods according to the 
first through fourth aspects, the control input is prefera- 
bly generated according to an optimum control algorithm. 
In the recording medium according to the first through 
fourth aspects, the algorithm of the control input generat- 
ing program preferably comprises an optimum control algo- 
rithm. 

With the above arrangement, according to the 
first through fourth aspects, it is possible to generate 
the control input for minimizing variations of the tempera- 
ture of the active element and the control input while 
holding in balance the ability of the temperature of the 
active element to converge to the target temperature and 
variations of the control input (variations of the heat 
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generating energy supplied to the heater) . As a result, 
the temperature of the active element can more stably be 
controlled at the target temperature (desired temperature) , 
and at the same time, the amount of the heat generating en- 
ergy supplied to the heater can be held to a required 
limit . 

When the control input is generated according to 
the optimum control algorithm in the first through fourth 
aspects of the present invention, it is preferable to de- 
termine in advance a model to be controlled which has, as 
state quantities to be controlled, the difference between 
the temperature of the active element and the target tem- 
perature, a change per given time of the difference (a rate 
of change of the difference) , and a change per given time 
of the temperature of the heater (a rate of change of the 
temperature of the heater) , for example, and which also 
has, as an input quantity, at least a change per given time 
of the control input (a rate of change of the control in- 
put) , and to construct the optimum control algorithm based 
on the model to be controlled. Alternatively, it is also 
preferable to determine in advance a model to be controlled 
which has, as state quantities to be controlled, the dif- 
ference between the temperature of the active element and 
the target temperature, a change per given time of the tem- 
perature of the active element (a rate of change of the 
temperature of the active element) , and a change per given 
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time of the temperature of the heater (a rate of change of 
the temperature of the heater) , for example, and which also 
has, as an input quantity, at least a change per given time 
of the control input (a rate of change of the control in- 
put) , and to construct the optimum control algorithm based 
on the model to be controlled. At any rate, there is de- 
termined a control input capable of minimizing an evaluat- 
ing function (which is represented as an integrated value 
(accumulated sum) of weighted sum of the squares of the 
above state quantities and the square of the change per 
given time of the control input) of the model to be con- 
trolled- With this arrangement, since the control input 
can be generated in a manner to not only eliminate the dif- 
ference between the temperature of the active element and 
the target temperature, but also suppress variations of the 
state quantities of the model to be controlled and the con- 
trol input while keeping them in balance, the stability of 
the temperature of the active element with respect to the 
target temperature can be increased. 

If the model to be controlled, which serves as a 
basis for the optimum control algorithm in the first 
through fourth aspects has, as state quantities to be con- 
trolled, the difference between the temperature of the ac- 
tive element and the target temperature, the change per 
given time of the difference, and the change per given time 
of the temperature of the heater, it is possible to gener- 
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ate the control input which includes the input component 
depending on the difference between the temperature of the 
active element and the target temperature and the input 
component depending on the temperature of the heater. If 
the model to be controlled has, as state quantities to be 
controlled, the difference between the temperature of the 
active element and the target temperature, the change per 
given time of the temperature of the active element, and 
the change per given time of the temperature of the heater, 
it is possible to generate the control input which includes 
the input component depending on the difference between the 
temperature of the active element and the target tempera- 
ture, the input component depending on the temperature of 
the heater, and the input component depending on the tem- 
perature of the active element. By including, as input 
quantities for the model to be controlled, a change per 
given time of the temperature of the exhaust gas and a 
change per given time of the target temperature for the 
temperature of the active element, in addition to the 
change per given time of the control input, it is possible 
to generate the control input including the input component 
depending on the temperature of the exhaust gas and the in- 
put component depending on the target temperature for the 
active element. 

In the temperature control apparatus according 
to the fifth through eighth aspects of the present inven- 
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tion, the heater control means preferably generates the 
control input according to an optimum control algorithm. 
Similarly, in temperature control methods according to the 
fifth through eighth aspects, the control input is prefera- 
bly generated according to an optimum control algorithm. 
In the recording medium according to the fifth through 
eighth aspects, the algorithm of the control input generat- 
ing program preferably comprises an optimum control algo- 
rithm. 

With the above arrangement, it is possible to 
generate the control input for minimizing variations of the 
temperature of. the heater and the control input while hold- 
ing in balance the ability of the temperature of the heater 
to converge to the target temperature and variations of the 
control input. As a result, the temperature of the heater 
can more stably be controlled at the target temperature and 
the temperature of the active element can more stably be 
controlled at a temperature (desired temperature) corre- 
sponding to the target temperature for the heater. At the 
same time, the amount of the heat generating energy sup- 
plied to the heater can be held to a required limit. 

When the control input is generated according to 
the optimum control algorithm in the fifth through eighth 
aspects of the present invention, it is preferable to de- 
termine in advance a model to be controlled which has, as 
state quantities to be controlled, the difference between 
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the temperature of the heater and the target temperature, a 
change per given time of the difference (a rate of change 
of the difference) , and a change per given time of the tem- 
perature of the active element (a rate of change of the 
temperature of the active element), for example, and which 
also has, as an input quantity, at least a change per given 
time of the control input (a rate of change of the control 
input) , and to construct the optimum control algorithm 
based on the model to be controlled. Alternatively, it is 
also preferable to determine in advance a model to be con- 
trolled which has, as state quantities to be controlled, 
the difference between the temperature of the heater and 
the target temperature, a change per given time of the tem- 
perature of the active element (a rate of change of the 
temperature of the active element), and a change per given 
time of the temperature of the heater (a rate of change of 
the temperature of the heater), for example, and which also 
has, as an input quantity, at least a change per given time 
of the control input (a rate of change of the control in- 
put) , and to construct the optimum control algorithm based 
on the model to be controlled. At any rate, there is de- 
termined a control input capable of minimizing an evaluat- 
ing function of the model to be controlled. With this ar- 
rangement, since the control input can be generated in a 
manner to not only eliminate the difference between the 
temperature of the heater and the target temperature, but 
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also suppress variations of the difference, the temperature 
of the active element, and the control input while keeping 
them in balance, the stability of the temperature of the 
heater with respect to the target temperature and the sta- 
bility of the temperature of the active element can be in- 
creased. 

If the model to be controlled, which serves as a 
basis for the optimum control algorithm in the fifth 
through eighth aspects has, as state quantities to be con- 
trolled, the difference between the temperature of the 
heater and the target temperature, the change per given 
time of the difference, and the change per given time of 
the temperature of the active element, it is possible to 
generate the control input which includes the input compo- 
nent depending on the difference between the temperature of 
the heater and the target temperature and the input compo- 
nent depending on the temperature of the active element. 
If the model to be controlled has, as state quantities to 
be controlled, the difference between the temperature of 
the heater and the target temperature, the change per given 
time of the temperature of the active element, and the 
change per given time of the temperature of the heater, it 
is possible to generate the control input which includes 
the input component depending on the difference between the 
temperature of the heater and the target temperature, the 
input component depending on the temperature of the heater, 
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and the input component depending on the temperature of the 
active element. By including, as input quantities for the 
model to be controlled, a change per given time of the tem- 
perature of the exhaust gas and a change per given time of 
the target temperature for the temperature of the heater, 
in addition to the change per given time of the control in- 
put, it is possible to generate the control input including 
the input component depending on the temperature of the ex- 
haust gas and the input component depending on the target 
temperature for the heater. 

When the optimum control algorithm according to 
the first through fourth aspects of the present invention 
and the predictive control algorithm for the temperature of 
the exhaust gas or the target temperature are combined, the 
control input is generated according to an optimum predic- 
tive control algorithm. When the control input is gener- 
ated according to an optimum predictive control algorithm, 
the temperature of the active element can be controlled at 
the target temperature with high stability while keeping in 
balance and suppressing variations of the state quantities 
of the model to be controlled, which serves as a basis for 
the optimum control algorithm, and the input quantity. The 
optimum predictive control algorithm can be constructed 
based on a model to be controlled which includes a change 
per given time of the temperature of the exhaust gas and/or 
a change per given time of the target temperature for the 
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active element, in addition to the change per given time of 
the control input, as the input quantities of the model to 
be controlled which has been given by way of example with 
respect to the first through fourth aspects, for example. 

The foregoing also applies to a combination of 
the optimum control algorithm according to the fifth 
through eighth aspects of the present invention and the 
predictive control algorithm for the temperature of the ex- 
haust gas or the target temperature. In this case, the 
control input is generated according to an optimum predic- 
tive control algorithm. When the control input is gener- 
ated according to an optimum predictive control algorithm, 
the temperature of the heater can be controlled at the tar- 
get temperature with high stability while keeping in bal- 
ance and suppressing variations of the state quantities of 
the model to be controlled, which serves as a basis for the 
optimum control algorithm, and the input quantity. The op- 
timum predictive control algorithm in this case can be con- 
structed based on a model to be controlled which includes a 
change per given time of the temperature of the exhaust gas 
and/or a change per given time of the target temperature, 
in addition to the change per given time of the control in- 
put, as the input quantities of the model to be controlled 
which has been given by way of example with respect to the 
fifth through eighth aspects, for example. 
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In the temperature control apparatus, the tem- 
perature control method, and the recording medium according 
to the present invention as described above in any of the 
aspects, the target temperature (a target temperature for 
the active element or the heater) in a period immediately 
after the internal combustion engine has started to operate 
until a third predetermined time elapses is preferably set 
to a temperature which is lower than the target temperature 
after elapse of the period immediately after the internal 
combustion engine has started to operate. 

Specifically, if the target temperature for the 
active element or the heater is high immediately after the 
internal combustion engine has started to operate, then 
when water in the exhaust passage is applied to the active 
element of the exhaust gas sensor, the active element tends 
to be damaged due to thermal stresses or the like which are 
caused by rapid heating of the active element. According 
to the present invention, the target temperature for the 
active element or the heater immediately after the internal 
combustion engine has started to operate is set to a tem- 
perature which is lower than the target temperature after 
elapse of the period immediately after the internal combus- 
tion engine has started to operate. In this manner, the 
active element of the exhaust gas sensor is prevented from 
being damaged by thermal stresses or the like. 
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In the temperature control apparatus according 
to the present invention as described above in any of the 
aspects, the heater comprises an electric heater for gener- 
ating heat when energized by a battery according to a pulse 
width control (PWM control) process, if the control input 
generated by the heater control means comprises a duty cy- 
cle (the ratio of a pulse duration to one period of a 
pulsed signal used in the pulse with control process) in 
the pulse width control process, the temperature control 
apparatus is preferably further characterized by means for 
correcting the duty cycle depending on the voltage of the 
battery. Similarly, in the temperature control method ac- 
cording to the present invention in any of the aspects, the 
heater comprises an electric heater for generating heat 
when energized by a battery according to a pulse width con- 
trol (PWM control) process, and the control input which is 
generated comprises a duty cycle in the pulse width control 
process, the temperature control method is preferably fur- 
ther characterized by the step of correcting the duty cycle 
depending on the voltage of the battery. Furthermore, in 
the recording medium according to the present invention in 
any of the aspects, the heater comprises an electric heater 
for generating heat when energized by a battery according 
to a pulse width control (PWM control) process, and the 
control input which is generated by the computer according 
to the control input program comprises a duty cycle in the 
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pulse width control process , wherein the temperature con- 
trol program preferably further includes a program for ena- 
bling the computer to perform a process of correcting the 
duty cycle generated by the control input generating pro- 
gram depending on the voltage of the battery. 

Specifically, if the voltage of the battery is 
substantially constant, then the duty cycle defines an 
amount of electric power as an amount of heat generating 
energy supplied to the heater. If the voltage of the bat- 
tery is varied by an alternator or the like, then the 
amount of electric power supplied to the heater is affected 
by not only the duty cycle but also the voltage of the bat- 
tery. Therefore, by correcting the duty cycle as the con- 
trol input depending on the voltage of the battery, varia- 
tions of the voltage of the battery in controlling the tem- 
perature of the active element or the heater at the target 
temperature can be compensated for. As a result, the tem- 
perature of the active element or the temperature of the 
heater can stably be controlled at the target temperature 
without being affected by variations of the voltage of the 
battery, and the stability of the temperature of the active 
element can be increased. 

According to the present invention, the exhaust 
gas sensor may comprise an 0 2 sensor disposed downstream of 
a catalytic converter for purifying the exhaust gas, for 
example. If the air-fuel ratio of the exhaust gas is con- 
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trolled to keep the output voltage of the 0 2 sensor at a 
predetermined level in order for the catalytic converter to 
perform its desired exhaust gas purifying capability, the 
temperature of the active element of the 0 2 sensor should 
preferably basically be controlled at a temperature equal 
to or higher than 750°C (e.g., 800°C) for better air-fuel 
ratio control. In this case, when the heater is to be con- 
trolled with a target temperature determined for the active 
element, the target temperature may be set to a temperature 
equal to or higher than 750°C (e.g., 800°C). When the 
heater is to be controlled with a target temperature deter- 
mined for the heater, the target temperature may be set to 
a temperature equal to or higher than 850°C (e.g., 900°C). 

Furthermore, immediately after the internal com- 
bustion engine has started to operate (in a period in which 
a third predetermined time (e.g., 15 seconds) elapses after 
the internal combustion engine has started to operate) , the 
temperature of the active element should preferably be con- 
trolled at a temperature (e.g., 600°C) lower than the above 
temperature in order to prevent the active element of the 
0 2 sensor from being damaged. In this case, when the 
heater is to be controlled with a target temperature deter- 
mined for the active element, the target temperature for 
the active element immediately after the internal combus- 
tion engine has started to operate may be set to 600 °C, for 
example. When the heater is to be controlled with a target 
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temperature determined for the heater, the target tempera- 
ture for the heater immediately after the internal combus- 
tion engine has started to operate may be set to 700°C, for 
example . 

Brief description of the drawings: 

FIG* 1 is a block diagram of an apparatus ac- 
cording to a first embodiment of the present invention; 
FIG. 2 is a fragmentary cross-sectional view showing a 
structure of an 0 2 sensor (exhaust gas sensor) in the appa- 
ratus shown in FIG. 1; FIG. 3 is a graph illustrative of 
the output characteristics of the 0 2 sensor shown in FIG. 
2; FIG. 4 is a block diagram showing a functional arrange- 
ment of a sensor temperature control means in the apparatus 
shown in FIG. 1; FIG. 5 is a cross-sectional view showing a 
processing operation of an exhaust temperature observer in 
the sensor temperature control means shown in FIG. 4; FIG. 
6 is a block diagram showing a functional arrangement of 
the exhaust temperature observer in the sensor temperature 
control means shown in FIG. 4; FIG. 7 is a block diagram 
showing a functional arrangement of a heater controller in 
the sensor temperature control means shown in FIG. 5; FIG. 
8 is a flowchart of an overall processing sequence of the 
sensor temperature control means in the apparatus shown in 
FIG. 1; FIGS. 9 through 11 are flowcharts of subroutines of 
the flowchart shown in FIG. 8; FIGS. 12 and 13 are graphs 
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showing the results of simulations of the first embodiment; 
and FIG. 14 is a block diagram showing a functional ar- 
rangement of a sensor temperature control means in an appa- 
ratus according to a second embodiment of the present in- 
vention. 

Best mode for carrying out the invention: 

A first embodiment of the present invention will 
be described below with reference to FIGS. 1 through 13. 
FIG. 1 shows in block form an overall arrangement of the 
apparatus according to the first embodiment of the present 
invention. In FIG. 1, an engine (an internal combustion 
engine) 1 mounted on an automobile, a hybrid vehicle, or 
the like combusts a mixture of fuel and air and generates 
an exhaust gas, which is discharged into the atmosphere 
through an exhaust passage 3 communicating with an exhaust 
port 2 of the engine 1. The exhaust passage 3 incorporates 
therein two catalytic converters 4, 5 disposed successively 
downstream for purifying the exhaust gas emitted from the 
engine 1 and flowing through the exhaust passage 3. The 
exhaust passage 3 includes a section upstream of the cata- 
lytic converter 4 (between the exhaust port 2 and the cata- 
lytic converter 4), a section between the catalytic con- 
verters 4, 5, and a section downstream of the catalytic 
converter 5. These sections of the exhaust passage 3 are 
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provided by respective exhaust pipes 6a, 6b, 6c each in the 
form of a tubular passage-defining member. 

Each of the catalytic converters 4, 5 contains a 
catalyst 7 (three-way catalyst in the present embodiment) . 
The catalyst 7 has a passage-defining honeycomb structure 
and allows the exhaust gas to flow therethrough. Though 
the catalytic converters 4, 5 may be of a unitary structure 
with two catalytic beds, each comprising a three-way cata- 
lyst, disposed respectively in upstream and downstream re- 
gions thereof. 

In the present embodiment, the air-fuel ratio in 
the exhaust gas emitted from the engine 1 is controlled in 
order for the upstream catalytic converter 4, in particu- 
lar, to have a good exhaust gas purifying capability (the 
ability of the catalytic converter 4 to purify CO, HC, and 
NOx) . For controlling the air-fuel ratio in the exhaust 
gas, an 0 2 sensor 8 is mounted on the exhaust passage 3 be- 
tween the catalytic converters 4, 5, i.e., on the exhaust 
passage defined by the exhaust pipe 6b, and a wide-range 
air-fuel ratio sensor 9 is mounted on the exhaust passage 3 
upstream of the catalytic converter 4, i.e., on the exhaust 
passage defined by the exhaust pipe 6a. When the catalytic 
converters 4, 5 take a unitary structure with two catalytic 
beds, the 0 2 sensor 8 may be mounted between the upstream 
catalytic bed and the downstream catalytic bed. 
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The 0 2 sensor 8 corresponds to an exhaust gas 
sensor according to the present invention. Basic struc- 
tural details and characteristics of the 0 2 sensor 8 will 
be described below. As shown in FIG. 2, the 0 2 sensor 8 
has an active element 10 (sensitive element) in the form of 
a hollow bottomed cylinder made primarily of a solid elec- 
trolyte permeable to oxygen ions, e.g., stabilized zirconia 
(Zr0 2 + Y 2 0 3 ) - The active element 10 has outer and inner 
surfaces coated with porous platinum electrodes 11, 12, re- 
spectively. The 0 2 sensor 8 also has a rod-shaped ceramic 
heater 13 inserted as an electric heater into the active 
element 10 for heating the active element 10 for activation 
and controlling the temperature of the active element 10. 
The active element 10 is filled with air containing oxygen 
at a constant concentration, i.e., under a constant partial 
pressure, in a space around the ceramic heater 13. The 0 2 
sensor 8 is placed in a sensor casing 14 mounted on the ex- 
haust pipe 6b such that the tip end of the active element 
10 has its outer surface positioned in contact with the ex- 
haust gas flowing in the exhaust pipe 6b. 

As shown in FIG. 2, the tip end of the active 
element 10 is covered with a tubular protector 15 which 
protects the active element 10 against the impingement of 
foreign matter thereon. The tip end of the active element 
10 which is positioned in the exhaust pipe 6b contacts the 
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exhaust gas through a plurality of holes (not shown) de- 
fined in the protector 15. 

The 0 2 sensor 8 thus constructed operates as 
follows: An electromotive force depending on the concen- 
tration of oxygen in the exhaust gas is generated between 
the platinum electrodes 11, 12 based on the difference be- 
tween the concentration of oxygen in the exhaust gas which 
is brought into contact with the outer surface of the tip 
end of the active element 10 and the concentration of oxy- 
gen in the air in the active element 10. The generated 
electromotive force is amplified by an amplifier (not 
shown) , and then produced as the output voltage Vout from 
the 0 2 sensor 8 . 

The output voltage Vout of the 0 2 sensor 8 has 
characteristics (output characteristics) with respect to 
the concentration of oxygen in the exhaust gas or the air- 
fuel ratio in the exhaust gas which is recognized from the 
concentration of oxygen, as represented by a solid-line 
curve "a" (so-called "Z curve") in FIG. 3. The solid-line 
curve "a" in FIG. 3 represents the output characteristics 
of the 0 2 sensor 8 when the temperature of the active ele- 
ment 10 is 800 °C. The relationship between the temperature 
of the active element 10 and the output characteristics of 
the 0 2 sensor 8 will be described later on. 

As indicated by the curve "a" in FIG. 3, the 
output characteristics of the 0 2 sensor 8 are generally of 
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such a nature that the output voltage Vout changes substan- 
tially linearly with high sensitivity with respect to the 
air-fuel ratio of the exhaust gas only when the air-fuel 
ratio represented by the concentration of oxygen in the ex- 
haust gas is present in a narrow air-fuel ratio range A 
near a stoichiometric air-fuel ratio. In the air-fuel ra- 
tio range A (hereinafter referred to as "high-sensitivity 
air-fuel ratio range A"), the gradient of a change in the 
output voltage Vout with respect to a change in the air- 
fuel ratio, i.e., the gradient of the curve of the output 
characteristics of the 0 2 sensor 8, is large. In an air- 
fuel ratio range richer than the high-sensitivity air-fuel 
ratio range A and an air-fuel ratio range leaner than the 
high-sensitivity air-fuel ratio range A, the gradient of a 
change in the output voltage Vout with respect to a change 
in the air-fuel ratio, i.e., the gradient of the curve of 
the output characteristics of the 0 2 sensor 8, is smaller. 

The wide-range air-fuel ratio sensor 9, which 
will not be described in detail below, comprises an air- 
fuel ratio sensor disclosed in Japanese laid-open patent 
publication No. 4-369471 by the applicant of the present 
application, for example. The wide-range air-fuel ratio 
sensor 9 is a sensor for generating an output voltage KACT 
which changes linearly with respect to the air-fuel ratio 
in the exhaust gas in the air-fuel ratio wider in range 
than the 0 2 sensor 8. The output voltage Vout of the 0 2 
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sensor 8 and the output voltage KACT of the wide-range air- 
fuel ratio sensor 9 will hereinafter be referred to as 
"output Vout" and "output KACT", respectively. 

The apparatus according to the present embodi- 
ment also has a control unit 16 for controlling the air- 
fuel ratio in the exhaust gas and controlling the tempera- 
ture of the active element 10 of the 0 2 sensor 8. The con- 
trol unit 16 comprises a microcomputer including a CPU, a 
RAM, and a ROM (not shown) . For carrying out a control 
process to be described later on, the control unit 16 is 
supplied with the outputs Vout and KACT from the 0 2 sensor 
8 and the wide-range air-fuel ratio sensor 9, and also with 
data of detected values representing the rotational speed 
NE of the engine 1, the intake pressure PB (specifically, 
the absolute pressure in the intake pipe of the engine 1), 
the atmospheric temperature TA, the engine temperature TW 
(specifically, the temperature of the coolant of the engine 
1), etc., from sensors (not shown) combined with the engine 
1. The control unit 16 is also supplied with data of de- 
tected values representing the voltage VB (hereinafter re- 
ferred to as "battery voltage VB") of a battery (not shown) 
as a power supply of electronic accessories including an 
ignition unit (not shown) of the engine 1, the control unit 
16, the ceramic heater 13, etc., from a sensor (not shown). 
The ROM of the control unit 16 corresponds to a recording 
medium according to the present invention. 
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The control unit 16 has as its functional means 
an air-fuel ratio control means 17 for controlling the air- 
fuel ratio in the exhaust gas emitted from the engine 1, 
and a sensor temperature control means 18 for controlling 
the temperature of the active element 10 of the 0 2 sensor 
8. 

The air-fuel ratio control means 17 controls the 
air-fuel ratio in the exhaust gas supplied from the engine 
1 to the catalytic converter 4 in order to achieve a good 
purifying ability (purification rate) of the catalytic con- 
verter 4 to purify CO (carbon monoxide), HC (hydrocarbon), 
and NOx (nitrogen oxide) . When the 0 2 sensor 8 of the 
above output characteristics is disposed downstream of the 
catalytic converter 4, a good purifying ability of the 
catalytic converter 4 to purify CO, HC, and NOx can be 
achieved irrespective of the deteriorated state of the 
catalytic converter 4 by controlling the air-fuel ratio in 
the exhaust gas supplied to the catalytic converter 4, 
i.e., the air-fuel ratio in the exhaust gas upstream of the 
catalytic converter 4, to settle the output Vout of the 0 2 
sensor 8 at a certain predetermined value Vop (see FIG. 3) . 

Specifically, the air-fuel ratio control means 
17 uses the predetermined value Vop as a target value for 
the output Vout of the 0 2 sensor 8, and controls the air- 
fuel ratio in the exhaust gas supplied from the engine 1 to 
the catalytic converter 4 in order to settle and keep the 
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output Vout of the 0 2 sensor 8 at the target value Vop. 
Such an air-fuel ratio control process is carried out by 
determining a target air-fuel ratio in the exhaust gas sup- 
plied to the catalytic converter 4 according to a feedback 
control process in order to converge the output Vout of the 
0 2 sensor 8 to the target value Vop, and adjusting the 
amount of fuel to be supplied to the engine 1 according to 
a feedback control process in order to converge the output 
KACT (a detected value of the air-fuel ratio) of the wide- 
range air-fuel ratio sensor 9 to the target air-fuel ratio. 
Specific details of the air-fuel ratio control process car- 
ried by the air-fuel ratio control means 17 do not consti- 
tute an essential feature of the present invention, and 
will not be described below. The air-fuel ratio control 
process carried by the air-fuel ratio control means 17 is 
carried out as described in paragraphs [0071] - [0362] in 
the specification of Japanese laid-open patent publication 
No. 11-324767 or U.S. patent No. 6,188,953, for example. 

The output characteristics of the 0 2 sensor 8 
change depending on the temperature of the active element 
10 thereof. In FIG. 3, the solid-line curve "a", a broken- 
line curve "b", a dot-and-dash-line curve ff b", and a two- 
dot-and-dash-line curve "d" represent the output character- 
istics of the 0 2 sensor 8 when the active element 10 has 
temperatures of 800°C, 750°C, 700°C, and 600°C, respec- 
tively. As can be seen from FIG . 3, if the temperature of 
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the active element 10 changes in a temperature range lower 
than 750°C, then the gradient (sensitivity) of a change in 
the output Vout of the 0 2 sensor 8 in the vicinity of the 
stoichiometric air-fuel ratio (the high-sensitivity air- 
fuel ratio range A) and the level of the output Vout at 
air-fuel ratios richer than the high-sensitivity air-fuel 
ratio range A tend to change. If the temperature of the 
active element 10 is 750 °C or higher, then a change in the 
output characteristics of the 0 2 sensor 8 with respect to a 
change in the temperature of the active element 10 is so 
small that the output characteristics of the 0 2 sensor 8 
are substantially constant. 

Since the output characteristics of the 0 2 sen- 
sor 8 change depending on the temperature of the active 
element 10 as described above, the control properties (sta- 
bility and quick response) of the air-fuel ratio control 
means 17 are likely to be lowered depending on the tempera- 
ture of the active element 10. This is because in control- 
ling the air-fuel ratio in the exhaust gas in order to keep 
the output Vout of the 0 2 sensor 8 at the target value Vop, 
the output characteristics of the 0 2 sensor 8 in the vicin- 
ity of the stoichiometric air-fuel ratio, i.e., the output 
characteristics of the 0 2 sensor 8 in the high-sensitivity 
air-fuel ratio range A, are liable to greatly affect those 
control properties. The target value Vop for the output 
Vout of the 0 2 sensor 8 to keep well the ability of the 
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catalyst 7 of the catalytic converter 4 to purify the ex- 
haust gas also changes depending on the temperature of the 
active element 10 in a temperature range lower than 750 °C. 
Therefore, it is preferable to keep the temperature of the 
active element 10 of the 0 2 sensor 8 basically at a con- 
stant level for the purpose of well controlling the air-air 
ratio with the air-fuel ratio control means 17, i.e., con- 
trolling the output Vout of the 0 2 sensor 8 at the target 
value Vop, and achieving a good purifying ability of the 
catalytic converter 4 . 

If the temperature of the active element 10 of 
the 0 2 sensor 8 is 750 °C or higher, then the output charac- 
teristics of the 0 2 sensor 8 are substantially constant and 
stable. According to the inventors' knowledge, if the tem- 
perature of the active element 10 is kept at a temperature 
equal or higher than 750°C, e.g., 800°C, then the target 
value Vop for the output Vout of the 0 2 sensor 8 to keep 
well the ability of the catalyst 7 of the catalytic con- 
verter 4 to purify the exhaust gas is present in an area 
denoted by Y on the curve "a" in FIG. 3, i.e., an inflec- 
tion point Y where the gradient of the curve "a" represent- 
ing the output characteristics of the 0 2 sensor 8 switches 
from a larger value to a smaller value as the air-fuel ra- 
tio becomes richer. At this time, the air-fuel ratio can 
be controlled to keep the output Vout of the 0 2 sensor 8 at 
the target value Vop. The reason for the above air-fuel 
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- f - ■ appears to be that the sensitivity of the out- 

er' \z Vout he 0 2 sensor 8 to the air-fuel ratio at the 

inflection ;. >int Y is neither excessively high nor small, 
but is appropriate. 

According to the present embodiment, the sensor 
temperature control means 18 controls the ceramic heater 13 
to keep the temperature of the active element 10 of the 0 2 
sensor 8 at a desired temperature which is basically equal 
to or higher than 750°C, e.g., 800°C. A control process 
carried out by the sensor temperature control means 18 will 
be described below. 

As shown in FIG. 4, the sensor temperature con- 
trol means 18 has as its major functions an exhaust tem- 
perature observer 19 for sequentially estimating an exhaust 
gas temperature Tgd in the exhaust passage 3 near the 0 2 
sensor 8, i.e., at an intermediate portion of the exhaust 
pipe 6b, an element temperature observer 20 (temperature 
estimating means) for estimating the temperature T 02 of the 
active element 10 of the 0 2 sensor 8 and the temperature 
Tht of the ceramic heater 13 using the estimated value of 
the exhaust gas temperature Tgd, a target value setting 
means 21 for setting a target value R (target temperature) 
for the temperature of the active element 10, and a heater 
controller 22 (heater control means) for controlling ener- 
gization of the ceramic heater 13, i.e., controlling the 
electric energy supplied to the ceramic heater 13, using 
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the estimated values of the temperature T02 of the active 
element 10 and the temperature Tht of the ceramic heater 
13, the target value R, and the estimated value of the ex- 
haust gas temperature Tgd. The exhaust temperature ob- 
server 19 is supplied with detected data of the rotational 
speed NE of the engine 1, the intake pressure PB, and the 
atmospheric temperature TA in order to estimate the exhaust 
gas temperature Tgd. The heater controller 22 and the ele- 
ment temperature observer 20 are supplied with the detected 
data of the battery voltage VB for performing their proc- 
essing sequences as described later* 

In the present embodiment , the ceramic heater 13 
is controlled for its energization (PWM control) by giving 
a pulsed voltage to a heater energization circuit (not 
shown) . The amount of electric energy supplied to the ce- 
ramic heater 13 can be controlled by adjusting the duty cy- 
cle DUT of the pulsed voltage (the ratio of the pulse dura- 
tion to one period of the pulsed voltage) . The heater con- 
troller 22 sequentially determines the duty cycle DUT of 
the pulsed voltage applied to the heater energization cir- 
cuit as a control input (manipulated variable) for control- 
ling the ceramic heater 13, and adjusts the duty cycle DUT 
to control the amount of electric energy supplied to the 
ceramic heater 13 and hence the amount of heat generated by 
the ceramic heater 13. The duty cycle DUT generated by the 
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heater controller 22 is also used In the processing se- 
quence of the element temperature observer 20. 

According to the present embodiment, the portion 
of the exhaust passage 3 which extends from the exhaust 
port 2 of the engine 1 to the position where the 0 2 sensor 
8 is located, i.e., the exhaust passage 3 upstream of the 
0 2 sensor 8, is divided into a plurality of (four in the 
present embodiment) partial exhaust passageways 3a, 3b, 3c, 
3d along the direction in which the exhaust passage 3 ex- 
tends, i.e., the direction in which the exhaust gas flows. 
The exhaust temperature observer 19 estimates, in a prede- 
termined cycle time (period), the temperature of the ex- 
haust gas at the exhaust port 2 (the inlet of the exhaust 
passage 3) and the temperatures of the exhaust gas in the 
respective partial exhaust passageways 3a, 3b, 3c, 3d, or 
specifically, the temperatures of the exhaust gas in the 
downstream ends of the respective partial exhaust passage- 
ways 3a, 3b, 3c, 3d, successively in the downstream direc- 
tion. Of the partial exhaust passageways 3a, 3b, 3c, 3d, 
the partial exhaust passageways 3a, 3b are two partial ex- 
haust passageways divided from the exhaust passage 3 up- 
stream of the catalytic converter 4, i.e., the exhaust pas- 
sage defined by the exhaust pipe 6a, the partial exhaust 
passageway 3c is a partial exhaust passageway extending 
from the inlet to outlet of the catalytic converter 4, 
i.e., the exhaust passage defined in the catalyst 7 in the 
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catalytic converter 4, and the partial exhaust passageway 
3d is a partial exhaust passageway extending from the out- 
let of the catalytic converter 4 to the position where the 
0 2 sensor 8 is located, i.e., the exhaust pipe 6b. The ex- 
haust temperature observer 19 has its algorithm constructed 
as follows: 

The temperature of the exhaust gas at the ex- 
haust port 2 of the engine 1 basically depends on the rota- 
tional speed NE and the intake pressure PB of the engine 1 
while the engine 1 is operating in a steady state in which 
the rotational speed NE and the intake pressure PB are kept 
constant. Therefore, the temperature of the exhaust gas at 
the exhaust port 2 can basically be estimated from detected 
values of the rotational speed NE and the intake pressure 
PB, which serve as parameters indicative of the operating 
state of the engine 1, based on a predetermined map which 
has been established by way of experimentation, for exam- 
ple. If the operating state (the rotational speed NE and 
the intake pressure PB) of the engine 1 varies, then the 
temperature of the exhaust gas at the exhaust port 2 suf- 
fers a time lag or delay in the response to the exhaust gas 
temperature determined by the map (hereinafter referred to 
as "basic exhaust gas temperature TMAP ( NE , PB ) " ) . 

According to the present embodiment, the exhaust 
temperature observer 19 determines, in a predetermined cy- 
cle time (processing period), the basic exhaust gas tem- 
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perature TMAP ( NE , PB ) from the detected values (latest de- 
tected values) of the rotational speed NE and the intake 
pressure PB of the engine 1 based on the map, and thereaf- 
ter sequentially estimates an exhaust gas temperature Texg 
at the exhaust port 2 as a value which follows, with a time 
lag of first order, the basic exhaust gas temperature 
TMAP ( NE , PB ) as expressed by the following equation (1): 
Texg(k) = ( 1-Ktex) • Texg ( k- 1 ) +Kt ex* TMAP (NE,PB) 

•••(1) 

where k represents the ordinal number of a processing pe- 
riod of the exhaust temperature observer 19, and Ktex a co- 
efficient (lag coefficient) predetermined by way of experi- 
mentation or the like (0 < Ktex < 1). In the present em- 
bodiment, the intake pressure PB of the engine 1 serves as 
a parameter representative of the amount of intake air in- 
troduced into the engine 1. Therefore, if a flow sensor is 
used for directly detecting the amount of intake air intro- 
duced into the engine 1, then the output of the flow sen- 
sor, i.e., a detected value of the amount of intake air, 
may be used instead of the detected value of the intake 
pressure PB. In the present embodiment, an initial value 
Texg(0) of the estimated value of the exhaust gas tempera- 
ture Texg is set to the atmospheric temperature TA detected 
by an atmospheric temperature sensor (not shown) or the en- 
gine temperature TW (the temperature of the coolant of the 
engine 1) detected by an engine temperature sensor (not 
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shown) when the engine 1 has started to operate (upon an 



engine startup) , as described later. 

Using the estimated value of the exhaust gas 
temperature Texg at the exhaust port 2, the temperatures of 
the exhaust gas in the respective partial exhaust passage- 
ways 3a, 3b, 3c, 3d are estimated as described below. For 
illustrative purpose, a general heat transfer that occurs 
when a fluid flows through a circular tube 23 (see FIG. 5) 
which extends in the direction of a Z-axis in the atmos- 
phere while exchanging heat with the tube wall of the cir- 
cular tube 23 will be described below. It is assumed that 
the fluid temperature Tg and the temperature Tw of the tube 
wall (hereinafter referred to as "circular tube temperature 
Tw") are functions Tg(t,z), Tw(t,z) of the time t and the 
position z in the direction of the Z-axis, the thermal con- 
ductivity of the tube wall of the circular tube 23 is infi- 
nite in the radial direction and nil in the direction of 
the Z-axis. It is also assumed that the heat transfer be- 
tween the fluid and the tube wall of the circular tube 23 
and the heat transfer between the tube wall of the circular 
tube 23 and the external atmosphere are proportional to 
their temperature differences according to the Newton law 
of cooling. At this time, the following equations (2-1), 
(2-2) are satisfied: 




• • • (2-1) 



Sw • pw • Cw • 



dTw 

dt 



= d-U-(Tg-Tw) + a2-U (TA-Tw)- • • (2-2) 
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where S g , p g , C g represent the density and specific heat of 
the fluid and the cross-sectional area of the fluid pas- 
sage, respectively, S w , p w , C w the density, specific heat, 
and cross-sectional area of the tube wall of the circular 
tube 23, respectively, V the speed of the fluid flowing 
through the circular tube 23, T A the atmospheric tempera- 
ture outside of the circular tube 23, U the inner circum- 
ferential length of the circular tube 23, <Xi the heat 
transfer coefficient between the fluid and the tube wall of 
the circular tube 23, and cc 2 the heat transfer coefficient 
between the tube wall of the circular tube 23 and the at- 
mosphere. It is assumed that the atmospheric temperature 
T A is kept constant around the circular tube 23. 

The above equations (2-1), (2-2) are modified 

into the following equations (3-1), (3-2): 
dTg dTg 

-^- = -V~ + a (Tw-Tg) - - . (3-1) 

9Tw 

=b-(Tg - Tw) + c-(TA - Tw) • - • (3-2) 

where a, b, c represent constants, a = a x • U/ (S g - p g - C g ) , b = 
cti-U/ (-Sw-pw-Cw) , c = a 2 -U/ (S w -p w -C w ) . 

The first term on the right side of the equation 
(3-1) is a shifting flow term representing a time-dependent 
rate of change of the fluid temperature Tg (a change in the 
temperature per unit time) depending on the temperature 
gradient in the flowing direction of the fluid and the 
speed of the fluid in a position z. The second term on the 
right side of the equation (3-1) is a heat transfer term 
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representing a time -dependent rate of change of the fluid 
temperature Tg (a change in the temperature per unit time) 
depending on the difference between the fluid temperature 
Tg and the circular tube temperature Tw in the position z, 
i.e., a time -dependent rate of change of the fluid tempera- 
ture Tg which is caused by the heat transfer between the 
fluid and the tube wall of the circular tube 23. There- 
fore, the equation (3-1) indicates that the time -dependent 
rate dTg/dt of change of the fluid temperature Tg in the 
position z depends on the temperature change component of 
the shifting flow term and the temperature change component 
of the heat transfer term, i.e., the sum of those tempera- 
ture change components. 

The first term on the right side of the equation 
(3-2) is a heat transfer term representing a time -dependent 
rate of change of the circular tube temperature Tw (a 
change in the temperature per unit time) depending on the 
difference between the circular tube temperature Tw and the 
fluid temperature Tg in the position z, i.e., a time- 
dependent rate of change of the circular tube temperature 
Tw which is caused by the heat transfer between the fluid 
and the tube wall of the circular tube 23 in the position 
z. The second term on the right side of the equation (3-2) 
is a heat radiation term representing a time -dependent rate 
of change of the circular tube temperature Tw (a change in 
the temperature per unit time) depending on the difference 
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between the circular tube temperature Tw and the atmos- 
pheric temperature Ta outside of the circular tube 23 in 
the position z, i.e., a time -dependent rate of change of 
the circular tube temperature Tw depending on the heat ra- 
diation from the tube wall of the circular tube 23 into the 
atmosphere in the position z. The equation (3-2) indicates 
that the time -dependent rate dTw/dt of change of the circu- 
lar tube temperature Tw in the position z depends on the 
temperature change component of the heat transfer term and 
the temperature change component of the heat radiation 
term, i.e., the sum of those temperature change components. 

According to the calculus of finite differences, 
the equations (3-1), (3-2) can be rewritten into the fol- 
lowing equations (4-1), (4-2): 

Tg(t+At,z) = Tg(t,z)-^7^ (Tg(t,z)-Tg(t,z- Az) ) 

Az •••(4-1) 

+ a- At (Tw(t,z)- Tg(t,z) ) 

Tw(t+At,z)=Tw(t,z)+b • At • (Tg(t,z)-Tw(t,z) ) 

•••(4-2) 

+c • At • (TA-Tw(t,z) ) ' 
The above equations (4-1), (4-2) indicate that 
if the fluid temperature Tg(t,z) and the circular tube tem- 
perature Tw(t,z) in the position z at the time t, and the 
fluid temperature Tg(t,z-Az) in a position z-Az which pre- 
cedes the position z (upstream thereof) at the time t are 
known, then the fluid temperature Tg(t+At,z) and the circu- 
lar tube temperature Tw(t+At,z) in the position z at a next 
time t+At can be determined, and that the fluid tempera- 
tures Tg and the circular tube temperatures Tw in succes- 
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sive positions z+Az, z+2Az, • • • can be determined by solv- 
ing the equations (4-1), (4-2) simultaneously in sequence 
for those positions. Specifically, if initial values of 
the fluid temperature Tg and the circular tube temperature 
Tw (initial values at t = 0) are given in the positions z, 
z+Az, z+2Az, and the fluid temperature Tg(t,0) at an 

origin (e.g., the inlet of the circular tube 23) in the di- 
rection of the Z-axis of the circular tube 23 is given (it 
is assumed that z-Az = 0), then the fluid temperatures Tg 
and the circular tube temperatures Tw in successive posi- 
tions z, z+Az, z+2Az, • • • at successive times t, t+At, 
t+2At, • • • can be calculated. 

The fluid temperature Tg(t,z) in the position z 
can be calculated by cumulatively adding (integrating), to 
the initial value Tg(0,z), the temperature change component 
depending on the fluid speed V and the temperature gradient 
in the position z (the temperature change component repre- 
sented by the second term of the equation (4-1)) and the 
temperature change component depending on the difference 
between the fluid temperature Tg and the circular tube tem- 
perature Tw in the position z (the temperature change com- 
ponent represented by the third term of the equation (4- 
1)), at each given time interval. The fluid temperatures 
in the other positions z+Az, z+2Az, • • • can similarly be 
calculated. The circular tube temperature Tw(t,z) in the 
position z can be calculated by cumulatively adding (inte- 
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grating), to the initial value Tw(0,z), the temperature 
change component depending on the difference between the 
fluid temperature Tg and the circular tube temperature Tw 
in the position z (the temperature change component repre- 
sented by the second term of the equation (4-2)) and the 
temperature change component depending on the difference 
between the circular tube temperature Tw and the atmos- 
pheric temperature Ta in the position z (the temperature 
change component represented by the third term of the equa- 
tion (4-2)), at each given time interval. 

In the present embodiment, the exhaust tempera- 
ture observer 19 uses the model equations (4-1), (4-2) and 
determines the temperatures of the exhaust gas in the re- 
spective partial exhaust passageways 3a, 3b, 3c, 3d as fol- 
lows : 

Of the partial exhaust passageways 3a, 3b, 3c, 
3d, each of the partial exhaust passageways 3a, 3b is de- 
fined by the exhaust pipe 6a. In order to estimate the 
temperatures of the exhaust gas in the partial exhaust pas- 
sageways 3a, 3b, the temperature changes depending on the 
speed of the exhaust gas and the temperature gradient 
thereof (the temperature gradient in the direction in which 
the exhaust gas flows), the heat transfer between the ex- 
haust gas and the exhaust pipe 6a, and the heat radiation 
from the exhaust pipe 6a into the atmosphere are taken into 



account In the same manner as described above with respect 
to the circular tube 23. 

An estimated value of the exhaust gas tempera- 
ture Tga in the partial exhaust passageway 3a and an esti- 
mated value of the temperature Twa (hereinafter referred to 
as "exhaust pipe temperature Twa") of the exhaust pipe 6a 
in the partial exhaust passageway 3a are determined by re- 
spective model equations (5-1), (5-2), shown below, in each 
cycle time of the processing sequence of the exhaust tem- 
perature observer 19. An estimated value of the exhaust 
gas temperature Tgb in the partial exhaust passageway 3b 
and an estimated value of the exhaust pipe temperature Twb 
in the partial exhaust passageway 3b are determined by re- 
spective model equations (6-1), (6-2), shown below, in each 
cycle time of the processing sequence of the exhaust tem- 
perature observer 19. More specifically, the exhaust gas 
temperature Tga and the exhaust pipe temperature Twa that 
are determined by the equations (5-1), (5-2) represent es- 
timated values of the temperatures in the vicinity of the 
downstream end of the partial exhaust passageway 3a. Like- 
wise, the exhaust gas temperature Tgb and the exhaust pipe 
temperature Twb that are determined by the equations (6-1), 
(6-2) represent estimated values of the temperatures in the 
vicinity of the downstream end of the partial exhaust pas- 
sageway 3b. 
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Tga(k+l) = Tga(k) - Vg — (Tga(k) - Texg(k) ) + Aadt (Twa(k) - Tga(k) ) 

LiSL 

• • • (5-1) 

Twa(k + 1) = Twa(k) + Badt(Tga(k) -Twa(k) ) + Ca-dt (TA(k) - Twa(k) ) 

• • • (5-2) 

dt 

Tgb(k + l)=Tgb(k)- Vg — (Tgb(k) -Tga(k) ) + Abdt (Twb(k) -Tgb(k) ) 

JjD 

• • • (6-1) 

Twb(k + 1) = Twb(k) + Bbdt(Tgb(k) - Twb(k) ) +Cb dt (TA(k) - Twb(k) ) 

• • • (6-2) 

In the equations (5-1), (5-2), (6-1), (6-2), dt 
represents the period (cycle time) of the processing se- 
quence of the exhaust temperature observer 19, and corre- 
sponds to At in the equations (4-1), (4-2). In the equa- 
tions (5-1), (6-1), La, Lb represent the respective lengths 
(fixed values) of the partial exhaust passageways 3a, 3b, 
and correspond to Az in the equation (4-1) . Aa, Ba, Ca in 
the equations (5-1), (5-2) and Ab, Bb, Cb in the equations 
(6-1), (6-2) represent model coefficients corresponding re- 
spectively to a, b, c in the equations (4-1), (4-2), and 
the values of those model coefficients are set (identified) 
in advance by way of experimentation or simulation. In the 
equations (5-1), (6-1), Vg represents a parameter (to be 
determined as described later on) indicative of the speed 
of the exhaust gas, and corresponds to V in the equation 
(4-1) . 

The exhaust gas temperature Texg(k) (the exhaust 
gas temperature at the exhaust port 2) which is required to 
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calculate a new estimated value Tga(k+1) of the exhaust gas 
temperature Tga according to the equation (5-1) is basi- 
cally of the latest value determined according to the equa- 
tion (1). Similarly, the exhaust gas temperature Tga(k) 
(the exhaust gas temperature in the partial exhaust pas- 
sageway 3a) which is required to calculate a new estimated 
value Tgb(k+1) of the exhaust gas temperature Tgb according 
to the equation (6-1) is basically of the latest value de- 
termined according to the equation (5-1). The atmospheric 
temperature Ta ( k) which is required in the calculation of 
the equations (5-2), (6-2) is of the latest value of the 
atmospheric temperature detected by an atmospheric tempera- 
ture sensor (in the present embodiment, a sensor on the en- 
gine 1 is used for this atmospheric temperature sensor) , 
not shown. In the present embodiment, the gas speed pa- 
rameter Vg which is required in the calculation of the 
equations (5-1), (6-1) is of a value which is calculated 
from latest detected values of the rotational speed NE and 
the intake pressure PB according to the following equation 
(7) : 

NE PB 

Vg = • . • ( 7 ) 

NEBASE PBBASE V ' 

where NEBASE, PBBASE represent a predetermined rotational 
speed and a predetermined intake pressure, which are set 
to, for example, the maximum rotational speed of the engine 
1 and 760 mmHg (« 101 kPa) , respectively. The gas speed 
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parameter Vg calculated according to the equation (7) is 
proportional to the speed of the exhaust gas, with Vg £ 1 . 

In the present embodiment, initial values 
Tga(O), Twa(O), Tgb(O), Twb(O) of the estimated values for 
the exhaust gas temperature Tga, the exhaust pipe tempera- 
ture Twa, the exhaust gas temperature Tgb, and the exhaust 
pipe temperature Twb are set to the atmospheric temperature 
TA which is detected by the atmospheric temperature sensor 
(not shown) or the engine temperature TW (the temperature 
of the coolant of the engine 1) detected by the engine tem- 
perature sensor (not shown) when the engine 1 has started 
to operated (upon an engine startup) . 

The partial exhaust passageway 3c is defined by 
the catalyst 7 in the catalytic converter 4. The catalyst 
7 generates heat by itself due to its own exhaust gas puri- 
fying action (specifically, an oxidizing/reducing action) , 
and the amount of heat (the amount of heat per unit time) 
generated by the catalyst 7 is substantially in proportion 
to the speed of the exhaust gas. This is because as the 
speed of the exhaust gas is higher, the exhaust gas compo- 
nents reacting with the catalyst 7 per unit time increase. 

According to the present embodiment, for esti- 
mating the exhaust gas temperature in the partial exhaust 
passageway 3c with high accuracy, the generation of heat by 
the catalyst 7 in the catalytic converter 4 as well as the 
temperature change depending on the speed and temperature 
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gradient of the exhaust gas, the heat transfer between the 
exhaust gas and the catalyst 7, and the heat radiation from 
the catalyst 7 into the atmosphere are taken into account. 

An estimated value of the exhaust gas tempera- 
ture Tgc in the partial exhaust passageway 3c and an esti- 
mated value of the temperature Twc (hereinafter referred to 
as "catalyst temperature Twc") of the catalyst 7 which de- 
fines the partial exhaust passageway 3c are determined by 
respective model equations (8-1), (8-2), shown below, in 
each cycle time of the processing sequence of the exhaust 
temperature observer 19. More specifically, the exhaust 
gas temperature Tgc and the catalyst temperature Twc that 
are determined by the equations (8-1), (8-2) represent es- 
timated values of the temperatures in the vicinity of the 
downstream end of the partial exhaust passageway 3c, i.e., 

in the vicinity of the outlet of the catalytic converter 4. 

dt 

Tgc(k + 1) =Tgc(k) - Vg— (Tgc(k) - Tgb(k) ) + Ac- dt(Twc(k) - Tgc(k) ) 

■Li C 

• • • (8-1) 

Twc(k + 1) = Twc(k) + Bcdt(Tgc(k) - Twc(k) ) 

+ Cc dt-(TA(k) -Twc(k) ) + Dcdt-Vg " " * {8 " 2) 

In the equation (8-1), Lc represents the length 
(fixed value) of the partial exhaust passageway 3c, and 
corresponds to Az in the equation (4-1) . Ac, Be, Cc in the 
equations (8-1), (8-2) represent model coefficients corre- 
sponding respectively to a, b, c in the equations (4-1), 
(4-2) , and the values of those model coefficients are set 
(identified) in advance by way of experimentation or simu- 
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lation. The fourth term on the right side of the equation 
(8-2) represents a temperature change component of the 
catalyst 7 in the catalytic converter 4 due to the heating 
of the catalyst 7 by Itself, i.e., the temperature change 
per period of the processing sequence of the exhaust tem- 
perature observer 19, and is proportional to the gas speed 
parameter Vg. As with Ac through Cc, Dc in the fourth term 
represents a model coefficient that is set (identified) in 
advance by way of experimentation or simulation. There- 
fore, the equation (8-2) corresponds to the combination of 
the right side of the equation (4-2) with a temperature 
change component due to the heating of a passage -defining 
member (the catalyst 7). 

dt, Vg in the equations (8-1), (8-2) have the 
same meanings and values as those in the equations (5-1) 
through (6-2). The value of Ta used in the calculation of 
the equation (8-2) is identical to those used in the equa- 
tion (5-2), (6-2). In the present embodiment, the initial 
values Tgc(O), Twc(O) of the exhaust gas temperature Tgc 
and the catalyst temperature Twc are equal to the detected 
value of the atmospheric temperature TA or the detected 
value of the engine temperature TW at the time the engine 1 
has started to operate, as with the equations (5-1) through 
(6-2) . 

The partial exhaust passageway 3d is defined by 
the exhaust pipe 6b similar to the exhaust pipe 6a which 
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define the partial exhaust passageways 3a, 3b. The exhaust 

gas temperature Tgd in the partial exhaust passageway 3d 

and the exhaust pipe temperature Twd of the exhaust pipe 

6b, or more specifically the temperature at the downstream 

end of the partial exhaust passageway 3d, are determined 

respectively by the following equations (9-1), (9-2) which 

are similar to the equations (5-1) through (6-2) : 

dt 

Tgd(k + 1) =Tgd(k) - Vg— (Tgd(k) - Tgc(k) ) + Ad-dt (Twd(k) - Tgd(k) ) 

JjQ 

• • • (9-1) 

Twd(k + l)=Twd(k) + Bddt(Tgd(k) -Twd(k) ) + Cd dt (TA(k) - Twd(k) ) 

• • • (9-2) 

In the equation (9-1), Ld represents the length 
(fixed value) of the partial exhaust passageway 3d, and 
corresponds to Az in the equation (4-1) . Ad, Bd, Cd in the 
equations (9-1), (9-2) represent model coefficients corre- 
sponding respectively to a, b, c in the equations (4-1), 
(4-2), and the values of those model coefficients are set 
(identified) in advance by way of experimentation or simu- 
lation. 

dt, Vg in the equations (9-1), (9-2) have the 
same meanings and values as those in the equations (5-1) 
through (6-2) . The value of Ta used in the calculation of 
the equation (9-2) is identical to those used in the equa- 
tion (5-2), (6-2), (8-2). The initial values Tgd(O), 
Twd(0) of the exhaust gas temperature Tgd and the catalyst 
temperature Twd are equal to the detected value of the at- 
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mospheric temperature TA or the detected value of the en- 
gine temperature TW at the time the engine 1 has started to 
operate, as with the equations (5-1) through (6-2). 

The processing sequence of the exhaust tempera- 
ture observer 19, as described above, determines estimated 
values of the exhaust gas temperatures Texe, Tga, Tgb, Tgc, 
Tgd in the exhaust port 2 of the engine 1 and the partial 
exhaust passageways 3a, 3b, 3c, 3d successively downstream 
in each cycle time. The estimated value of the exhaust gas 
temperature Tgd in the partial exhaust passageway 3d which 
is located most downstream corresponds to the temperature 
of the exhaust gas in the vicinity of the location of the 
0 2 sensor 8. The estimated value of the exhaust gas tem- 
perature Tgd is obtained as the estimated value of the ex- 
haust gas temperature in the vicinity of the location of 
the 0 2 sensor 8. In the present embodiment, the estimated 
value of the exhaust gas temperature Tgd corresponds to ex- 
haust gas temperature data according to the present inven- 
tion. 

The algorithm of the estimating process of the 
exhaust temperature observer 19 is shown in block form in 
FIG. 6. In FIG. 6, the model equation (1) is referred to 
as an exhaust port thermal model 24, the model equations 
(5-1), (5-2) and the model equations (6-1), (6-2) as pre- 
CAT exhaust system thermal models 25, 26, respectively, the 
model equations (8-1), (8-2) as an in-CAT exhaust system 
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thermal model 27, and the model equations (9-1) , (9-2) as a 
post-CAT exhaust system thermal model 28. As shown in FIG. 
6, each of the thermal models 24 through 28 is supplied 
with the detected values of the rotational speed NE and the 
intake pressure PB of the engine 1. The detected values of 
the rotational speed NE and the intake pressure PB which 
are supplied to the exhaust port thermal model 24 are used 
to determine the basic exhaust gas temperature TMAP, and 
the detected values of the rotational speed NE and the in- 
take pressure PB which are supplied to the exhaust system 
thermal models 25 through 28 are used to determine the 
value of the gas speed parameter Vg. Each of the thermal 
models 25 through 28 is also supplied with the detected 
value of the atmospheric temperature Ta. The pre-CAT ex- 
haust system thermal model 25, the pre-CAT exhaust system 
thermal model 2 6, the in-CAT exhaust system thermal model 
27, and the post-CAT exhaust system thermal model 28 are 
supplied with the estimated values of the exhaust gas tem- 
peratures Texg, Tga, Tgb, Tgc, respectively, which are out- 
putted from the higher-level thermal models 24, 25, 26, 27. 
The post-CAT exhaust system thermal model 28 eventually 
produces the estimated value of the exhaust gas temperature 
Tgd in the vicinity of the location of the 0 2 sensor 8. 

In the present embodiment, the detected value 
produced by the atmospheric temperature sensor on the en- 
gine 1 is used to estimate the temperatures of the passage- 
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defining members (the exhaust pipe 6a, the catalyst 7 in 
the catalytic converter 4, and the exhaust pipe 6b) which 
define the partial exhaust passageways 3a, 3b, 3c, 3d. 
However, the detected value produced by an atmospheric sen- 
sor which is disposed outside of the exhaust passage 3 may 
be used to estimate the temperatures of those passage- 
defining members* 

The element temperature observer 20 will be de- 
scribed below. The element temperature observer 20 esti- 
mates the temperature T Q2 of the active element 10 of the 0 2 
sensor 8 sequentially in given cycle times in view of the 
thermal transfer between the active element 10 and the ex- 
haust gas held in contact therewith, the heat radiation 
from the active element 10 into the air within the active 
element 10, and the thermal transfer between the active 
element 10 and the ceramic heater 13 (hereinafter referred 
to simply as "heater 13") which heats the active element 
10* The element temperature observer 20 also estimates the 
temperature Tht of the heater 13 in order to estimate the 
temperature T 02 of the active element 10. In estimating 
the temperature Tht of the heater 13, the element tempera- 
ture observer 20 takes into account the heat transfer be- 
tween the heater 13 and the active element 13 and the heat 
radiation from the heater 13 into the air within the active 
element 10, and also the heating of the heater 13 based on 
the electric energy supplied to the heater 13. The element 
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temperature observer 20 has an estimating algorithm for es- 
timating the temperature T 0 2 and the temperature Tht, which 
is constructed as follows: 

The element temperature observer 20 determines 
an estimated value of the temperature T Q2 of the active 
element 10 (hereinafter referred to as "element temperature 
T02") of the 0 2 sensor 8 and an estimated value of the tem- 
perature Tht of the heater 13 (hereinafter referred to as 
"heater temperature Tht") sequentially in given cycle times 
respectively according to the model equations (10-1) , (10- 
2) described below. The equation (10-1) is the equation of 
an element temperature model, and the equation (10-2) is 

the equation of a heater temperature model. 
T O2 (k + l) = 0J >2 (k) + Axdt(Tgd(k)-^ 2 (k)) 

+ Bxdt(Tht(k)-T 02 (k) ) - Ex • dt (T 02 (k) - TA'(k) ) 

• • • (10-1) 

Tht(k + 1) = Tht(k)-Cxdt(Tht(k) -T Q2 (k) ) 

VB ( k ) 2 

- Fx • dt • ( Tht ( k ) - TA '( k ) ) + Dx • dt • DUT ( k ) —f- 

NVB 

• • • (10-2) 

The equation (10-1) indicates that the tempera- 
ture change of the active element 10 in each cycle time de- 
pends on a temperature change component (the second term on 
the right side of the equation (10-1)) depending on the 
difference between the exhaust gas temperature Tgd in the 
vicinity of the location of the 0 2 sensor 8 (the exhaust 
gas temperature in the partial exhaust passageway 3d) and 
the element temperature T Q2 , i.e., a temperature change 
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component which is caused by the heat transfer between the 
active element 10 and the exhaust gas held in contact 
therewith, a temperature change component (the third term 
on the right hand of the equation (10-1)) depending on the 
difference between the element temperature T 02 and the 
heater temperature Tht, i.e., a temperature change compo- 
nent which is caused by the heat transfer between the ac- 
tive element 10 and the ceramic heater 13, and a tempera- 
ture change component (the fourth term on the right hand of 
the equation (10-1)) depending on the difference between 
the element temperature T 02 and the temperature TA' of the 
air in the active element 10, i.e., a temperature change 
component due to the heat radiation from the active element 
10 into the air therein, i.e., the sum of those temperature 
change components . 

The equation (10-2) indicates that the tempera- 
ture change of the heater 13 in each cycle time depends on 
a temperature change component (the second term on the 
right side of the equation (10-2)) depending on the differ- 
ence between the element temperature T 0 2 and the heater 
temperature Tht , i.e., a temperature change component which 
is caused by the heat transfer between the active element 
10 and the heater 13, a temperature change component (the 
third term on the right hand of the equation (10-2)) de- 
pending on the difference between the heater temperature 
Tht and the temperature TA f of the air in the active ele- 
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ment 10 , i.e., a temperature change component due to the 
heat radiation from the heater 13 into the air within the 
active element 10, and a temperature change component (the 
fourth term on the right side of the equation (10-2)) de- 
pending on the product of the duty cycle DUT (more accu- 
rately, the duty cycle DUT that is actually used for the 
heater controller 22 to control the energization of the 
heater 13) that is generated by the heat controller 22 as 
described later on and the square VB 2 of the battery volt- 
age VB, i.e., a temperature change component which is 
caused by the heating of the heater 13 based on the elec- 
tric energy supplied thereto, i.e., the sum of those tem- 
perature change components. 

In the equations (10-1), (10-2), Ax, Bx, Cx, Dx, 
Ex, Fx represent model coefficients whose values are set 
(identified) in advance by way of experimentation or simu- 
lation, and dt represents the period (cycle time) of the 
processing sequence of the element temperature observer 20. 
In the present embodiment, the period dt is set to the same 
value as the cycle time (represented by dt in the equations 
(5-1) through (9-2)) of the processing sequence of the ex- 
haust temperature observer 19. In the equation (10-2), NVB 
represents a predetermined reference value (e.g., 14 V) of 
the battery voltage VB. The reference value may basically 
be set as desired as a standard voltage (a voltage that can 
normally be employed) for the battery voltage VB. 
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The fourth term on the right side of the equa- 
tion (10-2) will supplementarily be described below. If 
the duty cycle of the PWM control for the heater 13 is con- 
stant and the resistance of the heater 13 upon its energi- 
zation is constant, then the electric power supplied to the 
heater 13 is proportional to the square of the voltage ap- 
plied to the heater 13, and the applied voltage is propor- 
tional to the battery voltage VB. The duty cycle DUT de- 
fines the period of time in which the heater 13 is ener- 
gized, per period of the pulse voltage for the PWM control. 
Therefore, the product of the duty cycle DUT and the square 
VB 2 of the battery voltage VB is proportional to the elec- 
tric power supplied to the heater 13. The battery voltage 
VB changes as an alternator for charging the battery is 
turned on and off, for example. In the equation (10-2), 
the duty cycle DUT and the square VB 2 of the battery volt- 
age VB are multiplied in order to obtain a temperature 
change component due to the heating of the heater 13 when 
it is supplied with the electric power. 

The duty cycle DUT(k) which is required in the 
calculation of the equation (10-2) is of the latest value 
of the duty cycle DUT that is actually used for the heater 
22 to control the energization of the heater 13 (PWM con- 
trol). In the present embodiment, the latest value of the 
atmospheric temperature TA detected by the atmospheric tem- 
perature sensor is substituted for the temperature TA' (k) 
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of the air in the active element 10 which is required in 
the calculation of the equations (10-1) , (10-2). In the 
present embodiment, therefore, TA' (k) = TA(k) . Further- 
more, the initial values T 0 2(0), Tht(0) of the element tem- 
perature T0 2 and the heater temperature Tht are equal to 
the detected value of the atmospheric temperature TA or the 
detected value of the engine temperature TW at the time the 
engine 1 has started to operate, as described later. 

The element temperature observer 20 sequentially 
calculates the estimated values of the element temperature 
T 0 2 and the heater temperature Tht according to the esti- 
mating algorithm described above. The estimated values of 
the element temperature T 0 2 and the heater temperature Tht 
correspond respectively to element temperature data and 
heater temperature data in the present invention. 

The heater controller 22 will be described be- 
low. Basically, the heater controller 22 sequentially gen- 
erates the duty cycle DUT as a control input (manipulated 
variable) for controlling the heater 13 according to an op- 
timum predictive control algorithm, and controls the elec- 
tric energy supplied to the heater 13 with the generated 
duty cycle DUT. 

According to the present embodiment, attention 
is paid to the difference between the element temperature 
T 0 2 and a target value therefor, a change per given time in 
the difference (corresponding to a rate of change of the 
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difference) , and a change per given time in the heater tem- 
perature Tht (corresponding to a rate of change of the 
heater temperature Tht), and model equations for an object 
to be controlled by the heater controller 22 are introduced 
using the above differences and changes as state quantities 
relative to the object to be controlled by the heater con- 
troller 22. The heater controller 22 has its algorithm 
(optimum predictive control algorithm) constructed as de- 
scribed below. In the present embodiment, a duty cycle 
(control input) generated by the optimum predictive control 
algorithm, to be described below, is more precisely a duty 
cycle required to control the element temperature T 0 2 at a 
target value providing the battery voltage VB is kept con- 
stant at the reference value NVB represented by the equa- 
tion (10-2) . The duty cycle generated by the optimum pre- 
dictive control algorithm which will be described below is 
referred to as a basic duty cycle SDUT or a basic control 
input SDUT. 

First, model equations for the object to be con- 
trolled by the heater controller 22 will be described be- 
low. Changes AT 0 2, ATht per given time in the element tem- 
perature T 0 2 and the heater temperature Tht are expressed 
by the following equations (11-1), (11-2) based on the re- 
spective model equations (10-1), (10-2) with respect to the 
element temperature observer 20. The equation (11-2) is 
derived from an equation produced by using VB(k) = NVB and 
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replacing DUT with SDUT in the equation (10-2) . For deriv- 
ing the equations (11-1), (11-2), the temperature TA ? in 
the active element 10 is kept constant, i.e., remains un- 
changed, i.e., TA' (k+1) = TA'(k), in at least one cycle 
time of the processing sequence of the heater controller 
22. 

A^ 2 (k +1 ) - AT G2 (k) + Ax- dt (ATgd(k) - AT 02 (k) ) 

+ Bx dt (ATht(k)- A^ 2 (k) )- Ex- dt • AT^ 2 (k) 
= (l-Axdt-Bxdt-Exdt)AT 02 (k) * * ' 

+ Ax • dt • ATgd ( k ) + Bx • dt • ATht ( k ) 

ATht(k + l) = ATht(k)- Cx- dt- (ATht (k) - A^ 2 (k) ) 

- Fx • dt • ATht ( k ) + Dx • dt • ASDUT ( k ) 
= (l-Cx-dt.Fx-dt) ATht(k) + Cx dt-AT 02 (k) * ' " (11 ~ 2) 

+ Dxdt-ASDUT(k) 

In the above equations (11-1), (11-2), AT 0 2(k) = 
T 02 (k+1) - T Q2 (k), ATht(k) = Tht(k+1) - Tht(k), ATgd(k) = 
Tgd(k+1) - Tgd(k), ASDUT ( k) = SDUT(k+l) - SDUT(k). 

A target value for the element temperature T 0 2 
is represented by R, and the difference e between the ele- 
ment temperature T Q2 and the target value R, i.e., the dif- 
ference in each cycle time (hereinafter referred to as 
"element temperature difference e"), is defined according 
to the following equation (12) : 

e(k) = T 02 (k) - R(k) 

• * ' (12) 

A change Ae in the element temperature differ- 
ence e in each cycle time (hereinafter referred to as "ele- 
ment temperature difference change Ae") is expressed by the 
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following equation (13) based on the above equations (11- 
1), (12): 

Ae(k+1)= AT 02 (k+l)-AR(k+l) 

= (1-Ax-dt-Bx-dt-Ex-dt) • Ae ( k) +Ax • dt • ATgd (k) 

+Bx-dt- ATht (k) - AR ( k+1 ) + ( 1-Ax • dt-Bx • dt-Ex • dt ) • AR ( k) 

• ' • (13) 

In the equation (13), Ae(k) = e(k+l) - e(k), 
AR(k) = R(k+1) - R(k). In deriving the equation (13), the 
equation AT 02 = Ae(k) + AR(k) (based on the equation (12)) 
is employed. 

The equation AT 02 = Ae(k) + AR ( k) is applied to 
the equation (11-2), and the resulting equation is modified 
into the following equation (14): 

ATht (k+1) =(l-Cx-dt-Fx-dt) • ATht (k) +Cx-dt • Ae (k) 
+Dx" dt • ASDUT ( k) +Cx-dt-AR(k) • • • (14) 

If a state quantity vector XO(k) = (e(k),Ae(k), 
ATht(k)) T (T represents a transposition) is introduced, 
then the following equation (15) is obtained from the equa- 
tions (14), (15) and the equation e(k+l) = e(k)+Ae(k): 
X0(k+1) = O-X0(k)+G-ASDUT(k)+Gd-ATgd(k)+Gr-R0(k+l) 

■ • • (15) 

where X0(k) = (e ( k) , Ae.( k) , ATht ( k) ) T , 
R0(k+1) = ( AR ( k+1 ) , AR ( k) ) T , 
G = (0,0,Dx-dt) T , 
Gd = (0,Ax-dt,0) T , 
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hi o 1 

O = lo 1- Axdt-Bxdt-Exdt Bx-dt I 

Lo Cxdt 1-Cxdt-FxdtJ 

To o 1 

Gr =1 -1 1- Ax dt-Bx dt-Ex dt I 
Lo Cxdt J 

In the equation (15), RO, G, Gd represent vectors defined 
in the above definition clause, and O, Gr represent ma- 
trixes defined in the above definition clause. 

The above equation (15) is a basic equation of 
the model of the object to be controlled by the heater con- 
troller 22. 

In the above description, the period of the con- 
trol process of the heater controller 22 is the same as the 
period dt of the processing sequences of the exhaust tem- 
perature observer 19 and the element temperature observer 
20. Therefore, the period dt is used in the vectors G, Gd 
and the matrixes <D, Gr in the equation (15) . It is pref- 
erable to carry out the processing sequences of the exhaust 
temperature observer 19 and the element temperature ob- 
server 20 in a relatively short period (e.g., a period of 
20 through 50 msec.) in order to increase the accuracy with 
which to estimate the temperatures. However, the period of 
the control process of the heater controller 22 may be 
longer than the period dt of the processing sequences of 
the exhaust temperature observer 19 and the element tem- 
perature observer 20 because the response speed of a change 
in the element temperature in response to the control input 



- 113 - 



(duty cycle DUT) is relatively low (several Hz in terms of 
frequencies) . According to an optimum predictive control 
process to be described later on, since future values of 
the target value R of the element temperature T 02 need to 
be stored and held for a certain time, the storage capacity 
of a memory for storing the target value R becomes large if 
the period of the control process of the heater controller 
22 is short. 

According to the present embodiment, the period 
(cycle time) of the control process of the heater control- 
ler 22 is set to a value dtc (e.g., 300 through 500 msec.) 
longer than the period dt of the processing sequences of 
the exhaust temperature observer 19 and the element tem- 
perature observer 20. 

In the present embodiment, the model equation of 
the object to be controlled by the heater controller 22 is 
rewritten from the equation (15) into the following equa- 
tion (16), using the period dtc of the control process of 
the heater controller 22: 

X0 (n+1) = <D-X0 (n) +G-ASDUT (n) +Gd- ATgd(n) +Gr-R0 (n+1) 

• • • (16) 

where X0 (n) = (e (n) , Ae (n) , ATht (n) ) T , 
R0(n+1) = (AR(n+l) , AR(n) ) T , 
G = (0,0,Dx-dtc) T , 
Gd = (0, Ax- cite, 0) T , 
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11 0 
<D= 0 1-Ax • dtc-Bx • dtc-Ex • dtc Bx • dtc 

0 Cx - dtc 1-Cx • dtc-Fx - dtc 

To o 1 

Gr = 1 -1 1- Ax- dtc-Bx- dtc- Ex- dtc I 
|_0 Cx-dtc J 

The equation (16) is a model equation of the ob- 
ject to be controlled which is actually used in the algo- 
rithm of the control process of the heater controller 22. 
In the equation (16) , n represents the ordinal number of 
the period dtc of the control process of the heater con- 
troller 22. 

Using the above model equation, the algorithm of 
the control process of the heater controller 22, i.e., the 
algorithm of the optimum predictive control process, is 
constructed as follows: It is assumed that the target 
value R of the element temperature T 0 2 is set for the fu- 
ture until after Mr steps (until after a multiple by Mr of 
the period dtc of the control process of the heater con- 
troller 22) , and the exhaust gas temperature Tgd which acts 
as a disturbance input is known in the future until after 
Md steps (until after a multiple by Md of the period dtc of 
the control process of the heater controller 22) . The 
value Mr will be referred to as a target value predicting 
time Mr, and the value Md as an exhaust gas temperature 
predicting time Md. These predicting times Mr, Md are rep- 
resented by integers whose unit is one period dtc of the 
control process of the heater controller 22. 
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A controller for generating a control input 
ASDUT for minimizing the value of an evaluating function JO 
according to the following equation (17) serves as an opti- 
mum predictive servo controller: 

CO 

J0= ]T [XO T (n) - QO • XO (n) +ASDUT T (n) • HO • ASDUT (n) ] • • • (17) 

n=-M+l 

where M represents a larger one of the target value pre- 
dicting time Mr and the exhaust gas temperature predicting 
time Md, i.e., M = max (Mr , Md) , and QO , HO weighted matrixes 
for adjusting the convergence of the state quantity vector 
XO and the power (size) of the control input ASDUT. QO 
represents a 3-row, 3-column diagonal matrix as XO is a cu- 
bic matrix, and HO is a Scalar quantity as ASDUT is a Sca- 
lar quantity. In the present embodiment, in order to re- 
duce the power consumption of the heater 13, QO is set to a 
unit matrix (a diagonal matrix whose all diagonal compo- 
nents are "1") and HO is set to a value (e.g., 1000) 
greater than the diagonal components of the matrix QO . The 
target value predicting time Mr is set to 20, for example, 
and the exhaust gas temperature predicting time Md is set 
to 10, for example, with the period of the control process 
of the heater controller 22 being in the range from 300 to 
500 msec. 

The control input ASDUT for minimizing the value 
of the evaluating function according to the equation (17) 
is expressed by the equation (18) given below. In the pre- 
sent embodiment, it is assumed that the exhaust gas tern- 
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perature Tgd is maintained at the present value in the fu- 
ture until after Md steps. 

Mr 

ASDUT(n)=FO • XO (n) [FrO (i) • R0(n+1) ] +Fdt • ATgd(n) • • • (18) 

i°0 

In the equation (18), FO in the first term on 
the right side represents a cubic row vector ( FsO, FeO, FxO) , 
FrO(i) (i = 1, 2, Mr) in the second term (the term of 

2) on the right side represent quadratic row vectors 
(Fr01(i), Fr02(i), and Fdt in the third term on the right 
side represents a Scalar quantity. They are expressed by 
the equations (19-1) through (19-3) given below. 

FO = (FsO, FeO, FxO) 

= - [HO + G T - P-GJ^-G 7 - P-O 

. . . (19-1) 

FrO(i) = (Fr01(i), Fr02(i)) (i = 1, 2, •••, Mr) 
= - [HO + G T -P-G]' 1 -G T - (^ T ) i_1 -P^Gr 

• • ■ (19-2) 

Md 

Fdt = X{-[HO + G T -P-G]" 1 -G T -(^Y-P Gd} • • - (19-3) 

1=0 

where P represents a matrix (a 3-row, 3-column matrix) sat- 
isfying the following Ricatti equation (20-1), and £ repre- 
sents a matrix (a 3-row, 3-column matrix) expressed by the 
following equation (20-2) : 
P = Q0 + <X> T -P-0> 

- O- P-G- [H0+G T - P-G] ' 1 -G T - P-O •••(20-1) 
£ = <X> +G- FO • • • (20-2) 

G, Gr, Gd, and <X> in the equations (19-1) 
through (19-3) and the equations (20-1), (20-2) are defined 
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in the definition clause for the equation (16) , and HO, QO 
in those equations represent weighted matrixes of the 
evaluating function JO according to the equation (17) (HO 
is a Scalar quantity) . 

The second term (the term of £) on the right 
side of the equation (18) is rewritten using the components 
of FrO, RO (see the definition clauses for the equations 
(19-2), (16)), and then modified into the following equa- 
tion (21) : 

Mr Mr 

2[Fr0(i)R0(n + l)]=2[Fr(i)AR(n + l)] • • • (21) 

where 

fFr02(l) :i = 0 

Fr ( i ) = I FrOl ( i ) + Fr02 ( i + 1 ) :i - 1 ,2 • -,Mr - 1 
|_Fr01(Mr) :i =Mr 

By putting the equation (21) into the equation 
(18) and rewriting the first term on the right side of the 
equation (18) using the components of FO, XO (see the defi- 
nition clauses for the equations (19-1), (16)), the equa- 
tion (18) is expressed by the following equation (22) : 
ASDUT(n) = FsOe(n) + FeOAe(n) + FxOATht(n) 

Mr 

+ £[Fr(i)AR(n +1) ] +Fdt- ATgd(n) 

1=0 

• - • (22) 

Since the basic control input SDUT(n) to be gen- 
erated by the heater controller 22 is represented by the 
sum of its initial value SDUT(O) and ASDUT(l), ASDUT(2), 
ASDUT(n) cumulatively added thereto, the following 
equation (23) is obtained from the above equation (22) : 
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SDUT(n)=FsO • ^ e(j)+FeO • e(n)+FxO • Tht (n) 

j-i 

Mr 

+]T [Fr (i) • R(n+i) ] +Fdt • Tgd(n) 

i=0 

-FeO • e (0) -FxO • Tht (0) 

Mr 

-]T [Fr (i) • R(0+i) ]-Fdt • Tgd (0)+SDUT (0) 

i-0 

• • • (23) 

By setting the initial value terms of the equa- 
tion (23), i.e., the sixth term (the term of FeO-e(O)) 
through the tenth term (SDUT(O)), to "0", the following 
equation (24) is obtained as an equation for calculating 
the basic control input SDUT(n) to be actually generated by 
the heater controller 22: 

n 

SDUT(n) = FsO ^e(j) + FeO-e(n) +FxO Tht(n) 

Mr 

+ £[Fr(i)R(n+l) ]+Fdt Tgd(n) 

i-0 

• • • (24) 

The equation (24) is a formula for calculating 
the basic control input SDUT(n) (basic duty cycle) for con- 
trolling the heater 13 with the heater controller 22 ac- 
cording to the optimum predictive control algorithm. The 
first through third terms (the term including £e(j) through 
the term including Tht(n)) of the equation (24) represent 
control input components (a feedback component which will 
hereinafter be referred to as "optimum F/B component 
Uopfb") depending on the element temperature difference e 
and the heater temperature Tht. Specifically, the first 
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and second terms represent a control input component de- 
pending on the element temperature difference e, and the 
third term represents a control input component depending 
on the heater temperature Tht. The fourth term (the term 
of £Fr (i) -R(n+1) ) on the right side of the equation (24) 
represents a control input component (a feed-forward compo- 
nent which will hereinafter be referred to as "optimum tar- 
get value F/F component Uopfr") depending on the target 
value R. The fifth term (the term including Tgd(n)) repre- 
sents a control input component (a feed-forward component 
which will hereinafter be referred to as "optimum distur- 
bance F/F component Uopfd") depending on the exhaust gas 
temperature Tgd (which functions as a disturbance on the 
object to be controlled) . 

The basic control input SDUT(n) (basic duty cy- 
cle) determined by the equation (24) is a control input 
(duty cycle) that is required to control the element tem- 
perature T 0 2 at the target value R in the case where the 
battery voltage VB is kept constant at the reference value 
SVB, as described above. Therefore, the heater controller 
22 sequentially calculates a control input DUT(n) which is 
capable of controlling the element temperature T 02 at the 
target value R independent of the battery voltage VB by 
correcting the basic control input SDUT(n), which has been 
calculated in each cycle time (period) of the control proc- 
ess according to the equation (24) , with the ratio 
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NVB 2 /VB(n) 2 of the square of the present value VB(n) of the 

battery voltage VB and the square of the reference value 

NVB, as indicated by the following equation (25) : 

NVB 2 

DUT(n)= VB(n)2 SDUT(n) • • • (25) 

The heater controller 22 limits the control in- 
put DUT(n), i.e., the duty cycle DUT(n), within a predeter- 
mined range (0 % ^ DUT(n) ^ 100 %), and then applies a 
pulse voltage having the duty cycle DUT(n) to a heater en- 
ergization circuit (not shown) to adjust the electric power 
supplied to the heater 13. 

The heater controller 22 which determines the 
duty cycle DUT as the control input according to the equa- 
tions (24), (25) is expressed in block form as shown in 
FIG. 7. 

In the present embodiment, since the exhaust gas 
temperature Tgd is maintained at the present value in the 
future until after Md steps, Fdt in the fifth term on the 
right side of the equation (24) is a Scalar quantity. If 
Tgd in each step in the future can be detected or esti- 
mated, then the control input DUT can be determined using 
those values of Tgd. In this case, Fdt is a vector com- 
prising elements ( (Md + 1) elements) within {} of the equa- 
tion (19-3) . More specifically, if the series-series data 
of the exhaust gas temperature Tgd from the present until 
after Md steps is represented by Tgd(n), Tgd(n+1), 
Tgd(n+Md), then the fifth term on the right side of the 
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equation (24) is expressed by the inner product (Scalar 
product) of the vector comprising elements ( (Md + 1) ele- 
ments) within {} of the equation (19-3) and a vector com- 
prising, as elements, the series-series data Tgd(n), 
Tgd(n+1), Tgd(n+Md) of the exhaust gas temperature 

Tgd. The inner product is in conformity with the fifth 
term on the right side of the equation (24) when Tgd(n) = 
Tgd(n+1) = ••• =Tgd(n+Md). 

FsO, FeO, FxO which are required to determine 
the control input DUT(n) according to the equation (24) are 
of values calculated in advance according to the equation 
(19-1). Fr(i) (i = 0, 1, Mr) is of values calculated 

in advance according to the equations (21), (19-2). Fdt is 
of a value calculated in advance according to the equation 
(19-3). These coefficients FsO, FeO, FxO, Fr(i), Fdt may 
not necessarily be of the values according to the defining 
equations, but may be of values adjusted by way of simula- 
tion or experimentation. Furthermore, the coefficients 
FsO, FeO, FxO, Fr(i), Fdt may be changed depending on the 
element temperature, the heater temperature, etc. 

The heater temperature Tht and the exhaust gas 
temperature Tgd which are required in the calculation of 
the equation (24) are of the latest estimated value of the 
heater temperature Tht determined by the element tempera- 
ture observer 20 and the latest estimated value of the ex- 
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haust gas temperature Tgd determined by the exhaust tem- 
perature observer 19. 

The element temperature difference e required in 
the calculation of the equation (24) is calculated from the 
latest estimated value of the element temperature T 02 de- 
termined by the element temperature observer 20 and the 
target value R which has been set in a cycle time prior to 
the target value predicting time Mr by the target value 
setting means 21. 

The target value setting means 21 basically sets 
a temperature (e.g., 800° C in the present embodiment) equal 
to or higher than 750° C at which the output characteristics 
of the 0 2 sensor 8 are stably good, as the target value R 
for the temperature of the active element 10 in the same 
cycle time as the cycle time (period) of the processing se- 
quence of the heater controller 22. In order to perform 
the processing sequence of the heater controller 22 accord- 
ing to the algorithm of the optimum predictive control 
process, the target value setting means 21 sets the target 
value R in each cycle time as a target value R(n+Mr) after 
the target value predicting time Mr from the present cycle 
time, and stores a series of target values R(n+Mr) for the 
target value predicting time Mr. Specifically, the target 
value setting means 21 stores Mr+1 target values R(n), 
R(n+1), R(n+Mr) while sequentially updating them. The 

target value R used to determine the element temperature 



- 123 - 



difference e that is required in the calculation of the 
equation (24) is the value R(n) set and stored by the tar- 
get value setting means 21 as described above in the cycle 
time prior to the target value predicting time Mr. The 
target values R(n), R(n+1), R(n+Mr) stored as de- 

scribed above are used to determine the value of the fourth 
term (the term of £ including R(n+i)) of the equation (24). 

If the target value R of the element temperature 
T 0 2 is set to a high temperature such as 800 °C from the 
start of operation of the engine 1, then the active element 
10 tends to be damaged due to stresses caused by quick 
heating if water is applied to the active element 10 of the 
0 2 sensor 8 when the engine 1 starts to operate. In the 
present invention, therefore, until a certain time (e.g., 
15 seconds) elapses from the start of operation of the en- 
gine 1, the target value setting means 21 sets the target 
value R of the element temperature T 0 2 to a temperature 
lower than 750°C, e.g., 600°C. In the present embodiment, 
if the atmospheric temperature TA is low (e.g., TA < 0°C) 
upon elapse of a predetermined period of time after the en- 
gine 1 has started to operate, the target value R for the 
element temperature T 0 2 is set to a temperature (750 °C < R 
< 800°C) slightly lower than the normal target value 
(800 °C) in order to prevent the heater 13 from being over- 
heated. 
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Overall operation of the apparatus, particu- 
larly, the sensor temperature control means 18, according 
to the present embodiment will be described below. 

When the engine 1 starts to operate (upon an en- 
gine startup) , the sensor temperature control means 18 sets 
initial values Texg(O), Tga(O), Tgb(O), Tgc(O), Tgd(O), 
Twa(O), Twb(O), Twd(O), T o2 (0), and Tht(O) of the estimated 
values of the exhaust gas temperatures Texg, Tga, Tgb, Tgc, 
Tgd, the exhaust pipe temperatures Twa, Twb, Twd, the cata- 
lyst temperature Twc, the element temperature T02, and the 
heater temperature Tht, as follows: In the present embodi- 
ment, while the engine 1 stops its operation, the shutdown 
time of the engine 1 is measured. When the engine 1 starts 
to operate, the sensor temperature control means 18 deter- 
mines whether the shutdown time immediately prior to the 
engine startup is in excess of a predetermined time (e.g., 
2 hours) or not. If the shutdown time > the predetermined 
time, then since the temperature in the exhaust passage 3 
and the tube wall thereof is considered to be substantially 
the same as the atmospheric temperature, the sensor tem- 
perature control means 18 sets the initial values Texg(O), 
Tga(O), Tgb(O), Tgc(O), Tgd(O), Twa(O), Twb(O), Twd(O), 
T O 2(0), and Tht(O) to the detected value of the atmospheric 
temperature TA upon the startup of the engine 1. If the 
shutdown time < the predetermined time, then since the tem- 
perature in the exhaust passage 3 and the tube wall thereof 
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is considered to be closer to the engine temperature TW 
(coolant temperature) than to the atmospheric temperature, 
because of the afterheat of the engine 1 after the engine 1 
has stopped its preceding operation, the sensor temperature 
control means 18 sets the initial values Texg(O), Tga(O), 
Tgb(O), Tgc(O), Tgd(O), Twa(O), Twb(O), Twd(O), TO2(0), and 
Tht(O) to the detected value of the engine temperature TW 
upon the startup of the engine 1. These initial values are 
thus set to a temperature close to the actual temperature. 

When the engine 1 starts to operate, the sensor 
temperature control means 18 executes a main routine shown 
in FIG. 8 in a predetermined cycle time. The period in 
which the main routine is executed is shorter than the pe- 
riod dt of the processing sequence of the exhaust tempera- 
ture observer 19 and the element temperature observer 20 
and hence shorter than the period dtc of the processing se- 
quence of the target value setting means 21 and the heater 
controller 22. 

The sensor temperature control means 18 acquires 
detected values of the rotational speed NE and the intake 
pressure PB of the engine 1, the atmospheric temperature Ta 
and the battery voltage VB in STEP1, and then determines 
the value of a countdown timer COPC for measuring the time 
dtc of one period of the processing sequence of the target 
value setting means 21 and the heater controller 22 in 
STEP2. The value of the countdown timer COPC has been ini- 
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tialized to "0" at the time when the engine 1 starts to op- 
erate . 

If COPC = 0, then the sensor temperature control 
means 18 newly sets the value of the countdown timer COPC 
to a timer setting time TM1 which corresponds to the period 
dtc of the control processes of the target value setting 
means 21 and the heater controller 22 in STEP3 . Thereaf- 
ter, the target value setting means 21 carries out a proc- 
ess of setting a target value R for the element temperature 
T02 of the 0 2 sensor 8 in STEP4, and the heater controller 
22 carries out a process of calculating a duty cycle DUT of 
the heater 13 in STEPS . If COPC * 0 in STEP2 , then the 
sensor temperature control means 18 counts down the value 
of the countdown timer COPC in STEP5, and skips the proc- 
essing in STEP4 and STEPS . Therefore, the processing in 
STEP4 and STEPS is carried out at the period dtc determined 
by the timer setting time TM1. 

The processing in STEP4 and STEPS is specifi- 
cally carried out as follows: First, the processing in 
STEP4 is carried out by the target value setting means 21 
as shown in FIG. 9. 

The target value setting means 21 compares the 
value of a parameter TSH representative of the time that 
has elapsed from the start of the engine 1 with a predeter- 
mined value XTM (e.g., 15 seconds) in STEP4-1. If TSH = 
XTM, i.e., if the engine 1 is in a state Immediately after 
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it has started to operate, then the target value setting 
means 21 sets the target value R for the element tempera- 
ture T02 to a low temperature (e.g., 600°C) in order to 
prevent damage to the active element 10 of the 0 2 sensor 8 
in STEP4-2. Specifically, the target value R that is set 
at this time is a target value R(n+Mr) after the target 
value predicting time Mr from the present. 

If TSH > XTM in STEP4-1, then the target value 
setting means 21 sets the target value R for the element 
temperature T 0 2 from the present detected value (acquired 
in STEP1 shown in FIG. 8) of the atmospheric temperature Ta 
based on a predetermined table in STEP4-3. The target 
value R that is set at this time is basically a predeter- 
mined value (800°C in the present embodiment) equal to or 
higher than 750 °C if the atmospheric temperature Ta is a 
normal temperature (e.g., Ta > 0°C) . When the atmospheric 
temperature Ta is low (e.g., Ta < 0°C) as when the engine 1 
is operating in a cold climate, if the target value R for 
the element temperature T 0 2 is a high temperature of 800 °C, 
the temperature of the heater 13 is liable to be exces- 
sively high. In the present embodiment, when the tempera- 
ture of the heater 13 becomes excessively high, the heater 
13 is forcibly de-energized by an overheating prevention 
process (described later on) to prevent itself from a fail- 
ure. In STEP4-3, according to the present embodiment, when 
the atmospheric temperature Ta is low (e.g., Ta < 0°C) , the 
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target value R for the element temperature T 02 is set to a 
value slightly lower than the normal value (e.g., 750°C < R 
< 800°C) . 

Specifically, as with the target value R set in 
STEP4-2, the target value R set in STEP4-3 is a target 
value R(n+Mr) after the target value predicting time Mr 
from the present. 

After having set the target value R (= R(n+Mr) ) 
in STEP4-2 or STEP4-3, the target value setting means 21 
updates the values of Mr+1 buffers RBF(O), RBF(l) , 
RBF(Mr) for storing target values R for the target value 
predicting time Mr in STEP4-4, STEP4-5. The processing in 
STEP4 is now finished. 

In STEP4-4, specifically, the Mr buffers RBF(j) 
(j = 0/ 1/ "••/ Mr-1) are updated from the values of RBF(j) 
to the values of RBF ( j +1) , and the value held in the buffer 
RBF(O) so far is erased. In STEP4-5, the buffer RBF (Mr) is 
updated to the target value newly set in STEP4-2 or STEP4- 
3. The values of the buffers RBF (0) , RBF ( 1 ) , RBF (Mr) 

thus updated correspond respectively to R(n), R(n+1), 
R(n+Mr) in the fourth term of the equation (24) . The val- 
ues of the buffers RBF (0 ) , RBF ( 1 ) , RBF (Mr) have been 
initialized to a predetermined value (e.g., the target 
value set in STEP4-2) at the time the engine 1 has started 
to operate. 
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The processing in STEP5 is carried out by the 
heater controller 22 as shown in FIG. 10. The heater con- 
troller 22 calculates an element temperature difference 
e(n) = T 0 2(n) - RBF(O) between the present estimated value 
To2(n) of the element temperature T02 and the value of the 
buffer RBF(O) (= R(n)), i.e., the target value R set by the 
target value setting means 21 prior to the target value 
predicting time Mr in STEP5-1. 

Then, the heater controller 22 determines the 
values of flags F/A, F/B in STEP5-2. The flag F/A is set 
to "0" or "1" in a limiting process (described later on) 
for limiting the duty cycle DUT. The flag F/A which is set 
to "1" means that the duty cycle DUT is forcibly limited to 
a predetermined upper or lower limit value, and the flag 
F/B which is set to "0" means that the duty cycle DUT is 
not limited to the predetermined upper or lower limit value 
(the upper limit value > DUT > the lower limit value) . The 
flag F/B is set to "1" when the heater 13 is forcibly de- 
energized by the overheating prevention process. The flags 
F/A, F/B are initially set to "0". 

If F/A = F/B = 0 in STEP5-2, then the heater 
controller 22 adds the present value of £e(j) in the first 
term of the equation (24) to the difference e(n) calculated 
in STEP5-1 in STEP5-3 . In this manner, the difference e(n) 
is cumulatively added (integrated) in each cycle time dtc 
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of the processing sequence of the heater controller 22. 
The Initial value of £e(j) is "0". 

If F/A = 1 or F/B = 1 in STEPS -2 , then since the 
present value of the duty cycle DUT is not a normal value, 
the heater controller 22 skips the processing in STEPS -3 , 
but goes to STEPS -4 , holding the present value of £e(j). 

Then, the heater controller 22 calculates the 
equations (24), (25) using the present value (latest value) 
of the element temperature difference e(n) determined in 
STEP5-1 and the present accumulated value of 5>(;j), thus 
calculating the present value DUT(n) of the control input 
DUT for the heater 13 in STEPS -4 . Specifically, the heater 
controller 22 calculates the basic duty cycle SDUT(n) ac- 
cording to the equation (24) from the present value of the 
difference e(n) determined in STEP5-1, the present accumu- 
lated value £e(;j), the present estimated value Tht(n) of 
the heater temperature Tht, the present values (= R(n), 
R(n+1), R(n+Mr)) of the buffers RBF(O), RBF(l), 

RBF(Mr), the present estimated value Tgd(n) of the exhaust 
gas temperature Tgd (the exhaust gas temperature at the lo- 
cation of the 0 2 sensor 8), and the values of predetermined 
coefficients FsO, FeO, FxO, Fr(i) (i = 0, 1, Mr), Fdt . 

The heater controller 22 then calculates the duty cycle 
DUT(n) by correcting the basic duty cycle SDUT(n) using the 
present value (the latest value acquired in STEP1 shown in 
FIG- 8) of the battery voltage VB. The initial value 
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Tht(O) of the estimated value of the heater temperature Tht 
and the initial value Tgd(O) of the estimated value of the 
exhaust gas temperature Tgd, which are required for the 
first processing of STEP5-4 (at a stage where the process- 
ing sequences of the exhaust temperature observer 19 and 
the element temperature observer 20 are not executed) after 
the engine 1 has started to operate, are set to the atmos- 
pheric temperature TA or the engine temperature TW when the 
engine 1 starts to operate (an engine startup) . Those ini- 
tial values Tht(0), Tgd(0) are used in the calculation of 
the equation (24) . After the processing sequences of the 
exhaust temperature observer 19 and the element temperature 
observer 20 are executed, the latest estimated values of 
the estimated values that are determined in the processing 
sequences of the exhaust temperature observer 19 and the 
element temperature observer 20 are used in the calculation 
of the equation (24), 

Then, the heater controller 22 carries out a 
limiting process for limiting the duty cycle DUT (n) calcu- 
lated in STEP5-4 in STEP5-5 through STEP5-11. Specifi- 
cally, the heater controller 22 determines whether the duty 
cycle DUT (n) is smaller than a predetermined lower limit 
value (e.g., "0%") or not in STEP5-5. If DUT (n) < the 
lower limit value, then the heater controller 22 forcibly 
sets the value of DUT (n) to the lower limit value in STEPS- 
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6. Thereafter, the value of the flag F/A (the flag used in 
STEP5-2) is set to "1" in STEP5-7. 

If DUT(n) ^ the lower limit value, then the 
heater controller 22 determines whether the duty cycle 
DUT(n) is greater than a predetermined upper limit value 
(e.g., 100 %) or not in STEP5-8 . If DUT(n) > the upper 
limit value, then the heater controller 22 forcibly sets 
the value of DUT (n) to the upper limit value in STEP5-9 , 
Thereafter, the value of the flag F/A is set to "1" in 
STEP5-10. If the lower limit value ^ DUT (n) ^ the upper 
limit value, then the heater controller 22 holds the value 
of DUT(n), and sets the flag F/A to "0" in STEP5-11. The 
processing in STEPS is now finished. 

Control then returns to the main routine shown 
in FIG. 8. The sensor temperature control means 18 carries 
out the processing in STEP7 through STEP13. The processing 
in STEP7 through STEP13 represents a process of preventing 
the heater 13 from being overheated. In STEP7, the sensor 
temperature control means 18 determines whether or not the 
present estimated value (latest value) of the heater tem- 
perature Tht is equal to or higher than a predetermined up- 
per limit value THTLMT (e.g., 930°C). In the present em- 
bodiment, if Tht ^ THTLMT, the sensor temperature control 
means 18 forcibly de-energizes the heater 13 to prevent the 
heater 13 from being damaged. However, the estimated value 
of Tht may temporarily rise to a value equal to or higher 
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than the upper limit value THTLMT due to a disturbance or 
the like. According to the present embodiment, therefore, 
the sensor temperature control means 18 forcibly de- 
energizes the heater 13 if the state in which Tht = THTLMT 
has continued for a predetermined time (e.g., 3 seconds, 
hereinafter referred to as "heater OFF delay time"). 

If Tht < THTLMT in STEP 7 , then the sensor tem- 
perature control means 18 sets the value of a countdown 
timer TMHTOFF for measuring the heater OFF delay time to a 
predetermined value TM2 corresponding to the heater OFF de- 
lay time in STEP8 . Since the sensor temperature control 
means 18 does not forcibly de-energize the heater 13 at 
this time, the sensor temperature control means 18 sets the 
value of the flag F/B (the flag used in STEPS -2 shown in 
FIG. 10) to "0" in STEP9 . 

If Tht = THTLMT in STEP7 . then the sensor tem- 
perature control means 18 counts down the value of the 
countdown timer TMHTOFF by "1" in STEP10. Then, the sensor 
temperature control means 18 determines whether the value 
of the countdown timer TMHTOFF is "0" or not, i.e., whether 
the heater OFF delay time TM2 has elapsed with Tht = THTLMT 
or not in STEP11. 

If TMHTOFF * 0, then the sensor temperature con- 
trol means 18 sets the flag F/B to "0" in STEP9. If 
TMHTOFF = 0, then the sensor temperature control means 18 
forcibly sets the present value of the duty cycle DUT to 
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"0" in STEP12 , and then sets the value of the flag F/B to 
"1" in STEP13 . 

When the flag F/B is set to "0" in STEP9, the 
sensor temperature control means 18 applies a pulsed volt- 
age to the heater energization circuit according to the 
present value of the duty cycle DUT (the latest value cal- 
culated in STEPS), energizing the heater 13 with the elec- 
tric energy depending on the duty cycle DUT. When the 
value of the flag F/B is set to "1" in STEP12, the sensor 
temperature control means 18 does not apply a pulsed volt- 
age to the heater energization circuit, thus de-energizing 
the heater 13. 

After having thus executed the processing in 
STEP7 through STEP13, i.e., the process of preventing the 
heater 13 from being overheated, the sensor temperature 
control means 18 determines the value of a countdown timer 
COBS for measuring the time dt of one period of the proc- 
essing sequences of the exhaust temperature observer 19 and 
the element temperature observer 20 in STEP14. The value 
of the countdown timer COBS is initially set to "0" when 
the engine 1 has started to operate. 

If COBS = 0, then the sensor temperature control 
means 18 newly sets the value of COBS to a timer setting 
time TM3 (shorter than TM1 in STEP3) which corresponds to 
the period dt of the processing sequences of the exhaust 
temperature observer 19 and the element temperature ob- 
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server 20 in STEP15 . Then, the exhaust temperature ob- 
server 19 carries out a process of estimating the exhaust 
gas temperature Tgd (the exhaust gas temperature in the vi- 
cinity of the location of the 0 2 sensor 8), and the element 
temperature observer 20 carries out a process of estimating 
the element temperature T 02 (including a process of esti- 
mating the heater temperature Tht) in STEP16 . If COBS * 0 
in STEP14 , the sensor temperature control means 18 counts 
down the value of COBC in STEP17 , skipping the processing 
in STEP15 and STEP 16 . The processing in STEP16 is there- 
fore carried out at a period dt which is determined by the 
timer setting time TM3 . The main routine shown in FIG. 8 
is now finished. 

The processing in STEP16 is specifically carried 
out as shown in FIG. 11. The exhaust temperature observer 
19 successively carries out the processing in STEP16-1 
through STEP16-6 to determine an estimated value of the ex- 
haust gas temperature Tgd in the vicinity of the location 
of the 0 2 sensor 8. In STEP16-1, the exhaust temperature 
observer 19 determines a gas speed parameter Vg according 
to the equation (7) using the present detected values (the 
latest values acquired in STEP1) of the rotational speed NE 
and the intake pressure PB of the engine 1. The gas speed 
parameter Vg is forcibly set to Vg = 1 if the result calcu- 
lated by the equation (7) exceeds "1" due to an excessive 
rotational speed of the engine 1 . 
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Then, the exhaust temperature observer 19 calcu- 
lates an estimated value of the exhaust gas temperature 
Texg at the exhaust port 2 of the engine 1 according to the 
equation (1) in STEP16-2. Specifically, the exhaust tem- 
perature observer 19 determines a basic exhaust gas tem- 
perature TMAP (NE, PB) from the present detected values of 
the rotational speed NE and the intake pressure PB of the 
engine 1 based on the predetermined map, and thereafter 
calculates the right side of the equation (1) using the ba- 
sic exhaust gas temperature TMAP (NE, PB) , the present esti- 
mated value Texg(k-l) (determined in STEP16-2 in the pre- 
ceding cycle time) of the exhaust gas temperature Texg, and 
the value of a predetermined coefficient Ktex, thus calcu- 
lating a new estimated value Texg(k) of the exhaust gas 
temperature Texg. In the present embodiment, while the en- 
gine 1 is idling and also while the supply of fuel to the 
engine 1 is being cut off, the basic exhaust gas tempera- 
ture TMAP used in the calculation of the equation (1) is 
set to predetermined values corresponding to the respective 
engine operating states. When the engine 1 starts to oper- 
ate (upon an engine startup) , the atmospheric temperature 
Ta or the engine temperature TW detected at this time is 
set as an initial value Texg(O) of the estimated value of 
the exhaust gas temperature Texg. When the equation (1) is 
calculated for the first time after the engine 1 has 
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started to operate, the initial value Texg(O) is used as 
the value of Texg(k-l). 

Then, the exhaust temperature observer 19 calcu- 
lates an estimated value of the exhaust gas temperature Tga 
and an estimated value of the exhaust pipe temperature Twa 
in the partial exhaust passageway 3a according to the re- 
spective equations (5-1), (5-2) in STEP16-3. Specifically, 
the exhaust temperature observer 19 determines a new esti- 
mated value Tga(k+1) of the exhaust gas temperature Tga by 
calculating the right side of the equation (5-1) using the 
present estimated value Tga(k) (determined in STEP16-3 in 
the preceding cycle time) of the exhaust gas temperature 
Tga, the present estimated value (determined in STEP16-3 in 
the preceding cycle time) of the exhaust pipe temperature 
Twa, the present estimated value of the exhaust gas tem- 
perature Texg previously calculated in STEP16-2, the pre- 
sent value of the gas speed parameter Vg calculated in 
STEP16-1, the value of the predetermined model coefficient 
Aa, and the value of the period dt of the processing se- 
quence of the exhaust temperature observer 19. 

The exhaust temperature observer 19 calculates a 
new estimated value Twa(k+1) of the exhaust pipe tempera- 
ture Twa by calculating the right side of the equation (5- 
2) using the present estimated value Tga(k) (determined in 
STEP16-3 in the preceding cycle time) of the exhaust gas 
temperature Tga, the present estimated value (determined in 
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STEP16-3 in the preceding cycle time) of the exhaust pipe 
temperature Twa, the values of the predetermined model co- 
efficients Ba, Ca, and the value of the period dt of the 
processing sequence of the exhaust temperature observer 19. 

When the engine 1 starts to operate (upon an en- 
gine startup) , the atmospheric temperature Ta or the engine 
temperature TW detected at this time is set as initial val- 
ues Tga(O), Twa(O) of the estimated values of the exhaust 
gas temperature Tga and the exhaust pipe temperature Twa. 
When the equations (5-1), (5-2) are calculated for the 
first time after the engine 1 has started to operate, these 
initial values Tga(O), Twa(O) are used as the respective 
values of Tga(k-l), Twa(k-l). 

Then, the exhaust temperature observer 19 calcu- 
lates an estimated value of the exhaust gas temperature Tgb 
and an estimated value of the exhaust pipe temperature Twb 
in the partial exhaust passageway 3b according to the re- 
spective equations (6-1), (6-2) in STEP16-4. Specifically, 
the exhaust temperature observer 19 determines a new esti- 
mated value Tgb(k+1) of the exhaust gas temperature Tgb by 
calculating the right side of the equation (6-1) using the 
present estimated value Tgb(k) (the latest value determined 
in STEP16-4 in the preceding cycle time) of the exhaust gas 
temperature Tgb, the present estimated value (the latest 
value determined in STEP16-4 in the preceding cycle time) 
of the exhaust pipe temperature Twb, the present estimated 
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value of the exhaust gas temperature Tga previously calcu- 
lated in STEP16-3, the present value of the gas speed pa- 
rameter Vg calculated in STEP16-1, the value of the prede- 
termined model coefficient Ab, and the value of the period 
dt of the processing sequence of the exhaust temperature 
observer 19. 

The exhaust temperature observer 19 calculates a 
new estimated value Twb(k+1) of the exhaust pipe tempera- 
ture Twb by calculating the right side of the equation (6- 
2) using the present estimated value Tgb(k) (the latest 
value determined in STEP16-4 in the preceding cycle time) 
of the exhaust gas temperature Tgb, the present estimated 
value (the latest value determined in STEP16-4 in the pre- 
ceding cycle time) of the exhaust pipe temperature Twb # the 
values of the predetermined model coefficients Bb, Cb, and 
the value of the period dt of the processing sequence of 
the exhaust temperature observer 19. 

When the engine 1 starts to operate, the atmos- 
pheric temperature Ta or the engine temperature TW detected 
at this time is set as initial values Tgb(O), Twb(O) of the 
estimated values of the exhaust gas temperature Tgb and the 
exhaust pipe temperature Twb. When the equations (6-1), 
(6-2) are calculated for the first time after the engine 1 
has started to operate, these initial values Tgb(O), Twb(O) 
are used as the respective values of Tgb(k-l), Twb(k-l). 
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Then, the exhaust temperature observer 19 calcu- 
lates an estimated value of the exhaust gas temperature Tgc 
and an estimated value of the catalyst temperature Twc in 
the partial exhaust passageway 3c according to the respec- 
tive equations (8-1), (8-2) in STEP16-5. Specifically, the 
exhaust temperature observer 19 determines a new estimated 
value Tgc(k+1) of the exhaust gas temperature Tgc by calcu- 
lating the right side of the equation (8-1) using the pre- 
sent estimated value Tgc(k) (the latest value determined in 
STEP16-5 in the preceding cycle time) of the exhaust gas 
temperature Tgc, the present estimated value (the latest 
value determined in STEP16-5 in the preceding cycle time) 
of the catalyst temperature Twc, the present estimated 
value of the exhaust gas temperature Tgb previously calcu- 
lated in STEP16-4, the present value of the gas speed pa- 
rameter Vg calculated in STEP16-1, the value of the prede- 
termined model coefficient Ac, and the value of the period 
dt of the processing sequence of the exhaust temperature 
observer 19. 

The exhaust temperature observer 19 calculates a 
new estimated value Twc(k+1) of the catalyst temperature 
Twc by calculating the right side of the equation (8-2) us- 
ing the present estimated value Tgc(k) (the latest value 
determined in STEP16-5 in the preceding cycle time) of the 
exhaust gas temperature Tgc, the present estimated value 
(the latest value determined in STEP16-5 in the preceding 
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cycle time) of the catalyst temperature Twc, the present 
value of the gas speed parameter Vg calculated in STEP16-1, 
the values of the predetermined model coefficients Be, Cc, 
Dc, and the value of the period dt of the processing se- 
quence of the exhaust temperature observer 19. 

When the engine 1 starts to operate, the atmos- 
pheric temperature Ta or the engine temperature TW detected 
at this time is set as initial values Tgc(O), Twc(O) of the 
estimated values of the exhaust gas temperature Tgc and the 
catalyst temperature Twc. When the equations (8-1), (8-2) 
are calculated for the first time after the engine 1 has 
started to operate, these initial values Tgc(O), Twc(O) are 
used as the respective values of Tgc(k-l), Twc(k-l). 

Then, the exhaust temperature observer 19 calcu- 
lates an estimated value of the exhaust gas temperature Tgd 
and an estimated value of the exhaust pipe temperature Twd 
in the partial exhaust passageway 3d (near the location of 
the 0 2 sensor 8) according to the respective equations (9- 
1), (9-2) in STEP16-6. Specifically, the exhaust tempera- 
ture observer 19 determines a new estimated value Tgd(k+1) 
of the exhaust gas temperature Tgd by calculating the right 
side of the equation (9-1) using the present estimated 
value Tgd(k) (the latest value determined in STEP16-6 in 
the preceding cycle time) of the exhaust gas temperature 
Tgd, the present estimated value (the latest value deter- 
mined in STEP16-6 in the preceding cycle time) of the ex- 
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haust pipe temperature Twd, the present estimated value of 
the exhaust gas temperature Tgc previously calculated in 
STEP16-5, the present value of the gas speed parameter Vg 
calculated in STEP16-1, the value of the predetermined 
model coefficient Ad, and the value of the period dt of the 
processing sequence of the exhaust temperature observer 19. 

The exhaust temperature observer 19 calculates a 
new estimated value Twd(k+1) of the exhaust pipe tempera- 
ture Twd by calculating the right side of the equation (9- 
2) using the present estimated value Tgd(k) (the latest 
value determined in STEP16-6 in the preceding cycle time) 
of the exhaust gas temperature Tgd, the present estimated 
value (the latest value determined in STEP16-6 in the pre- 
ceding cycle time) of the exhaust pipe temperature Twd, the 
values of the predetermined model coefficients Bd, Cd, and 
the value of the period dt of the processing sequence of 
the exhaust temperature observer 19. 

When the engine 1 starts to operate, the atmos- 
pheric temperature Ta or the engine temperature TW detected 
at this time is set as initial values Tgd(O), Twd(O) of the 
estimated values of the exhaust gas temperature Tgd and the 
exhaust pipe temperature Twd. When the equations (9-1), 
(9-2) are calculated for the first time after the engine 1 
has started to operate, these initial values Tgd(O), Twd(O) 
are used as the respective values of Tgd(k-l), Twd(k-l). 
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Then, the element temperature observer 20 exe- 
cutes the processing in STEP16-7 to determine estimated 
values of the element temperature T02 of the 0 2 sensor 8 and 
the heater temperature Tht according to the equations (10- 
1), (10-2). Specif ically, the element temperature observer 
20 determines a new estimated value T 0 2(k+1) of the element 
temperature T Q2 by calculating the right side of the equa- 
tion (10-1) using the present estimated value T 0 2(k) (the 
latest value determined in STEP16-7 in the preceding cycle 
time) of the element temperature T 02 , the present estimated 
value Tht(k) (the latest value determined in STEP16-7 in 
the preceding cycle time) of the heater temperature Tht, 
the present estimated value Tgd(k) of the exhaust gas tem- 
perature Tgd previously calculated in STEP16-6, the present 
value TA(k) (the latest value acquired in STEP1 shown in 
FIG. 8) of the detected value of the atmospheric tempera- 
ture TA as the temperature TA' of the air in the active 
element 10, the values of the predetermined model coeffi- 
cients Ax, Bx, and the value of the period dt (= the period 
of the processing sequence of the exhaust temperature ob- 
server 19) of the processing sequence of the element tem- 
perature observer 20. 

Then, the element temperature observer 20 deter- 
mines a new estimated value Tht(k+1) of the heater tempera- 
ture Tht by calculating the right side of the equation (10- 
2) using the present estimated value T 0 2(k) (the latest 
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value determined in STEP16-7 in the preceding cycle time) 
of the element temperature T 0 2e the present estimated value 
Tht(k) (the latest value determined in STEP16-7 in the pre- 
ceding cycle time) of the heater temperature Tht, the pre- 
sent value TA(k) (the latest value acquired in STEP1 shown 
in FIG. 8) of the detected value of the atmospheric tem- 
perature TA as the temperature TA 1 of the air in the active 
element 10, the present value DUT(k) of the duty cycle DUT, 
the values of the predetermined model coefficients Cx, Dx, 
and the value of the period dt of the processing sequence 
of the element temperature observer 20. 

When the engine 1 starts to operate, the atmos- 
pheric temperature TA or the engine temperature TW detected 
at this time is set as initial values T O 2(0), Tht(0) of the 
estimated values of the element temperature T Q 2 and the 
heater temperature Tht. When the equations (10-1) , (10-2) 
are calculated for the first time after the engine 1 has 
started to operate, these initial values T O 2(0), Tht(0) are 
used as the respective values of T 0 2(k-1), Tht(k-l). The 
duty cycle DUT(k) used in the equation (10-2) is basically 
of the latest value determined by the heater controller 22 
in STEP5. However, if the value of the duty cycle DUT is 
limited in STEP12 to "0", i.e., to de-energize the heater 
13, then the limited value of the duty cycle DUT is used in 
the equation (10-2) . 
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The above processing sequence of the sensor tem- 
perature control means 18 controls the electric energy sup- 
plied to the heater 13 of the 0 2 sensor 8 in order to keep 
the element temperature T 0 2 of the 0 2 sensor 8 at the target 
value R. Except immediately after the engine 1 has started 
to operate and when the atmospheric temperature Ta is con- 
siderably low, the target value R is normally set to 800°C. 
As a result, the output characteristics of the 0 2 sensor 8 
can be maintained stably as the characteristics suitable 
for controlling the air-fuel ratio of the engine 1, i.e., 
for controlling the air-fuel ratio thereof for the cata- 
lytic converter 4 to perform a better exhaust purifying ca- 
pability, and the air-fuel ratio of the engine 1 can well 
be controlled to allow the catalytic converter 4 to perform 
a better exhaust purifying capability. 

According to the present embodiment, the duty 
cycle DUT as a control input for the heater 13 includes a 
control input component (the first and second terms of the 
equation (24) ) depending on the difference (element tem- 
perature difference) e between the estimated value of the 
element temperature T 02 and the target value R, and also 
includes a control input component (the third term of the 
equation (24) depending on the estimated value of the 
heater temperature Tht. Therefore, when the element tem- 
perature T02 varies with respect to the target value R, it 
is possible to converge the element temperature T Q2 
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smoothly to the target value R while preventing the duty 
cycle DUT of the heater 13 from being excessively changed. 
According to the present embodiment, since ATht of the 
heater temperature Tht is a state quantity of the model to 
be controlled, the control input component depending on the 
estimated value of the heater temperature Tht is meant to 
be a feedback component . 

According to the present embodiment, further- 
more, the duty cycle DUT also includes a control input com- 
ponent depending on the estimated value of the exhaust gas 
temperature Tgd which acts as a disturbant factor for vary- 
ing the element temperature T 02f i.e., the optimum distur- 
bance F/F component Uopfd. The coefficient Fdt relative to 
the optimum disturbance F/F component Uopfd is determined 
according to a predictive control algorithm on the assump- 
tion that the present exhaust gas temperature will continue 
until after the exhaust gas temperature predicting time Md. 
Consequently, it is possible to control the element tem- 
perature T02 at the target value R while suppressing varia- 
tions of the element temperature T G 2 due to variations of 
the exhaust gas temperature Tgd. In particular, as the co- 
efficient Fdt relative to the optimum disturbance F/F com- 
ponent Uopfd is determined according to a predictive con- 
trol algorithm, any variations of the element temperature 
T 0 2 due to variations of the exhaust gas temperature Tgd 
can be minimized. As a result, the stability of the proc- 
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ess of controlling the element temperature T02 at the tar- 
get value R is effectively increased, and hence the stabil- 
ity of the output characteristics of the 0 2 sensor 8 is 
also effectively increased. 

Furthermore, the control input DUT includes the 
control input component depending on the target value R for 
the element temperature T 02 , i.e., the optimum target value 
F/F component Uopfr. In addition, the optimum target value 
F/F component Uopfr is turned into a control input compo- 
nent depending on the target value R from the present until 
after the target value predicting time Mr by the predictive 
control algorithm. Therefore, when the target value R 
changes from a low temperature (600°C) immediately after 
the engine 1 has started to operate to a normal high tem- 
perature (750°C through 800°C) in particular, the control 
input DUT is prevented from becoming temporarily large ex- 
cessively, and the element temperature T 02 is prevented 
from overshooting with respect to the target value R. The 
stability of the output characteristics of the 0 2 sensor 8 
is also effectively increased. 

FIGS. 12 and 13 show the results of simulations 
of the present embodiment. FIG. 12 illustrates how the 
heater temperature Tht, the element temperature T 02 , and 
the duty cycle DUT change when the target value R changes 
from a low-temperature target value (600°C) to a high- 
temperature target value (800°C) at a time tl while a vehi- 
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cle carrying the engine 1 is running at a constant speed 
(in a steady operating state of the engine 1) . In FIG . 12, 
the solid-line curves e, g, i indicate how the heater tem- 
perature Tht, the element temperature T 0 2, and the duty cy- 
cle DUT, respectively, change when the heater 13 is con- 
trolled according to the present embodiment. The broken- 
line curves f, h, j indicate how the heater temperature 
Tht, the element temperature T 02 , and the duty cycle DUT, 
respectively, change in a comparative example . According 
to the comparative example, the duty cycle DUT is deter- 
mined using an equation similar to the equation (24) except 
that the optimum target value F/F component Uopfr (the 
fourth term of the equation (24)) depending on the target 
value R is removed from the equation (24). 

It can be seen from FIG. 24 that according to 
the present embodiment, since the duty cycle DUT includes 
the optimum target value F/F component Uopfr depending on 
the target value R, when the target value R switches, the 
duty cycle DUT is prevented from becoming excessively 
large, and transient overshooting of the heater temperature 
Tht and the element temperature T 0 2 is reduced. 

FIG. 13 shows how the element temperature T 0 2 
changes while a vehicle carrying the engine 1 is running at 
varying vehicle speeds as indicated by a lower curve in 
FIG. 13 (in various changing operating states of the engine 
1) . In FIG. 13, the solid-line curve p indicates how the 
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element temperature T02 changes when the heater 13 is con- 
trolled according to the present embodiment. The broken- 
line curve q indicates how the element temperature T 0 2 
changes in a comparative example. According to the com- 
parative example, the duty cycle DUT is determined using an 
equation similar to the equation (24) except that the opti- 
mum disturbance F/F component Uopfd (the fifth term of the 
equation (24)) depending on the exhaust gas temperature Tgd 
is removed from the equation (24) . 

It can be seen from FIG. 13 that according to 
the present embodiment, since the duty cycle DUT includes 
the optimum disturbance F/F component Uopfd depending on 
the exhaust gas temperature Tgd, a range of variations of 
the element temperature T 02 due to changes in the exhaust 
gas temperature Tgd is reduced. 

A second embodiment of the present invention 
will be described below with reference to FIG. 14. The 
second embodiment is partly different in arrangement or 
function from the first embodiment described above, and 
those structural or functional parts of the second embodi- 
ment which are identical to those of the first embodiment 
are denoted by identical reference characters, and will not 
be described in detail below. 

According to the present embodiment, as shown in 
the block diagram of FIG. 14, the sensor temperature con- 
trol means 18 of the control unit 16 shown in FIG. 1 com- 
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prises, as functional means, an exhaust temperature ob- 
server 19, an element temperature observer 20, a target 
value setting means 31, and a heater controller 32. The 
exhaust temperature observer 19 and the element temperature 
observer 20 are identical to those of the first embodiment. 
In the present embodiment, the target value setting means 
31 and the heater controller 32 have their processing peri- 
ods identical to the processing periods of the target value 
setting means 21 and the heater controller 22 according the 
first embodiment. 

The target value setting means 31 serves to set 
a target value R f for the heater temperature Tht of the 0 2 
sensor 8. According to the inventors' knowledge, the 
heater temperature Tht is relatively highly correlated to 
the element temperature T 0 2 and tends to be higher than the 
element temperature T 02 by a constant temperature in a 
steady state. According to the present embodiment, the 
target value setting means 31 sets, as the target value R 1 
for the heater temperature Tht, a value R + DR which is 
higher than the target value R for the element temperature 
T 02 that is set as described in the first embodiment (the 
target value R set by the processing sequence shown in FIG. 
9), by a predetermined value DR (e.g., 100°C) . As with the 
first embodiment, the target value R 1 that is set by the 
target value setting means 31 in each cycle time of its 
processing sequence is a target value after the target 
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value predicting time Mr, and the target value R 1 in the 
period of the target value predicting time Mr is sequen- 
tially updated and stored. 

The heater controller 32 sequentially generates 
the duty cycle DUT as a control input in order to keep the 
heater temperature Tht at the target value R' . In the pre- 
sent embodiment, as with the first embodiment, the heater 
controller 32 calculates a basic duty cycle SDUT according 
to an optimum predictive control algorithm, and corrects 
the basic duty cycle SDUT depending on the battery voltage 
VB according to the equation (25) thereby to generate the 
duty cycle DUT. 

More specifically, according to the present em- 
bodiment, attention is paid to the difference e f between 
the heater temperature Tht and a target value R' therefor, 
a change Ae 1 per given time in the difference e 1 (corre- 
sponding to a rate of change of the difference e'), and a 
change AT 02 per given time in the element temperature T 02 
(corresponding to a rate of change of the element tempera- 
ture T02) / and a model equation for an object to be con- 
trolled by the heater controller 32 is introduced using the 
above differences and changes as state quantities relative 
to the object to be controlled by the heater controller 32. 

If the difference e 1 (hereinafter referred to as 
"heater temperature difference e f ") is defined as e'(n) = 
Tht(n) - R'(n), then the model equation is given as the 
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following equation (26) based on the above equations (11- 
1) , (11-2) according to the same idea as with the first em- 
bodiment : 

Xl(n+1) = <I> * • XI (n) +G ' • ASDUT (n) +Gd ' • ATgd (n) +Gr ' • Rl (n+1 ) 

• • • (26) 

where XI (n) = (e • (n) , Ae ' (n) , AT G2 (n) ) T , 
Rl(n+1) = (AR' (n+1) ,AR' (n) ) T , 
G' = (0, Dx-dtc, 0) T , 
Gd' = (0,0,Ax-dtc) T , 

ri i o i 

0' = ! 0 l-Cxdtc-Fxdtc Cxdtc I 

|_0 Bxdtc 1- Ax dtc- Bx dtc-Ex dtcJ 

l~0 0 1 

Gr'=l-1 1- Cx dtc- Fx dtc I 
Lo Bxdtc J 

In the present embodiment, the basic control in- 
put SDUT (a control input at the time the battery voltage 
VB is equal to the reference value NVB) to be determined by 
the heater controller 22 is given by the equation (28) 
shown below as having integrated ASDUT (a control input on 
the model according to the equation (26) ) which minimizes 
an evaluating function Jl according to the following equa- 
tion (27) : 

Jl- 2 [Xl T (n) QO Xl(n) + ASDUT T (n) HO ASDUT(n) ] • • • (27) 

n=M+l 

where M = max (Mr, Md) 



- 153 - 

i 

n 

SDUT (n)=Fsl+£ e' (j)+Fel • e' (n)+Fxl • % 2 (n) 

j-i 

Mr 

+ J] [Fr* (i) • R' (n+i) ] +Fdt *• Tgd(n) 

i-0 

• ■ • (28) 

The coefficients Fsl, Fel, Fxl in the first 
through third terms, the coefficient Fr' (i) (i = 0, 1, 
Mr) in the fourth term, and the coefficient Fdt ' in the 
fifth term on the right side of the eguation (28) are coef- 
ficients given respectively by the following equations (2 9- 
1) through (29-3) : 

Fl ■ (Fsl, Fel, Fxl) 

= - [HO + G ' T • P ' •G , ]~ 1, G ,T *P' • <X> ' 

• • • (29-1) 

rFrl2(l) :i = 0 

Fr'(i) = l Frll(i) + Frl2(i +1) :i = 1 ,2 , -,Mr- 1 
[Frll(Mr) :i=Mr 

• • • (29-2) 

Md 

Fat , = ]^{-[HO + G ,T •P , •G , ]" 1 •G ,T •(£ ,T ) i •P , •Gd , } • • • (29-3) 

1=0 

where P' = QO + <D ,T -P''0' 

- <J>' • P ' • G ' • [HO + G ,T -P' -G' ]" 1 -G ,T - P' •<£* 
£ ' = <D • + g ' • Fl 
(Frll(i), Frl2(i)) 

= - [HO + G ,T - P' -G' ] _1 -G ,T - (£ ,T ) 1-1 • P' -Gr' 

(i = 1, 2, • • •, Mr) 
In- the present embodiment, the weighted matrixes 
QO, HO with respect to the evaluating function Jl, the tar- 
get value predicting time Mr, and the exhaust gas tempera- 
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ture predicting time Md are identical to those in the first 
embodiment. However, they may be set to values different 
from those in the first embodiment. The coefficients Fsl, 
Fel, Fxl, Fr'(i)/ Fdt ' in the equation (28) may not neces- 
sarily be of the values according to the defining equations 
(29-1) through (29-3), but may be of values adjusted by way 
of simulation or experimentation. Furthermore, the coeffi- 
cients Fsl, Fel, Fxl, Fr'(i), Fdt 1 may be changed depending 
on the element temperature, the heater temperature, etc. 
In the present embodiment, as with the first embodiment, 
the exhaust gas temperature Tgd is maintained at the pre- 
sent value in the future until after Md steps. However, if 
Tgd at each time in the future can be detected or esti- 
mated, then the control input DUT may be determined using 
those values (in this case, Fdt' is a vector). 

The above equation (28) is a formula for sequen- 
tially calculating a basic control input SDUT (n) with which 
the heater controller 32 controls the heater 13 in the pre- 
sent embodiment. Specifically, the heater controller 32 
sequentially calculates the basic control input SDUT(n) in 
each cycle time (period) of the control processing of the 
heater controller 32 according to the equation (28), and 
corrects the basic control input SDUT(n) depending on the 
battery voltage VB according to the equation (25), thereby 
determining the duty cycle DUT(n). The terms on the right 
side of the equation (28) have the same meanings as those 
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in the first embodiment. Specifically, the first through 
third terms (the term including Ee'Cj) through the term in- 
cluding T 0 2(n)) on the right side represent a control input 
component (a feedback component based on an optimum control 
algorithm) depending on the heater temperature difference 
e f and the element temperature T 0 2- More specifically, the 
first and second terms represent a control input component 
depending on the heater temperature difference e', and the 
third term represents a control input component depending 
on the element temperature T 0 2- The fourth term (the term 
of £Fr 1 ( i ) • R 1 (n+1 ) ) on the right side of the equation (28) 
and the fifth term (the term including Tgd(n)) on the right 
side thereof represent control input components (feed- 
forward components based on a predictive control algorithm) 
depending on the target value R and the exhaust gas tem- 
perature Tgd, respectively. 

As the element temperature T Q 2 and the exhaust 
gas temperature Tgd which are required to determine the ba- 
sic control input SDUT (n) according to the equation (28), 
there are employed, respectively, the latest value of the 
estimated value of the element temperature T 0 2 determined 
by the element temperature observer 20 and the latest value 
of the estimated value of the exhaust gas temperature Tgd 
determined by the exhaust temperature observer 19. 

The heater temperature difference e' required 
for the calculation according to the equation (28) is cal- 
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culated from the latest value of the estimated value of the 
heater temperature Tht determined by the element tempera- 
ture observer 20 and the target value R f that has been set 
in a cycle time before the target value predicting time Mr 
by the target value setting means 31. 

The other processing details than those de- 
scribed above are identical to those according to the first 
embodiment. In the present embodiment, the electric power 
supplied to the heater 13 of the 0 2 sensor 8 is controlled 
in order to maintain the heater temperature Tht of the 0 2 
sensor 8 at the target value R ' . In this case, except im- 
mediately after the engine 1 starts to operate or when the 
atmospheric temperature TA is considerably low (TA < 0°C) , 
the target value R ' is usually set to a temperature (900°C) 
which is higher than a preferred target temperature of 
800 °C for the active element 10 by a predetermined value DR 
(100°C in the present embodiment) . As a result, the tem- 
perature T 0 2 of the active element 10 of the 0 2 sensor 8 is 
indirectly controlled substantially at the temperature of 
800°C. Therefore, as with the first embodiment, the output 
characteristics of the 0 2 sensor 8 can stably be kept as 
characteristics suitable for controlling the air-fuel ratio 
of the engine 1 (for controlling the air-fuel ratio to keep 
a good purifying capability of the catalytic converter 4), 
and hence the air-fuel ratio is controlled well to reliably 
keep a good purifying capability of the catalytic converter 
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4 . During a predetermined period of time immediately after 
the engine 1 has started to operate, the target temperature 
R f for the heater 13 is set to a temperature (700°C) which 
is higher than a low temperature (600°C) as the target tem- 
perature R for the active element 10 than the predetermined 
value DR, for thereby preventing the active element 10 from 
being damaged by stresses due to abrupt heating. If the 
atmospheric temperature T A is low (TA < 0°C) , then inasmuch 
as the target value R for the active element 10 is set to a 
value in the range of 750°C ^ R < 800°C, the target value 
R 1 for the heater 13 is set to a value in the range of 
850°C ^ R ' < 900°C to prevent the heater 13 from being 
overheated. 

According to the present embodiment, the duty 
cycle DUT as a control input to the heater 13 includes the 
control input component (the first term (the term including 
£e'(j)) and the second term (the term including e f (n)) of 
the equation (28) ) depending on the difference between the 
estimated value of the heater temperature Tht an the target 
temperature R', and the control input component (the third 
term of the equation (28)) depending on the estimated value 
of the element temperature T 0 2/ as feedback components. In 
addition, according to the present embodiment, a predictive 
control algorithm is also applied, and the duty cycle DUT 
includes the control input component (a feed-forward compo- 
nent of the fifth term on the right side of the equation 
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(28)) depending on the exhaust gas temperature Tgd, and the 
control input component (a feed-forward component of the 
fourth term on the right side of the equation (28)) depend- 
ing on the target value R 1 . As a result, the present em- 
bodiment provides the same operation and advantages as with 
the first embodiment. Thus, the heater temperature Tht can 
reliably be controlled stably at the desired target value 
R', and the element temperature T 0 2 can be controlled 
stably at a desired temperature. 

A third embodiment of the present invention will 
be described below. The present embodiment is different 
from the first embodiment only as to the processing of the 
heater controller (specifically, the processing in STEP5-4 
shown in FIG. 10) , and those structural or functional parts 
of the third embodiment which are identical to those of the 
first embodiment are denoted by identical figures and ref- 
erence characters, and will not be described in detail be- 
low. 

According to the present embodiment, the heater 
controller 22 sequentially generates DUT as a control input 
according to an optimum predictive control algorithm, using 
the difference e (element temperature difference e) between 
the element temperature T 02 and the target value R, the 
element temperature T 0 2, and the heater temperature Tht as 
state quantities to be controlled. An algorithm of the 
processing sequence of the heater controller 22 is con- 
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structed as follows: The equations (11-1), (11-2) are 
brought together, providing the following equations (30-1), 
(30-2) : 

Xz(k+1) = Az-Xz(k)+Bz-ASDUT(k)+Ez-ATgd(k) •••(30-1) 

AT 02 (k) = Cz'Xz(k) •••(30-2) 

where Xz(k) = (AT 02 (k), Tht(k)) T 

[1- Axdt -Bxdt -Ex dt Bxdt ] 

AZ= [ Cx dt 1-Cx dt -Fx dtJ 

Bz = (0, Dx-dt) T 
Ez = (Ax-dt, 0) T 
Cz = (1, 0) 

From the equations (30-1), (30-2), there is ob- 
tained the following equation (31): 
AT 02 (k) = Cz-Xz(k+1) 

= Cz« Az'Xz(k)+Cz«Bz'ASDUT(k)+Cz'Ez-ATgd(k) 

•••(31) 

The equation (31) is the same as the equation 
(11-2), and is produced by rewriting the equation (11-2) 
using the matrix Az and the vectors Bz, Ez, Cz defined by 
the definition clauses of the equations (22-1), (22-2). 

In the present embodiment, a change (difference) 
Ae per each given time in the element temperature differ- 
ence e defined by the equation (12) is defined by Ae(k+1) = 
e(k+l) - e(k). The following equation (32) is obtained 
from this definition formula and the equation (31): 
e(k+l) = e(k) + Ae(k+1) 

- e(k) + (AT 02 (k+l) - AR(k+l)) 
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= e(k) + Cz-Az-Xz(k) + Cz • Bz • ASDUT ( k) 

+ Cz-Ez- ATgd(k) - AR(k+l) ---(32) 

The equation (32) and the above equation (30-1) 

are put together, providing the following equation (33) : 

X2(k+1) = (|>2 • X2 ( k) +G2 • ASDUT ( k) 

+ Gd2- ATgd(k) +Gr2- AR(k+l) •••(33) 

where X2(k) = (e(k), AT G2 (k), Tht(k)) T 

l~l 1- Ax-dt-Bx-dt-Ex-dt Bx-dt 1 

02 = I 0 1-Ax-dt-Bx-dt-Ex-dt Bx-dt I 

Lo Cxdt 1-Cxdt-FxdtJ 

G2 - (0, 0, Dx*dt) T 

Gd2 = (Ax-dt, Ax-dt, 0) T 

Gr2 = (-1, 0, 0) T 

The equation (33) is a basic formula of the 
model to be controlled by the heater controller 22 accord- 
ing to the present embodiment. In this model to be con- 
trolled, the state quantity to be controlled is a state 
quantity vector X2(k) = (e(k), AT 02 (k), ATht(k)) T comprising 
the element temperature difference e, the change AT C 2 per 
given time in the element temperature T Q 2r and the change 
ATht per given time in the heater temperature Tht . 

In the present embodiment, for the reasons which 
are the same as with the first embodiment, the cycle time 
of the cycle time of the processing sequence of the heater 
controller 22 is longer than the period of the processing 
sequence of the element temperature observer 20 and the ex- 
haust temperature observer 19. Consequently, the equation 
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of the model to be controlled which is actually used in the 
present embodiment is the following equation (34) using the 
period dtc of the processing sequence of the heater con- 
troller 22 and the ordinal number n of the processing pe- 
riod thereof: 

X2(n+1) = <|>2-X2 (n)+G2- ASDUT (n) 

+ Gd2-ATgd(n) + Gr2 - AR (n+1 ) —-(34) 

where X2(n) = (e(n), AT 0 2(n), Tht(n)) T 

fl 1- Ax- dtc- Bx-dtc- Ex- dtc Bx-dtc 1 

02= 0 1- Ax- dtc- Bx-dtc- Ex- dtc Bx-dtc I 

0 Cx-dtc 1-Cx-dtc-FxdtcJ 

G2 - (0, 0, Dx-dtc) T 

Gd2 = (Ax -dtc, Ax -dtc, 0) T 

Gr2 = (-1, 0, 0) T 

The basic control input SDUT (a control input at 
the time the battery voltage VB is equal to the reference 
value NVB) to be determined by the heater controller 22 ac- 
cording to the optimum predictive control algorithm based 
on the model equation (34) is given by the equation (36) 
shown below as having integrated ASDUT (a control input on 
the model according to the equation (34)) which minimizes 
an evaluating function J2 according to the following equa- 
tion (35) : 

OO 

J2= J] [X2 T {n) • QO • X2 (n) +ASDUT T (n) • HO • ASDUT (n) ] 

n — M+l 

• • • (35) 
where M = max (Mr, Md) 



- 162 - 



SDUT(n)=Fs2+2 e(j)+Fx2 - ^ 2 (n)+Fx3 • Tht(n) 
j-i 

Mr 

+ ]T [Fr2 (i) • R(n+I) ] +Fdt2 • Tgd(n) 

i-l 

• • • (36) 

The coefficients Fs2, Fx2, Fx3 in the first 
through third terms, the coefficient Fr2(i) (i = 0, 1, 
Mr) in the fourth term, and the coefficient Fdt2 in the 
fifth term on the right side of the equation (36) are coef- 
ficients given respectively by the following equations (37- 
1) through (37-3) : 

F2 ■ (Fs2, Fx2, Fx3) 

= - [HO + G2 T -P2-G2] -1 -G2 T -P2-02 •••(37-1) 
Fr2(i) = - [HO + G2 T -P2-G2] _1 -G2 T - (^2 T ) i_1 -P2-Gr2 

(i = 1, 2, • • • , Mr) • • • (37-2) 

Md 

Fdt2 = £{-[H0 + G2 T -P2-G2]- 1 -G2 T -(<-2 T ) 1 -P2-Gd2}- . . (37-3) 

i=0 

where P2 = Q0 + <X>2 T • P2 • <J>2 

- 02 • P2 • G2 • [HO + G2 T -P2-G2]" 1 -G2 T -P2-<D2 
C,2 = 0>2 + G2-F2 

In the present embodiment, the weighted matrixes 
Q0, HO with respect to the evaluating function J2, the tar- 
get value predicting time Mr, and the exhaust gas tempera- 
ture predicting time Md are identical to those in the first 
embodiment. However, they may be set to values different 
from those in the first embodiment. The coefficients Fs2, 
Fx2, Fx3, Fr2(i), Fdt2 in the equation (36) may not neces- 
sarily be of the values according to the defining equations 
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(37-1) through (237-3) , but may be of values adjusted by 
way of simulation or experimentation. Furthermore, the co- 
efficients Fs2, Fx2, Fx3, Fr2(i) f Fdt2 may be changed de- 
pending on the element temperature, the heater temperature, 
etc. In the present embodiment, as with the first embodi- 
ment, the exhaust gas temperature Tgd is maintained at the 
present value in the future until after Md steps. However, 
if Tgd at each time in the future can be detected or esti- 
mated, then the control input DUT may be determined using 
those values (in this case, Fdt2 is a vector) . 

The above equation (36) is a formula for sequen- 
tially calculating a basic control input SDUT (n) with which 
the heater controller 22 controls the heater 13 in the pre- 
sent embodiment. Specifically, the heater controller 22 
sequentially calculates the basic control input SDUT(n) in 
each cycle time (period) of the control processing of the 
heater controller 22 according to the equation (36), and 
corrects the basic control input SDUT(n) depending on the 
battery voltage VB according to the equation (25), thereby 
determining the duty cycle DUT(n), as with the first em- 
bodiment. The first through third terms (the term includ- 
ing £e(j) through the term including Tht (n) ) on the right 
side of the equation (36) represent a control input compo- 
nent (a feedback component based on an optimum control al- 
gorithm) depending on the heater temperature difference e, 
the element temperature T 0 2, and the heater temperature 
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Tht. More specif ically, the first term represents a con- 
trol input component depending on the heater temperature 
difference e, the second term represents a control input 
component depending on the element temperature T 02 , and the 
third term represents a control input component depending 
on the heater temperature Tht. The fourth term (the term 
of £Fr2 (i) * R (n+i) ) on the right side of the equation (36) 
and the fifth term (the term including Tgd(n)) on the right 
side thereof represent control input components (feed- 
forward components based on a predictive control algorithm) 
depending on the target value R and the exhaust gas tem- 
perature Tgd, respectively. 

As the element temperature T 0 2 and the heater 
temperature Tht which are required to determine the basic 
control input SDUT(n) according to the equation (36), there 
are employed, respectively, the latest value of the esti- 
mated value of the element temperature T 0 2 determined by 
the element temperature observer 20 and the latest value of 
the estimated value of the heater temperature Tht. As the 
exhaust gas temperature Tgd in the equation (36) , there is 
employed the latest value of the estimated value of the ex- 
haust gas temperature Tgd determined by the exhaust tem- 
perature observer 19. 

The heater temperature difference e required for 
the calculation according to the equation (36) is calcu- 
lated from the latest value of the estimated value of the 
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element temperature T02 determined by the element tempera- 
ture observer 20 and the target value R(n) that has been 
set in a cycle time before the target value predicting time 
Mr by the target value setting means 21. The fourth term 
on the right side of the equation (36) is calculated using 
time-series data R(n+1), R(n+2), R(n+Mr) of the target 

value R which have been set from a cycle time which is one 
cycle time after the cycle time before the target value 
predicting time Mr to the present time by the target value 
setting means 21. 

The other processing details than those de- 
scribed above are identical to those according to the first 
embodiment. More specifically, in the present embodiment, 
the heater controller 22 determines DUT(n) according to the 
equations (36), (37) in STEPS -4 shown in FIG. 10. The 
other processing details than the processing in STEPS -4 are 
exactly the same as those according to the first embodi- 
ment. In the present embodiment, as with the first embodi- 
ment, the element temperature T 02 is controlled at the tar- 
get value R. In this case, the basic control input SDUT 
does not include an input component that is directly pro- 
portional to the element temperature difference e, but in- 
cludes a component (the second term of the equation (36)) 
depending on the element temperature T 02 . unlike the first 
embodiment. Therefore, immediately after the engine 1 
starts to operate or when the target value R is changed by 
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the target value setting means 21 from a low- temperature 
target value (600°C) to a high -temperature target value 
(750°C through 800° C), as the element temperature T 02 ap- 
proaches the target value R # the element temperature T 0 2 is 
effectively prevented from overshooting with respect to the 
target value R. The advantages produced by the basic con- 
trol input SDUT and the duty cycle DUT which include input 
components depending on the heater temperature Tht, the 
target value R, and the exhaust gas temperature Tgd are the 
same as those according to the first embodiment. 

A fourth embodiment of the present invention 
will be described below. The present embodiment is differ- 
ent from the second embodiment only as to the processing of 
the heater controller, and those structural or functional 
parts of the fourth embodiment which are identical to those 
of the second embodiment (including those of the first em- 
bodiment which are incorporated in the second embodiment) 
are denoted by identical figures and reference characters, 
and will not be described in detail below. 

The present embodiment is different from the 
second embodiment only as to the processing sequence of the 
heater controller 32 (FIG. 14) described in the second em- 
bodiment, and determines the duty cycle DUT according to an 
optimum predictive control algorithm in the form described 
in the third embodiment. According to the present embodi- 
ment, specifically, the state quantity of an object to be 
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controlled by the heater controller 32 is a state quantity 
vector X3(n) = (e'(n), AT 0 2(n), ATht(n))T comprising the 
difference (heater temperature difference) e 1 between the 
heater temperature Tht and the target value R' , a change 
AT02 per given time in the element temperature T 0 2/ and a 
change ATht per given time in the heater temperature Tht. 
That is, the state quantity of the object to be controlled 
uses the heater temperature difference e'(n), rather than 
the heater temperature difference e(n) of the state quan- 
tity vector X2(n) = (e(n), AT 0 2(n), ATht(n))T in the third 
embodiment. If a change Ae ' per given time in the heater 
temperature difference e 1 is redefined as Ae'(n+1) = 
e 1 (n+1) -e 1 (n) , then a model to be controlled is expressed 
by the following equation (38) according to the same con- 
cept as with the third embodiment: 
X3(n+1) = <t>3 • X3 (n) +G3 • ASDUT (n) 

+ Gd3-ATgd(n)+Gr3-AR' (n+1) -—(38) 

where X3 (n) = (e'(n), AT 02 (n), Tht(n)) T 

Tl Cx • dt c 1 - Cx • dt c - Fx • dt cl 

0)3 = 1 0 1- Ax-dtc- Bxdtc-Exdtc Bx dtc 1 

[o Cxdtc 1-Cxdtc-FxdtcJ 

G3 = (Dx-dtc f 0, Dx;dtc) T 

Gd3 = (0, Ax-dtc, 0) T 

Gr3 = (-1, 0, 0) T 

The basic control input SDUT (a control input at 
the time the battery voltage VB is equal to the reference 
value NVB) to be determined by the heater controller 32 ac- 



- 168 - 



cording to the present embodiment is given by the equation 
(40) shown below as having integrated ASDUT (a control in- 
put on the model according to the equation (38)) which 
minimizes an evaluating function J3 according to the fol- 
lowing equation (39) : 

oo 

J3== 2 [X3 T (n) . QO • X3 (n) +ASDUT T (n) • HO - ASDUT (n) ] • • • (39) 

n=-M+l 

where M = max (Mr, Md) 

n 

SDUT(n)=Fs3+5] e' (j)+Fx4 • 1^ 2 (n)+Fx5 • Tht(n) 

•••(40) 

+ 2 [Fr3(i) • R' (n+i) ]+Fdt3 • Tgd(n) 

i-1 

The coefficients Fs3, Fx4, Fx5 in the first 
through third terms, the coefficient Fr3(i) (i = 0, 1, 
Mr) in the fourth term, and the coefficient Fdt3 in the 
fifth term on the right side of the equation (40) are coef- 
ficients given respectively by the following equations (41- 
1) through (41-3) : 

F3 = (Fs3, Fx4, Fx5) 

= - [HO + G3 T -P3-G3] _1 -G3 T -P3-03 •••(41-1) 
Fr3(i) = - [HO + G3 T -P3-G3]" 1 -G3 T - (<;3 T ) i - 1 -P3-Gr3 

(i = 1, 2, • • • , Mr) • • • (41-2) 

Md 

Fdt3 = X{-[H0 + G3 T P3 G3]- 1 G3 T (^3 T ) i P3 Gd3}- • • (41-3) 

1-0 

where P3 = QO + <t>3 T • P3 • 03 

- 03- P3-G3- [HO + G3 T -P3-G3]- 1 -G3 T -P3-0)3 
C3 = <J>3 + G3-F3 

In the present embodiment, the weighted matrixes 
QO, HO with respect to the evaluating function J3, the tar- 
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get value predicting time Mr, and the exhaust gas tempera- 
ture predicting time Md are identical to those in the first 
embodiment. However, they may be set to values different 
from those in the first embodiment. The coefficients Fs3, 
Fx4, Fx5, Fr3(i), Fdt3 in the equation (40) may not neces- 
sarily be of the values according to the defining equations 
(41-1) through (41-3) , but may be of values adjusted by way 
of simulation or experimentation. Furthermore, the coeffi- 
cients Fs3, Fx4, Fx5, Fr3 (i) , Fdt3 may be changed depending 
on the element temperature, the heater temperature, etc. 
In the present embodiment, as with the first embodiment, 
the exhaust gas temperature Tgd is maintained at the pre- 
sent value in the future until after Md steps. However, if 
Tgd at each time in the future can be detected or esti- 
mated, then the control input DUT may be determined using 
those values (in this case, Fdt3 is a vector) . 

The above equation (40) is a formula for sequen- 
tially calculating a basic control input SDUT(n) with which 
the heater controller 32 controls the heater 13 in the pre- 
sent embodiment. Specifically, the heater controller 32 
sequentially calculates the basic control input SDUT(n) in 
each cycle time (period) of the control processing of the 
heater controller 32 according to the equation (40), and 
corrects the basic control input SDUT(n) depending on the 
battery voltage VB according to the equation (25), thereby 
determining the duty cycle DUT(n), as with the first em- 
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bodiment. The first through third terms (the term includ- 
ing £e' (j) through the term including Tht (n) ) on the right 
side of the equation (40) represent a control input compo- 
nent (a feedback component based on an optimum control al- 
gorithm) depending on the heater temperature difference e f , 
the element temperature T 02 , and the heater temperature 
Tht. More specifically the first term represents a con- 
trol input component depending on the heater temperature 
difference e 1 , the second term represents a control input 
component depending on the element temperature T 02 , and the 
third term represents a control input component depending 
on the heater temperature Tht- The fourth term (the term 
of £ Fr 3 (i) . R 1 (n+i) ) on the right side of the equation (40) 
and the fifth term (the term including Tgd(n)) on the right 
side thereof represent control input components (feed- 
forward components based on a predictive control algorithm) 
depending on the target value R 1 and the exhaust gas tem- 
perature Tgd, respectively. 

As the element temperature T 0 2 and the heater 
temperature Tht which are required to determine the basic 
control input SDUT(n) according to the equation (40), there 
are employed, respectively, the latest value of the esti- 
mated value of the element temperature T 02 determined by 
the element temperature observer 20 and the latest value of 
the estimated value of the exhaust gas temperature Tgd de- 
termined by the exhaust temperature observer 19. The 
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heater temperature difference e' required for the calcula- 
tion according to the equation (40) is calculated from the 
latest value of the estimated value of the heater tempera- 
ture Tht determined by the element temperature observer 20 
and the target value R 1 (n) that has been set in a cycle 
time before the target value predicting time Mr by the tar- 
get value setting means 31. The fourth term on the right 
side of the equation (40) is calculated using time-series 
data R'(n+1), R 1 (n+2) , R 1 (n+Mr) of the target value R' 

which have been set from a cycle time which is one cycle 
time after the cycle time before the target value predict- 
ing time Mr to the present time by the target value setting 
means 31. 

The other processing details than those de- 
scribed above are identical to those according to the sec- 
ond embodiment. In the present embodiment , as with the 
second embodiment, the heater temperature Tht is controlled 
at the target value R', and the element temperature T 0 2 is 
controlled at a temperature corresponding to the target 
value R 1 . In this case, the basic control input SDUT does 
not include an input component that is directly propor- 
tional to the element temperature difference e 1 , but in- 
cludes a component (the third term of the equation (40)) 
depending on the heater temperature Tht, unlike the second 
embodiment. Therefore, immediately after the engine 1 
starts to operate or when the target value R 1 is changed by 
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the target value setting means 31 from a low-temperature 
target value (700 °C in the present embodiment) to a high- 
temperature target value (850°C through 900°C) , as the 
heater temperature Tht approaches the target value R 1 , the 
heater temperature Tht is effectively prevented from over- 
shooting with respect to the target value R' . The element 
temperature T 0 2 is thus allowed to converge smoothly to a 
temperature corresponding to the target value R' for the 
heater temperature Tht. The advantages produced by the ba- 
sic control input SDUT and the duty cycle DUT which include 
input components depending on the heater temperature T 02 , 
the target value R' , and the exhaust gas temperature Tgd 
are the same as those according to the second embodiment. 

In the first through fourth embodiments de- 
scribed above, the exhaust gas temperature Tgd is esti- 
mated. However, an exhaust gas sensor may be disposed in 
the vicinity of the 0 2 sensor 8, and the exhaust gas tem- 
perature Tgd may be detected by the exhaust gas sensor. 
Even if the exhaust gas sensor is disposed in a location 
remote from the 0 2 sensor 8, if the temperature detected by 
the exhaust gas sensor is substantially the same as the ex- 
haust gas temperature Tgd near the 0 2 sensor 8 due to the 
layout of the exhaust system or the like, then the detected 
temperature may be substituted for the exhaust gas tempera- 
ture Tgd near the 0 2 sensor 8. If the temperature detected 
by the exhaust gas sensor is used, then the temperature 
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(the latest value) detected by the exhaust gas sensor is 
used as the value of the exhaust gas temperature Tgd in the 
equation (10-1) to estimate the element temperature T 0 2 and 
the heater temperature Tht . Furthermore, the temperature 
(the latest value) detected by the exhaust gas sensor is 
used as the value of the exhaust gas temperature Tgd in the 
equation (24), the equation (28) , the equation (36), or the 
equation (40) to calculate the basic duty cycle SDUT and 
then calculate the duty cycle DUT according to the equation 
(25) . 

If the exhaust gas sensor is disposed in a loca- 
tion remote from the 0 2 sensor 8 and if the temperature de- 
tected by the exhaust gas sensor is not necessarily be same 
as the exhaust gas temperature Tgd near the 0 2 sensor 8, 
then it is possible to estimate the exhaust gas temperature 
Tgd near the 0 2 sensor 8 using the temperature detected by 
the exhaust gas sensor. For example, if an exhaust gas 
temperature sensor is provided which is capable of detect- 
ing the exhaust gas temperature Tgd in the partial exhaust 
passageway 3b, then the exhaust gas temperature Tgd near 
the 0 2 sensor 8 can be estimated by the calculations of the 
equations (8), (9) by using the detected temperature as the 
exhaust gas temperature Tgb in the equation (8-1) . In this 
case, the calculations of the equations (1), (5), (6) are 
not required. 
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In the first and second embodiments, both the 
element temperature T02 and the heater temperature Tht are 
estimated. However, either one of them may be detected di- 
rectly by a temperature sensor. If the element temperature 
T 0 2 is detected, then the heater temperature Tht may be es- 
timated using the detected value of the element temperature 
T 0 2 as the value of the element temperature T02 in the equa- 
tion (10-2) . The duty cycle DUT may be calculated using 
the estimated value as the heater temperature Tht in the 
calculation of the equation (24), the equation (28), the 
equation (36) , or the equation (40) , and also using the de- 
tected value of the element temperature T Q 2 as the value of 
the element temperature T 0 2- If the heater temperature Tht 
is detected, then the element temperature T Q 2 may be esti- 
mated using the detected value of the heater temperature 
Tht as the value of the heater temperature Tht in the equa- 
tion (10-1) . The duty cycle DUT may be calculated using 
the estimated value as the element temperature T 0 2 in the 
calculation of the equation (24), the equation (28), the 
equation (36), or the equation (40), and also using the de- 
tected value of the heater temperature Tht as the value of 
the heater temperature Tht. If both the element tempera- 
ture T 0 2 and the heater temperature Tht are detected by the 
temperature sensors, then the detected temperatures may be 
used as the values of the element temperature T Q 2 and the 
heater temperature Tht in the calculation of the equation 
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(24), the equation (28) f the equation (36), or the equation 
(40) to calculate the duty cycle DUT. 

In the first and second embodiments, the element 
temperature T 0 2 of the 0 2 sensor 8 or the heater temperature 
Tht is controlled at the target value R or R' according to 
the optimum predictive control algorithm. However, the 
present invention is not limited to the optimum predictive 
control algorithm. 

For example, the control input DUT may be deter- 
mined according to an ordinary optimum control algorithm 
which includes no predictive control algorithm. In this 
case, the control input DUT may sequentially be calculated 
according to an equation which is produced by removing the 
fourth term (the term including R(n+i)) and the fifth term 
(the term including Tgd(n)) from the right side of the 
equation (24) or the equation (36) or by removing the 
fourth term (the term including R' (n+i) ) and the fifth term 
(the term including Tgd(n)) from the right side of the 
equation (28) or the equation (40) . According to this 
modification, the heater controller for determining the 
control input DUT is an optimum servocontroller for deter- 
mining the control input DUT in order to minimize the value 
of the evaluating function JO or Jl or J2 or J3 where M = 0 
in the equation (17), the equation (27), the equation (35) 
or the equation (39) . 
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Alternatively, the control Input DUT may calcu- 
lated according to an equation which is produced by remov- 
ing one or two of the third term (the component depending 
on the heater temperature Tht) # the fourth term (the term 
including R(n+i)), and the fifth term (the term including 
Tgd(n)) from the right side of the equation (24). Further 
alternatively, the control input DUT may calculated accord- 
ing to an equation which is produced by removing one or two 
of the third term (the component depending on the element 
temperature T02), the fourth term (the term including 
R'(n+i)), and the fifth term (the term including Tgd(n)) 
from the right side of the equation (28). Further alterna- 
tively, the control input DUT may calculated according to 
an equation which is produced by removing one or two of the 
third term (the component depending on the heater tempera- 
ture Tht), the fourth term (the term including R(n+i)), and 
the fifth term (the term including Tgd(n) ) from the right 
side of the equation (36). Further alternatively, the con- 
trol input DUT may calculated according to an equation 
which is produced by removing one or two of the second term 
(the component depending on the element temperature T 02 ) . 
the fourth term (the term including R'(n+i)), and the fifth 
term (the term including Tgd(n)) from the right side of the 
equation (40). The component depending on the element tem- 
perature difference e in the equation (24) or the equation 
(36) and the component depending on the heater temperature 
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difference e f in the equation (28) or the equation (40) may 
be determined according to a PI control law or a PID con- 
trol law. 

In each of the embodiments described above, the 
element temperature T 0 2 of the 0 2 sensor 8 is controlled. 
However, the present invention is also applicable to an ex- 
haust gas sensor other than the 0 2 sensor 8 (e.g., the 
wide-range air-fuel ratio sensor 9 or a humidity sensor for 
generating an output signal representative of the water 
content of the exhaust gas) . 

The internal combustion engine to which the pre- 
sent invention is applicable may be an ordinary port- 
injected internal combustion engine, a spark-ignition in- 
ternal combustion engine with direct fuel injection into 
cylinders, a diesel engine, an internal combustion engine 
for use as an outboard engine on a boat, etc. 

Industrial applicability: 

As described above, the present invention is 
useful as a technology for appropriately controlling the 
temperature of an exhaust gas sensor disposed in the ex- 
haust system of an internal combustion engine mounted on an 
automobile, a hybrid vehicle, an outboard engine assembly, 
or the like, at a desirable temperature for stabilizing the 
output characteristics of the exhaust gas sensor. 



